M=[4 -2 1; 1089 33 1; 1 1 1] ,b=[-1500;-900;1321] ,sol=inv(M)*b; a2=sol(1),a1=sol(2),a0=sol(3), f=sprintf("[%f.*x.^2+%f.*x+%f, 0.*x]",a2,a1,a0); fplot(f,[-100, 100]), grid on,titulo=sprintf(" ECUACION CUADRATICA %f*x^2+%f*x+%f",a2,a1,a0),title(titulo), print -djpeg grafica.jpg;
