[求助-----help]
可以用什么方法求解啊？只要一组解就可以了。周一要交作业啦，救命！谢谢！方程组：
b1+b2+b3+b4=1
b1*c1+b2*c2+b3*c3+b4*c4=1/2
b1*c1^2+b2*c2^2+b3*c3^2+b4*c4^2=1/3
b1*c1^3+b2*c2^3+b3*c3^3+b4*c4^3=1/4
b3*a32*c2+b4*a42*c2+b4*a43*c3=1/6
b3*a32*c2*c3+b4*a42*c2*c4+b4*a43*c3*c4=1/8
b3*a32*c2^2+b4*a42*c3^2+b4*a43*c4^2=1/12
b4*a43*a32*c2^2=1/24
其中b1,b2,b3,b4,c1,c2,c3,c4,a32,a42,a43均为未知数，且都为四阶的显式Runge-Kutta方法中的系数
四阶的显式Runge-Kutta方法为：
function E=rk4(a,b,ya,m)
% a,b区间两端点值
% ya为给定某一点的函数值
% m最大迭代次数
h=(b-a)/m;
T=zeros(1,m-1);
Y=zeros(1,m-1);
T=a:h:b;
Y(1)=ya;
for y=1:m-1;
k1=f(T(j),Y(j));
k2=f(T(j)+c2*h,Y(j)+h*a21*k1);
k3=f(T(j)+c3*h,Y(j)+h*(a31*k2+a32*k2));
k4=f(T(j)+c4*h,Y(j)+h*(a43*k3+a42*k2+a41*k1));
Y(j+1)=Y(j)+h*(b1*k1+b2*k2+b3*k3+b4*k4);
end
E=;
页:
[1]