Initial condition

> alpha[1,3]:=BesselJZeros(1,3);

> u:=(r,theta,t)->BesselJ(1,alpha[1,3]*r)*cos(theta)*cos(alpha[1,3]*t);

>

Plot of Initial condition

Animation of Membrane