x1=input("Dame el valor de x1: ")
b1=input("Dame el Valor del P(x) =: ");
printf("P(%f) = a0 + %f a1 + %f a2 = %f\n",x1,x1,x1^2,b1)

x2=input("Dame el valor de x2: ")
b2=input("Dame el Valor del P(x) =: ");
printf("P(%f) = a0 + %f a1 + %f a2 = %f\n",x2,x2,x2^2,b2)

x3=input("Dame el valor de x3: ")
b3=input("Dame el valor del P(x) =: ");
printf("P(%f) = a0 + %f a1 + %f a2 = %f\n",x3,x3,x3^2,b3)

M=[1 x1 x1^2; 1 x2 x2^2; 1 x3 x3^2]
b=[b1;b2;b3]

inv(M)

incognitas=inv(M)*b

a0=incognitas(1);
a1=incognitas(2);
a2=incognitas(3);

x=[-10:.2:10];
px=a0+a1.*x+a2.*x.^2;

plot(x,px,'r*-');
grid on
title('Ecuación cuadrática P(x)')
xlabel("Valores de X")
ylabel("P(x)")

print -djpeg "grafica.jpg"




