?? feuler.m
字號:
%改進的euler法解一維微分方程的通用子函數
%函數的輸入變量為:T為x的取值上限,n為把[0,T]區間等份的數目,u0為方程的初始值
function u=feuler(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
for i=1:n-1
u(i+1)=u(i)+h*f(h*i,u(i));
u(i+1)=u(i)+(h/2)*(f(h*(i+1),u(i+1))+f(h*i,u(i)));
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -