
h = (b-a)/N;
x = zeros(1,N+1);
y = zeros(1,N+1);
x(1) = a;
y(1) = y0;
for n=1:N
x(n+1)=x(n)+h;
y(n+1)=y(n)+h*feval(f,x(n),y(n))
end
fun.m
function f = fun(x,y)
%FUN 此处显示有关此函数的摘要
% 此处显示详细说明
f = 8-3*y;
end