f=inline('cos(x)+sin(y)','x','y'); %微分方程的右边项 dx=0.05; %x方向步长 xleft=pi/2; %区域的左边界 xright=3*pi/2; %区域的右边界 xx=xleft:dx:xright; %一系列离散的点 n=length(xx); %点的个数 y0=0; %%(1)欧拉法
