📄 demodifeuler.m
字号:
function y = DEModifEuler(f, h,a,b,y0,varvec)
format long;
N = (b-a)/h;
y = zeros(N+1,1);
y(1) = y0;
x = a:h:b;
var = findsym(f);
for i=2:N+1
v1 = Funval(f,varvec,[x(i-1) y(i-1)]);
t = y(i-1) + h*v1;
v2 = Funval(f,varvec,[x(i) t]);
y(i) = y(i-1)+h*(v1+v2)/2;
end
format short;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -