?? naeuler2s.m
字號:
function [x,y]=naeuler2s(dyfun,xspan,y0,h)
%用途: 2階改進euler法解常微分方程組y'=f(x,y),y(x0)=y0
%格式: [x,y]=naeuler2s(dyfun,xspan,y0,h),dyfun為函數f(x,y),xspan為求解區間[x0,xN]
% y0為初值向量y(x0),h為步長,x返回節點,y返回數值解
x=xspan(1):h:xspan(2);
y=zeros(length(y0),length(x));
y(:,1)=y0(:);
for n=1:length(x)-1
k1=feval(dyfun,x(n),y(:,n));
y(:,n+1)=y(:,n)+h*k1;
k2=feval(dyfun,x(n+1),y(:,n+1));
y(:,n+1)=y(:,n)+h*(k1+k2)/2;
end
x=x';y=y';
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -