献丑了。
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function x=mysolve3(a,b)
%%求解x^3+a*x^2-b=0的解
p=(-a*a)/3;
q=(2*a^3+27*(-b))/27;
D=(p/3)^3+(q/2)^2;
if (D>=0)
temp=-q/2+sqrt(D);
if (temp>0) u=(temp)^(1/3); else u=-(-temp)^(1/3); end
temp=-q/2-sqrt(D);
if (temp>0) v=(temp)^(1/3); else v=-(-temp)^(1/3); end
y(1)=u+v;
y(2)=-(u+v)/2+i*sqrt(3)*(u-v)/2;
y(3)=-(u+v)/2-i*sqrt(3)*(u-v)/2;
else
temp=acos(-q/2*(-p/3)^-1.5);
y(1)=2*sqrt(-p/3)*cos(temp/3);
y(2)=-2*sqrt(-p/3)*cos((temp+pi)/3);
y(3)=-2*sqrt(-p/3)*cos((temp-pi)/3);
end
x=y-a/3;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
mysolve3(1,3)
得到三个解
你可以通过vpa(solve('x^3+x^2-3'),5)进行检验
温馨提示:内容为网友见解,仅供参考