如何用Runge-Kutta迭代求解二阶常微分方程组

http://wapwenku.baidu.com/view/667bfc264431b90d6c85c79c.html?ssid=0&from=1002253n&uid=0&pu=sz%401320_480%2Ccuid%400O-jiY8uHtgjaS8m_P-j8la2SaYJa-aNlu-1ilax2a87a2tggPvC8_u6v8jxa2fHA%2Ccua%40_a-qi4uq-igBNE6lI5me6NNy2I_UCvh4SdNqA%2Ccut%405tVgirktSh_Uh2I4gIvwiyNqmoi5pQqAC%2Cosname%40baiduboxapp%2Cctv%402%2Ccfrom%401002482e%2Ccen%40cuid_cua_cut%2Ccsrc%40searchf_box_txt%2Cta%40zbios_2_6.0_6_7.6%2Cusm%400%2Cvmgdb%400020100228y&bd_page_type=1&baiduid=EE189F6EFEBB76B49175343889EE42BF&tj=wenkuala_3_0_10_l1
这个PPT文档的第5页

您好,希望以下回答能帮助您
一.用matlab 中的solve函数
>>syms x y; %定义两个符号变量；
>>[x ,y]=solve('y=2*x+3','y=3*x-7');%定义一个 2x1 的数组,存放x,y
>>x
>>x=10.0000
>>y
>>y=23.0000
二.用matlab 中的反向斜线运算符（backward slash）
分析：
方程组可化为
2*x-y=-3；
3*x-y=7；
AX=B (*)
A=[2,-1;3,-1]; B=[-3,7];
X=A\B %可以看成将（*）式左边都除以系数矩阵A
>>A=[2,-1;3,-1];
>>B=[-3,7];
>>X=A\b
X =
10.0000 % x = 10.0000
23.0000 % y = 23.0000