Initial condition

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

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

>

Plot of Initial condition

Animation of Membrane