作者：夏云木子 1、（1） >> clear
>> y=dsolve('D2y=y','y(0)=0,Dy(0)=1'); >> disp('y(t)='),pretty(y)
（2） >> clear
>> y=dsolve('D2y=1+(1+x^2)*y','y(0)=1,Dy(0)=3','x'); >> disp('y(x)='),pretty(y)
（3）、 >> clear
>> Ti_713=@(x,Y)[Y(2);-Y(1)-(Y(1)^2-1)*Y(2)]; >> ode23(Ti_713,[0,5],[0,0.25]),grid
(4)、 >> clear
>> y=dsolve('D3y=-y','y(0)=1,Dy(0)=D2y(0)=0'); >> disp('y(t)='),pretty(y)
（5）、 >> clear
>> Ti_715=@(x,Y)[Y(2);Y(3);2*Y(3)+3*Y(2)+3*exp(2*x)]; >> ode23(Ti_715,[0,0.15],[1,10,30]),grid
2、（1） >> clear
>> y=dsolve('Dy=2*x-y','y(1)=1','x')
>> syms x
>> ezplot(2*x+exp(-x)*exp(1)-2,[1,2]),grid
（2）、
>> clear
>> y=dsolve('Dy-y-x/y','y(0)=1','x')
>> syms x
>> ezplot(((3*exp(2*x))/2-x-1/2)^(1/2),[1,2]),grid
（3）、 >> clear
>> y=dsolve('Dy=-x*y^2','y(0)=2','x')
(4)、 >> clear
>> y=dsolve('Dy=3*y/(1+x)-y','y(0)=1','x')
>> ezplot(y,[0,5]),grid
（5）、
>> clear
>> y=dsolve('x^2+y+(x-2*y)*Dy=0','y(0)=1','x')
>> ezplot(y,[0,5]),grid
3、（1） >> clear
>> [x,y]=dsolve('Dx=y,Dy=-x')
（2）、 >> clear
>> [f,g]=dsolve('Df=3*f+4*g,Dg=-4*f+3*g','f(0)=0,g(0)=1','x')
计算方法及其MATLAB实现第七章作业
