一、四阶龙格库塔法
u
四阶龙格库塔原理:
u
y
(i+1)=y(
i
)+k*h //h
为步长,
K
为平均系数
u
a=0;b=20;n=1000;h=(b-a)/
n;T
=
a:h:b;X
=zeros(1,n+1);
u
Y=zeros(1,n+1);f=@(
t,x,y
)(1*y);g=@(
t,x,y
)(-2*x-3*y+1);X(1)=1;Y(1)=0;
u
for i=1:n;f1=f(T(i),X(i),Y(i)); g1=g(T(i),X(i),Y(i));
u
f2=f(T(
i
)+h/2,X(
i
)+f1*h/2,Y(
i
)+g1*h/2);
u
g2=g(T(
i
)+h/2,X(
i
)+f1*h/2,Y(
i
)+g1*h/2);
u
f3=f(T(
i
)+h/2,X(
i
)+f2*h/2,Y(
i
)+g2*h/2);
u