函数文件1: function b=F(x0,h,u,tau) b(,)=x0()-u(); b(,)=x0()-u()+*h*1e8*cos(tau)*x0(); 函数文件2: function g=Jacobian(x0,h,tau) g(,)=; g(,)=; g(,)=; g(,)=+*h*1e8*cos(tau); 函数文件3: function x=Euler(h,x0,u,tau) % x0 表示迭代初值 % u 表示上一节点数值 tol=1e-; x1=x0-Jacobian(x…