文件说明
1、使用说明
1)加载f.m文件至matlab
2)首先在命令行使用solve函数解得方程x^5=-1在复平面上的5个根;
>> solve('x^5=-1')
-1
1/4 + (2^(1/2)*(5^(1/2) + 5)^(1/2)*i)/4 - 5^(1/2)/4
1/4 - (2^(1/2)*(5^(1/2) + 5)^(1/2)*i)/4 - 5^(1/2)/4
(2^(1/2)*(5 - 5^(1/2))^(1/2)*i)/4 + 5^(1/2)/4 + 1/4
5^(1/2)/4 - (2^(1/2)*(5 - 5^(1/2))^(1/2)*i)/4 + 1/4
3)然后将这5个根代入f.m文件中函数f(A,B,C,D,E)中;
>> f(1/4 + (2^(1/2)*(5^(1/2) + 5)^(1/2)*i)/4 - 5^(1/2)/4, 1/4 - (2^(1/2)*(5^(1/2) +
5)^(1/2)*i)/4 - 5^(1/2)/4,(2^(1/2)*(5 - 5^(1/2))^(1/2)*i)/4 + 5^(1/2)/4 + 1/4, 5^(1/2)/4 -
(2^(1/2)*(5 - 5^(1/2))^(1/2)*i)/4 + 1/4, -1)
4)等待matlab运行完毕即可得到如下牛顿初始点收敛图
2、代码说明
function f(A,B,C,D,E)