        DIMENSION XMU(30),DI(30),BX(30),BX1(30),BX2(30),BX3(30),BX4(30)
        DIMENSION XIP(30),XIPP(30),EXPP(30),X(30)
        DX=315.0
        E=3.0E+07
        N=1
        L=6300.0
        W=400.0
        RHOW=9.355E-05
        BETN=(2*N+1.0)*1.5708/L
        S1=0
        S2=0
        S3=0
        DI(0)=5.0E+07
        DI(1)=5.0E+07
        DI(2)=8.3E+07
        DI(3)=10.0E+07
        DI(4)=12.4E+07
        DI(5)=13.0E+07
        DI(6)=14.0E+07
        DI(7)=15.0E+07
        DI(8)=16.0E+07
        DI(9)=16.4E+07
        DI(10)=16.6E+07
        DI(11)=16.0E+07
        DI(12)=15.5E+07
        DI(13)=14.5E+07
        DI(14)=13.5E+07
        DI(15)=12.5E+07
        DI(16)=10.5E+07
        DI(17)=8.3E+07
        DI(18)=7.5E+07
        DI(19)=6.0E+07
        DI(20)=5.8E+07
        XMU(0)=6.0
        XMU(1)=6.0
        XMU(2)=8.6
        XMU(3)=13.0
        XMU(4)=17.2
        XMU(5)=21.0
        XMU(6)=16.0
        XMU(7)=16.0
        XMU(8)=18.0
        XMU(9)=21.0
        XMU(10)=25.0
        XMU(11)=17.0
        XMU(12)=17.0
        XMU(13)=16.8
        XMU(14)=12.0
        XMU(15)=11.0
        XMU(16)=8.6
        XMU(17)=11.6
        XMU(18)=11.6
        XMU(19)=8.0
        XMU(20)=4.0
        X(0)=0 
        X(1)=315.0
        X(2)=630.0
        X(3)=945.0
        X(4)=1260.0
        X(5)=1575.0
        X(6)=1890.0
        X(7)=2205.0
        X(8)=2520.0
        X(9)=2835.0
        X(10)=3150.0
        X(11)=3465.0
        X(12)=3780.0
        X(13)=4095.0
        X(14)=4410.0
        X(15)=4725.0
        X(16)=5040.0
        X(17)=5355.0
        X(18)=5670.0
        X(19)=5985.0
        X(20)=6300.0
        M1=2
        M2=19
        M3=1 
        DO 10 M=M1,M2,M3 
        A=BETN*X(M)               
        BX(M)=EXP(-A)+COS(A)-SIN(A)
C        PRINT *,BX(M)
        BX1(M)=-BETN*(EXP(-A)+SIN(A)+COS(A))
        BX2(M)=(BETN**2)*(EXP(-A)-COS(A)+SIN(A))
        BX3(M)=(BETN**3)*(-EXP(-A)+SIN(A)+COS(A))
        BX4(M)=(BETN**4)*BX(M)
        XIP(M)=(DI(M+1)-DI(M-1))/(2.0*DX)
        XIPP(M)=(DI(M+1)+DI(M-1)-2.0*DI(M))/(DX**2)
        Q=E*DI(M)*BX4(M)
        R=2.0*E*XIP(M)*BX3(M)
        S=E*XIPP(M)*BX2(M)
        EXPP(M)=Q+R+S
C        PRINT *,S
        S1=S1+BX(M)*EXPP(M)*DX
C        PRINT *,S1
        S2=S2+(BX(M)**2)*XMU(M)*DX
C        PRINT *,S2
        S3=L
10      CONTINUE
        S3=L
        T=S3*0.5*W*RHOW*SQRT(W*L)
        PRINT *,S1,S2,T
        OMEGA =SQRT(S1/(S2+T))
        PRINT *,N,OMEGA
        END
           
  

