2Â¥: Originally posted by hsdwlxhx at 2015-08-10 20:09:31
5άµü´ú£¬ÄãҲ̫bugÁË

3Â¥: Originally posted by RK1224 at 2015-08-10 20:43:33
ÎÒÊDzËÄñ¡£¡£¡£Ï£ÍûÄãÖ¸µãÒ»ÏÂ...

5Â¥: Originally posted by dingd at 2015-08-10 23:28:57
1£ºµÚÒ»¸ö¹«Ê½×óÓÒÀ¨ºÅ²»Æ¥Å䣻
2£ºµÚ2¡¢µÚ3¹«Ê½Öеġ°e¡±ÊÇʲô£¿Ó¦¸ÃÊÇ¡°e+0¡±°É£¡

6Â¥: Originally posted by RK1224 at 2015-08-11 13:32:07
ллÄãµÄ»Ø¸´£¬È·ÊµÊÇÎÒûдÇå³þ£¬ÎÒ¸ÄÁËÒ»ÏÂÈçÏ£º
(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296=0;

38708*x1-38943*x2+1174.2=0;

-23679*x2+47828*x3-24149*x4=0;

(9.2333e-08)*x ...

6Â¥: Originally posted by RK1224 at 2015-08-11 13:32:07
ллÄãµÄ»Ø¸´£¬È·ÊµÊÇÎÒûдÇå³þ£¬ÎÒ¸ÄÁËÒ»ÏÂÈçÏ£º
(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296=0;

38708*x1-38943*x2+1174.2=0;

-23679*x2+47828*x3-24149*x4=0;

(9.2333e-08)*x ...

function solve_eq
clear all;clc
format long

x0=[0 0 0 0 0];

options=optimset('TolFun',1e-9);
[x fval]=fsolve(@fun,x0,options);
[x;fval]'

function f=fun(u)



x1=u(1);
x2=u(2);
x3=u(3);
x4=u(4);
x5=u(5);

f(1)=(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296;

f(2)=38708*x1-38943*x2+1174.2;

f(3)=-23679*x2+47828*x3-24149*x4;

f(4)=(9.2333e-08)*x1^4+23914*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8;

f(5)=5.3514*x3-83.6514*x5+1566;