sigma=40*10**6 # in Pa Tensile stress
E=69*10**9 #Modulus of elaticity in pa
Ys=0.3 #Specific surface energy in N/m**2
a=2*E*Ys/(math.pi*sigma**2)
print"Maximum lemgth of a surface flaw without fracture is ",round(a*10**6,1),"micro m"
T=800.0+273.0 # Temperature in K
import math
L_M=24*10**3
t=math.pow(10,(L_M/T)-20)
print"Time to rupture is ",round(t,0),"h"