利用dsolve（）函数，可求得常微分方程的初值问题 （1+x^2)y''=2xy'的解析解。
实现代码
syms y(x),D2y=diff(y,2);Dy=diff(y,1);
disp('常微分方程的解析解')
y=dsolve((1+x^2)*D2y==2*x*Dy,y(0)==1,Dy(0)==3)
