x=var('x')
t=var('t')
c=var('c') # An arbitary Constant
u=function('u',x,t)
u=sin(9*t)*sin((1/4)*x)
pde=diff(u,t,t)-c^2*diff(u,x,x)
print solve(pde,c)
#Thus, given 'u' is a solution of PDE since we can find a suitable constant 'c' for which equation is satisfied.
plot3d(lambda x, t: sin(9*t)*sin((1/4)*x), (-2,2), (-2,2))

SageDays/Kreyszig-11.1-3 (last edited 2010-08-12 09:33:42 by Aashita)