Tengo un pequeño problema, tengo un polinomio de la forma:
p(x) = a(1) + (x+a(2)) + (x+a(3))^2 + (x+a(4))^3 + (x+a(5))^4 + (x+a(6))^5);
y un vector fila a=[a(1) a(2) a(3) a(4) a(5) a(6)]
Quiero obtener un fichero de funcion que haga esto en la ventana de comandos:
>> a=[1, 2, 3, 4, 5,6]; -->Introdusco un vector a de 6 elementos
>> raices=prueba2(a) ---> llamo a la funcion prueba2 para obtener las raices del polinomio.
Luego que me arroje las raices del polinomio p(x) evaluado en el vector "a" . He tenido problemas por la forma del polinomio. Agradeceria alguna ayuda con los codigos que debo utilizar. Gracias.
Hasta ahora he intentado esto:
function raices=pregunta2(a)
syms x
if nargin >=6
error('Demasiados argumentos de entrada')
end
a=input();
syms x
p=a(1)+(x+a(2))+(x+a(3))^2 +(x+a(4))^3 +(x+a(5))^4+(x+a(6))^5;
m=expand(p);
n=coeffs(m);
w=sort(n);
raices=roots(w);
dip(raices)