matlab解多元高次方程组A B C D E为变量的5元次方程组,我一直调试不出结果,一运行就是busy,是不是程序哪里写错了?syms A B C D E;eq1='1/4*(110+220*A)*(3/25+B)*(69/50+A)*(11/100+C)-1377380*B*A*(11/400+5/4*C)*(1/2+A)-305/4*A*(11/400+5/4*C)*(3/25+B)*(69/50+A)=0';eq2='(271698404112853371/4611686018427387904
2019-04-12
matlab解多元高次方程组
A B C D E为变量的5元次方程组,我一直调试不出结果,一运行就是busy,是不是程序哪里写错了?
syms A B C D E;
eq1='1/4*(110+220*A)*(3/25+B)*(69/50+A)*(11/100+C)-1377380*B*A*(11/400+5/4*C)*(1/2+A)-305/4*A*(11/400+5/4*C)*(3/25+B)*(69/50+A)=0';
eq2='(271698404112853371/46116860184273879040+33962300514106665/1152921504606846976*D+18411/31250*A)*(3/25+B)*(69/50+A)*(9/10+B)*((49/2+49*E)*(1/2+E)+81/625)-10901624/125*B*A*(6712548807494023/18446744073709551616+4214072214526101626067/360287970189639680000000*D+6498/390625*A)*(9/10+B)*((49/2+49*E)*(1/2+E)+81/625)-3715146/78125*B*(6712548807494023/18446744073709551616+4214072214526101626067/360287970189639680000000*D+6498/390625*A)*(3/25+B)*(69/50+A)=0';
eq3='4516*B*A*(3/20+C)+(21/1250+7/50*B)*(69/50+A)*(3/20+C)-9/50*C*(3/25+B)*(69/50+A)-1/500*C*(3/25+B)*(69/50+A)*(3/20+C)=0';
eq4='6156/625*B*(1/20+D)-17/200*D*(1/20+D)*(9/10+B)*((49+98*E)*(1+2*E)+81/625)+9/50*D*(9/10+B)*((49+98*E)*(1+2*E)+81/625)=0';
eq5='17/200*D-7/100*E=0';
[A,B,C,D,E]=solve(eq1,eq2,eq3,eq4,eq5,'A,B,C,D,E')
A B C D E为变量的5元6次方程组,用的是MATLAB7.0
问题是我busy之后没有结果啊!正常的话你把结果发上来!
优质解答
用Newton-Raphson法试了下：
给定初值[1,1,1,1,1],得到
-0.43992713717969
-0.00010699797645
-0.16462315624105
0.00000678235923
0.00000823572192
给定初值[1,2,3,4,5],得到另一个
-0.01494259126252
-0.00029926103702
-0.46941875159537
2.06762484056609
2.51068730640168
用Newton-Raphson法试了下：
给定初值[1,1,1,1,1],得到
-0.43992713717969
-0.00010699797645
-0.16462315624105
0.00000678235923
0.00000823572192
给定初值[1,2,3,4,5],得到另一个
-0.01494259126252
-0.00029926103702
-0.46941875159537
2.06762484056609
2.51068730640168