2021-01-07Matlab数値解析通常の微分方程式初期境界値問題数値解法標準RungeKutta4次4段階式オイラー法

Matlab数値解析通常微分方程式初期境界値問題数値解法標準RungeKuta4次4段式オイラー法

1.標準ルンゲクタ4次4段式

function y=rk4(fun,a,b,y0,n)
h=(b-a)/n;
y(1)=y0;
for k=1:n
  x=a+(k-1)*h;
  k1=h*feval(fun,x,y(k));
  k2=h*feval(fun,x+h/2,y(k)+k1/2);
  k3=h*feval(fun,x+h/2,y(k)+k2/2);
  k4=h*feval(fun,x+h,y(k)+k3);
  y(k+1)=y(k)+(k1+2*k2+2*k3+k4)/6;
end

関数ファイル

 

function u=frk4(x,y)u=y-2*x/y;

呼び出しプログラム

y=rk4(@frk4,0,1,1,10)

2オイラー法

function y=euler(f,a,b,y0,h)
n=(b-a)/h;
y(1)=y0;
for i=1:n
  x(i)=a+(i-1)*h;
  y(i+1)=h*feval(f,x(i),y(i));
end

関数ファイル

 

function u=feuler(x,y)u=x^3-y/x;

呼び出しプログラム

y=euler(@feuler,1,2,0.4,0.2)