n85chapt9ex5euler.m

来自「MATLAB 说明您提供的文件包含的内容和具体功能(至少要20个」· M 代码 · 共 20 行

function E=N85chapt9ex5Euler(a,b,N,ya)
%学号 20042374  姓名 钟健升   作业号  85
%p264 第5题Euler方法
%f是原方程化为微分方程后的右端表达式
%g是微分方程的解的表达式
%a,b是区间[a,b]中的a,b
%N是等分区间的个数
syms x y 
f=exp(x^2);
h=(b-a)/N;
y=zeros(1,N+1);
x=zeros(1,N+1);
y(1)=ya;
x=a:h:b;
for i=1:N
    y(i+1)=y(i)+h*exp(x(i)^2);
end 
E=[x',y'];