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...