Initial condition

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

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

>

Plot of Initial condition

Animation of Membrane