if __name__ == "__main__":
start = time.time()
neutron_cycle = 5000;
cycle = 10;
r = 0.41;
b = 0.627;
pin_num = 3;
tot_num = pin_num ** 2;
width = b * pin_num;
fuel = Geometry('circle');
fuel.circle(r, 2*b, pin_num);
keff = np.zeros((cycle,1),float);
en_gp = 3;
## h1 = Isotope('h1', 1, [0,0.2,20], [0,0.00004,4],6.7358*1e22); # calculate macro cross-section
## o16 = Isotope('o16', 16, [0,0.0001,4], [0,0.00000003,3],3.3679*1e22);
## u235 = Isotope('u235', 235, [530,99,10], [1,0.09,4],9.3472*1e20);
## u238 = Isotope('u238', 238, [0.00002,2,9], [0.3,0.07,5],2.1523*1e22);
## pu239 = Isotope('pu239', 239, [748,269,8], [2,0.05,5],1*1e18);
h1 = Isotope('h1', 1, [0,0.294,28.966],[0,0.149,240.23], [0,0.00004,4],6.7358*1e22); # calculate macro cross-section
o16 = Isotope('o16', 16, [0,0.0001,4], [0,0.0001,4],[0,0.00000003,3],3.3679*1e22);
u235 = Isotope('u235', 235, [504.8,86.7,15.98], [271.5,131.97,152.82],[1.219,0.095,6.3],9.3472*1e20);
u238 = Isotope('u238', 238, [0.00001,2.4,9.4], [0.002,277.7,319.06],[0.3,0.07,4.825],2.1523*1e22);
pu239 = Isotope('pu239', 239, [699.34,274.32,7.90], [289.36,184.06,155.87], [1.369,1.800,0.065],0);
