program tawami
        implicit real*8 (a-h, o-z)
c アクリル梁の寸法（cm）
c t:板厚, h:桁高, ell:長さ
         t=0.3d0
         h=3.d0
         write(*,*) '梁の長さ(cm）を入れて下さい'
         read(*,*) ell
c 断面二次モーメント：xi
         xi=t*h**3/12.d0
c ヤング率:e （kgf/cm^2)
         e=3.d4
        open(7,file='tawami.d')
c 荷重を1kgf ずつ増やした時の
c 片持ち梁の先端のたわみδ=PL^3/(3EI)を求める
        do 10 p=0.d0,100.d0,  1.d0
         d=p*ell**3/3.d0/e/xi
         write(7,100) d, p
10      continue
        close(7)
100     format(1p2d11.3)
        end