implicit real*8(a-h,o-z)
      dimension u(0:4,0:200),t(0:200)
      real m1,m2,l1,l2,g
      real h,a1,a2,a3,alpha,betha,gama
      integer n
      n=50
      h=0.1
      l1=0.3
      l2=0.3
      g=9.8
      m1=0.1
      m2=0.1
      alpha=(m1*l1*l2)/(g*(m1+m2))
      betha=(m1*l1*(m2/m1+1.))/(m1+m2)+l2
      gama=(m2/m1)*g
      a1=12./(12.*alpha+h**2.*betha)
      a2=betha/(12.*alpha/h**2.+betha)
      a3=gama/(12.*alpha/h**4.+betha/h**2.)
c      print *, n,h,l1,l2,g,alpha,betha,gama,a1,a2,a3
      t(0)=0.
      do i=1,n
      t(i)=t(i-1)+0.05
      enddo
      do 1 il=0,0
      sl=il
      u(il,0)=5.
      u(il,1)=-4.5
      do 2 i=1,n-1
      u(il,i+1)=a1*(2.*u(il,i)-2.*u(il,i)+u(il,i+1))
     & +a2*(2.*u(il,i)-u(il,i+1)-u(il,i-1))-a3*u(il,i) 

      print 28, t(i),u(il,i+1)
  28  format(f8.2,e15.4)
   2  continue
   1  continue
      end