?? delgkt4_qt.m
字號(hào):
function y = DELGKT4_qt(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
K1 = Funval(f,varvec,[x(i-1) y(i-1)]);
K2 = Funval(f,varvec,[x(i-1)+h/3 y(i-1)+K1*h/3]);
K3 = Funval(f,varvec,[x(i-1)+2*h/3 y(i-1)-K1*h/3+K2*h]);
K4 = Funval(f,varvec,[x(i-1)+h y(i-1)+h*(K1-K2+K3)]);
y(i) = y(i-1)+h*(K1+3*K2+3*K3+K4)/8;
end
format short;
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -