Question: Why do I get different solutions for harmonic oscilators with forcing Dirac for standard and laplace methos in dsolve?


 

> 

restart

> 

eq := diff(y(t), t, t)+9*y(t) = Dirac(t)

diff(diff(y(t), t), t)+9*y(t) = Dirac(t)

(1)
> 

with(inttrans):

> 

icon := y(0) = 0, (D(y))(0) = 0;

y(0) = 0, (D(y))(0) = 0

(2)
> 

dsolve({eq, icon}, y(t), method = laplace);

y(t) = (1/3)*sin(3*t)

(3)
> 

dsolve({eq, icon}, y(t));

y(t) = (1/3)*(Heaviside(t)-1/2)*sin(3*t)

(4)

``


 

Download Dirac_for_Primes.mw

Please Wait...