double precision pi_ common /mathcons/pi_ double precision L common /plankcons/L double precision ksol,f1R,f1I,f3,f4 common /planksol/ksol,f1R,f1I,f3,f4 integer Xpoints parameter (Xpoints = 10000) integer Modes parameter (Modes = 50) double precision Norm(Modes) double precision Coef(Modes,Modes) double precision Unit(Modes,Modes) double precision alphaR(Modes),alphaI(Modes),beta(Modes) double precision Vl(Modes,Modes),Vr(Modes,Modes) double precision kmode(Modes) double precision f1Rmode(Modes),f1Imode(Modes), + f3mode(Modes),f4mode(Modes) common /plankmat/Norm,Coef,Unit,alphaR,alphaI,beta,Vl,Vr, + kmode,f1Rmode,f1Imode,f3mode,f4mode real evalue(Modes) integer iorder(Modes) common /plankmodes/iorder,evalue