Question

ind 2 imaginary roots of equation 1/(x*3+xA2+1), mullers method, matlab nswer only if sure

Answer 1

Hi,

The code to your question is

clc
e=0.0001;
x1=-1;x2=0;x3=1;
f=inline('1/((x^3)+x^2+1)','x');
N=0;
while (1)
if abs(f(x3))>e
c=f(x3);
t=[((x1-x3)^2),(x1-x3);((x2-x3)^2),(x2-x3)];
p=[f(x1)-f(x3);f(x2)-f(x3)];
L=inv(t)*p;
a=L(1);
b=L(2);
if b>0
x=x3-((2*c)/(b+sqrt(b^2-4*a*c)));
else
x=x3-((2*c)/(b-sqrt(b^2-4*a*c)));
end
x1=x2;
x2=x3;
x3=x;

else
break;
end
end
x3
conj(x3)

Revert in case of any queries

Thanks.