Question

怎么用MATLAB写龙格库塔四阶方程

Answer 1

四阶龙格-库塔是求解常微分方程（常微分方程组）精度最高的一种数值方法。四阶龙格-库塔迭代公式为

根据四阶龙格-库塔迭代公式，可用matlab软件编写得到如下计算代码：

调用方法：

clc;close all;

y0=1;

h=0.1;a=0;b=1;

[x,y] = runge_kutta(@func,y0,h,a,b);

disp('        x        y')

[x',y']

x=0:0.1:1; 

plot(x,y,':'),grid on

xlabel('x'),ylabel('y（x）/y’(x)');

title('函数图');

function f = func(x,y)

f=sqrt(x+y); 

end