from __future__ import division
import math
#B1=t_D*b/6(2+sigma_2/sigma_1)
Mx=200*10**6 # (N*mm)
y=[381,352,269.5,145.8,0,-145.8,-269.5,-352,-381,-352,-269.5,-145.8,0,145.8,269.5,352]
sigma_2=sigma_16=352
sigma_1=381
A=100
B1=A+0.8*149.6/6*(2+sigma_2/sigma_1)*2
Ixx=2.52*10**8
def sigma_z(y):
return Mx/Ixx*y
print "Stringer/Boom \t y(mm) \t \tsigma_z(N/mm^2)"
for i in range (0,16):
print "%2.0f \t\t %3.1f\t\t\t%3.1f "%(i+1,y[i],sigma_z(y[i]))