（2） 设方程式为y’=-40y,y(0)=2用欧拉数值方法与ode45求方程式的数值解 请大神用MATLAB编写啊

欧拉方法
第一步：建立并保存名为 funfcn1.m 的 m 文件函数；
function f=funfcn1(x,y)
f=-40*y;
第二步：建立并保存名为Qeuler1.m的m文件函数；
function[X,Y,n]=Qeuler1(funfcn,x0,y0,b,h)
x=x0;n=(b-x0)/h;X=zeros(n,1);y=y0;Y=zeros(n,1);k=1;X(k)=x;Y(k)=y;
for k=2:n+1
fxy=feval(funfcn,x,y);
x=x+h;y=y+h*fxy;
X(k)=x;Y(k)=y';
end
plot(X,Y,'rp');grid,xlabel('自变量X'),ylabel('因变量Y')
title('用向前欧拉公式计算dy/dx=f(x,y),y(x0)=y0在[x0,b]上的数值解')
第三步：建立如下程序并保存名为jisuan1_1.m的m文件
clc,clear
x0=0;y0=3;b=2;h1=0.1;%赋初始值
[x1,Y1 ,n1]=QEuler1(@funfcn1,x0,y0,b,h1);%调用函数