--Copyright The Numerical Algorithms Group Limited 1994.
)clear all
y := operator y
deq := differentiate(y x, x, 2) + 2*w[0]*differentiate(y x, x) + w[0]**2*y x
sol:= solve(deq = sin (w*x), y, x=0,[0,0])
work:= sol *sin(w*x)
integrate(work,x)
D(%,x)-work
simplify %

