!縺薙▲縺｡繧剃ｽｿ縺�
!�堤ｨｮ譚先侭縺ｮ縺帙ｓ譁ｭ陬懈ｭ｣菫よ焚繧定ｪｿ縺ｹ繧九�
!迚�戟縺｡譴√�繝�ぅ繝｢縺ｮ逅�ｫ門､繧貞�縺吶→縺�
implicit real*8(a-h, o-z)
!蜈･蜉帙☆繧九ヵ繧｡繧､繝ｫ縺ｮ蜷榊燕
open(8,file='tawami.dat')
!蜃ｺ蜉帙☆繧九ヵ繧｡繧､繝ｫ縺ｮ蜷榊燕
open (7,file='2syukekka.dat')
b1=0.840d0!隗呈攝縺ｮ蟷�
b2=0.009d0!驪ｼ譚ｿ縺ｮ蟷�
h1=0.0!隗呈攝縺ｮ鬮倥＆
h2=0.500d0!驪ｼ譚ｿ縺ｮ鬮倥＆
do i2=100,330,10
h1=h1+1.0d-2
h3=h2-h1*2.0d0!逵溘ｓ荳ｭ縺ｮ蜷域攸縺ｮ髟ｷ縺�
ell=10.00d0!繧ｹ繝代Φ
e=7.694d9!繝､繝ｳ繧ｰ邇�
g=e/15d0!縺帙ｓ譁ｭ蠑ｾ諤ｧ菫よ焚
poa2=0.30d0
e2=206.0d9
g2=e2/((1.0d0+poa2)*2.0d0)
pkata=1.000d3!闕ｷ驥�(迚�戟縺｡)
kakun=1.0d0!隗呈攝縺ｮ謨ｰ(險育ｮ励〒菴ｿ縺�°繧我ｸ蠢懷ｮ滓焚)
koun=2.0!d0驪ｼ譚ｿ縺ｮ謨ｰ(險育ｮ励〒菴ｿ縺�°繧我ｸ蠢懷ｮ滓焚)
dmenseki=b1*h1*kakun*2.0d0+b2*h2*koun!譁ｭ髱｢遨�
ea=b1*h1*kakun*2.0d0*e+b2*h2*koun*e2!謠帷ｮ玲妙髱｢遨�
!yz
yz=(e2*(b2*(h2-h1)*(h2/2-(h2-h1)/2)))/(e2*(b2*(h2-h1))+e*2*kakun/2*b1*h1)
!譛ｨ譚宣Κ蛻�妙髱｢2谺｡繝｢繝ｼ繝｡繝ｳ繝�
dnijimoku=b1*kakun*(h2**3)/12.0d0-b1*kakun*(h3**3)/12.0d0
!驪ｼ譚ｿ驛ｨ蛻�妙髱｢2谺｡繝｢繝ｼ繝｡繝ｳ繝�
dnijikouhan=b2*koun*(h2**3)/12.0d0
!譖ｲ縺貞央諤ｧ
ei=dnijimoku*e+dnijikouhan*e2
hakob=b1*kakun+b2!bank縺ｮb
hakoh=h2-h1!bank縺ｮh
hakom=hakob*h1/(hakoh*b2)!cowper縺ｮm
hakom2=hakom**2!bank縺ｮm**2
hakom3=hakom**3!bank縺ｮm**3
hakon=hakob/hakoh!bank縺ｮn
hakon2=hakon**2!bank縺ｮn**2
a=e2/e!bank縺ｮﾎｱ
a2=a**2!banka縺ｮﾎｱ**
a3=a**3!banka縺ｮﾎｱ*3
hakok=(20.0d0*(a+3.0d0*hakom)**2)&
&/((e/g)*(60.0d0*hakom2*hakon2+60.0d0*a*hakom*hakon2)&
&+(e/g2)*(180.0d0*hakom3+300.0d0*a*hakom2+144.0d0*a2*hakom+24.0d0*a3)&
&+poa2*(30.0d0*hakom2+6.0d0*a*hakom-4.0d0*a2))!bank縺ｮk*
kea=hakok*ea!bank縺ｮk
syotoukata=(pkata*ell**3)/(3.0d0*ei)!蛻晉ｭ画｢∫炊隲門､(迚�戟縺｡)
timokata=syotoukata+&
&(pkata*ell)/kea!繝�ぅ繝｢繧ｷ繧ｧ繝ｳ繧ｳ逅�ｫ門､(迚�戟縺｡)
timowarusyotoukata=timokata/syotoukata
read(8,*) atusa,ccxtawami
ccxk=(pkata*ell)/((ccxtawami-syotoukata)*ea)
sigumamoku=((e*pkata*ell/2.0d0)*(h2/2.0d0+yz))/ei
sotaigosa=((timokata-ccxtawami)/timokata)*100
!print*,"ccxk"
!print*,ccxk
!print*,"sigumamoku"
!print*,sigumamoku
!print*,"yz"
!print*,yz
!print*,"hakom"
!print*,hakom
!print*,"hakon"
!print*,hakon
!print*,"譁ｭ髱｢遨�"
!print*,dmenseki
!print*,"譖ｲ縺貞央諤ｧ"
!print*,ei
!print*,"謠帷ｮ励○繧捺妙陬懈ｭ｣菫よ焚bank"
!print*,hakok
!print*,"縺帙ｓ譁ｭ陬懈ｭ｣菫よ焚bank"
!print*,kea
!print*,"蛻晉ｭ画｢∫炊隲門､(迚�戟縺｡)"
!print*,syotoukata
!print*,"繝�ぅ繝｢繧ｷ繧ｧ繝ｳ繧ｳ逅�ｫ門､(迚�戟縺｡)"
!print*,timokata
!print*,"繝�ぅ繝｢繧ｷ繧ｧ繝ｳ繧ｳ/蛻晉ｭ画｢�(迚�戟縺｡)"
!print*,timowarusyotoukata
!print*,"逶ｸ蟇ｾ隱､蟾ｮ繝�ぅ繝｢繧ｷ繧ｧ繝ｳ繧ｳ(%)"
!print*,sotaigosa

write(7,'(1pe13.5,a,1pe13.5,a,1pe13.5,a,1pe13.5)')atusa&
&,' ',sotaigosa,' ',syotoukata,' ',timokata
end do
close(8)
end