微分方程

limit(f,x,3,'left')
f=sin(4*x)/x; limit(f)
g=(x^2-3*x)/x;
ezplot(1/g)
hold on
plot(double(solve(g))*[1,1],[-2,2],'--');
hold off

for i=0:10:30
s=dsolve('Dx=-2*x+8')
f=subs(s,'C1',i);
ezplot(f,[-1,1])
hold on
end

g=dsolve('D2y-Dy-2*y-2*t=0','y(0)=0','Dy(0)=0');

s=dsolve('D2x+Dx=0,Dp=-p-x','x(0)=4,Dx(0)=0,Dp(0)=0');
s=dsolve('D2x=-x,Dp=-p-x','x(0)=4,Dx(0)=0,Dp(0)=0');

注意Dx和x分两边放

s=dsolve('D2x+2*Dx=-x,Dp=x','x(0)=1,Dx(0)=0,p(0)=0');
ezplot(s.x,s.p);