So basically I have to: Firstly calculate a certain value 'e'
be = b/(n+1)
h=np.arange(1,800,1)
e1 = (h*(h+t1))/(2*(be[1]+h))
e2 = (h*(h+t1))/(2*(be[2]+h))
'e' continues up to 25. As you can imagine it is rather slow to define a new variable e#n every time. I need all this values 'e' to put in in a different formula:
Iyy1= ((1/12)*be[1]*t1**3 + be[1]*t1*e1**3) + ((1/12)*t1*h**3 + 0.25*t1*h**2 - t1*e1*h**2 + 0.25*(h**2)*(t1**2) + h*t1*e1**2 - 0.5*e1*h*t1**2) + ((2/15)*h*t1**3 + ((2/5)*h*t1)*((t1+0.5*h)+(0.5*h-(e1-0.5*t1)))**2)
Here I need to do the same thing with a range of 'e' instead of 'be'. So in the end i would have to get 25 arrays Iyy. I've tried some thing with loops, but i couldn't solve this (I'm quite a beginnner). Can Anyone help out?