implicit real*8(a-h, o-z)
p=50.d3
s=p/3.d0
ew=7.85d9 !s7
gw=0.450d9 !s7
hs=0.d0 !s7
es=206.d9
gs=es/2.d0/(1.d0+0.3d0)
b=80.d-3
ts=6.d-3
bn=(b-ts)+es/ew*ts
h=200.d-3
yw=h/2.d0-hs
si=ts*(h**3-(2.d0*yw)**3)/12.d0
wi=b*h**3/12.d0-si
tyz=-ew*s/b/(ew*wi+es*si)*(bn/2.d0*(yw**2-h**2/4.d0)-b*yw**2/2.d0)
t23=3.d0*s/2.d0/b/h
open(7,file='tyz')
open(8,file='t23')
open(9,file='ryz')
do i=1,9
 read*, e90
 read*, e45
 read*, e00
 eyz=e45-(e90+e00)/2.d0
 eyz=-eyz *1.d-6
 write(9,*) 0.2+0.25*(i-1), eyz*2.d0*gw
end do
do x=0.d0, 0.2d0, 0.1d0
 write(7,*) x, 0.d0 
 write(8,*) x, 0.d0 
end do
do x=0.2d0, 0.7d0, 0.1d0
 write(7,*) x, tyz
 write(8,*) x, t23
end do
do x=0.7d0, 1.7d0, 0.1d0
 write(7,*) x, -tyz
 write(8,*) x, -t23
end do
do x=1.7d0, 2.2d0, 0.1d0
 write(7,*) x, tyz
 write(8,*) x, t23
end do
do x=2.2d0, 2.4d0, 0.1d0
 write(7,*) x, 0
 write(8,*) x, 0
end do
!
close(7)
close(8)
close(9)
end




