求y（t）数值解并绘制其图像，与解析解对比精度。

YNhhU

微分方程y^'' (x)-(2-1/x) y^' (x)+(1-1/x)y(x)=x^2 e^(-5x)，满足边界条件y(1)=π，y(π)=1，求y（t）数值解并绘制其图像，与解析解对比精度。

Zedd

你这是带边界条件的微分方程，采用有限差分法求解数值解，采用dsolve求解解析解。

sharkN

典型的边值微分方程问题（BVP ODE），Matlab可用bvp4c命令试试，下面使用更简单的1stOpt求解：
3 m9 D$ T. Z  N! Q
" s) g- [8 f& p5 I8 {* z! MVariable x=[1,pi],y=[pi,1];
Plot x[x],y,y'[y2];
. t& G' T( r* H6 _; _ODEFunction y''=x^2*exp(-5*x)-(1-1/x)*y+(2-1/x)*y';