% 廢話不多說，直接復制就可以運行，文章底部有參考例子，不明白就看看

function [ yy ] = R_K_4( f , y0 , x0 , xn , hh )

% f 是inline function的句柄

% y0 是初始值

% [x0 xn]是計算范圍

% hh 是步長

%% ===輸入判斷===

if (4 == nargin) % == 當沒有給步長的時候使用默認步長==

h = 0.1 ;

elseif (5 == nargin)

h = hh ;

end

%% ===參數初始化(可以修改，可從網上找相關參數的值)，以下是通用參數===

c = [1/6 1/3 1/3 1/6] ;

r = [0.5 0.5 1] ;

u = [0.5 0 0.5 0 0 1] ;

%% ===

x = x0:h:xn ; % == 計算節點==

y = y0 ;

N = length(x) ; % == 節點數==

for k = 1:N-1 % == 以下就直接按書本公式寫的，應該容易理解== kk(1) = f( x(k) , y(k) ) ;

kk(2) = f( x(k)+r(1)*h , y(k) + u(1)*h*kk(1) ) ;

kk(3) = f( x(k)+r(2)*h , y(k) + u(2)*h*kk(1) + u(3)*h*kk(2) ) ;

kk(4) = f( x(k)+r(3)*h , y(k) + u(4)*h*kk(1) + u(5)*h*kk(2) + u(6)*h*kk(3) ) ;

y(k+1) = y(k) + h* ( c(1)*kk(1) + c(2)*kk(2) + c(3)*kk(3) + c(4)*kk(4) ) ;

end

總結

以上是生活随笔為你收集整理的四阶龙格库塔matlab计算例题,四阶龙格库塔法matlab实现的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。