please solve in matlab and include the phase portrait

a) Find all equilibrium points and study their stability using linearization. b) Using quadratic Lyapunov functions, estimate the regions of attraction of each asymptotically stable equilibrium point. Try to make the estimates as large as you can. c) Draw the phase portnait of the system to find the exact regions of attraction and compare them with your estimates. x1=(x1+x13)+2x2x2=2x1(x2+x23) x1=x1x13+x2x2=x13x2x1=x2x2=x1x2+2x2321x25x1=x2x2=2x1+3x2+2sat(3x2)