> |
> | with(DEtools):with(plots):NumericEventHandler(invalid_operation=exception): |
Warning, the name changecoords has been redefined
Block Spike
> | de:=2*diff(y(t),t$2)+diff(y(t),t)+2*y(t)=Heaviside(t-5)-Heaviside(t-20); |
> | plot(Heaviside(t-5)-Heaviside(t-20),t=0..30); |
> | sol:=solve(dsolve({de,y(0)=0,D(y)(0)=0}),y(t)): |
> | plot(sol,t=0..40); |
Increasing input
> | g(t):=(Heaviside(t-5)-Heaviside(t-10))*(t-5)/5+Heaviside(t-10); |
> | plot(g(t),t=0..20); |
> | de:=diff(y(t),t$2)+4*y(t)=g(t); |
> | sol:=solve(dsolve({de,y(0)=0,D(y)(0)=0}),y(t)); |
> | plot(sol,t=0..16); |
Diracspike
> | de:=2*diff(y(t),t$2)+diff(y(t),t)+2*y(t)=Dirac(t-5); |
> | sol:=solve(dsolve({de,y(0)=0,D(y)(0)=0}),y(t)); |
> | plot(sol,t=4..25); |
> |
> |