(a)

> SOLN := dsolve( { MODEL, IC }, y(t) );

> SOLN := subs( ln(-1+t^2) = ln(1-t^2) + I*Pi, SOLN );

>