internship formulas

`(delomega_x)/(delt) = [(2*a^2)/(a^2+c^2)-(2*a^2)/(a^2+b^2)]*omega_z*omega_y+P_osin(t)*(2*a^2)/(a^2+b^2)*omega_Z+(2*a^2)/(a^2+c^2)*omega_y+P_osin(t)+ L*tau_x`

`(delomega_y)/(delt) = [(2*b^2)/(a^2+b^2)-(2*b^2)/(b^2+c^2)]*omega_x*omega_z+P_ocos(t)*(2*b^2)/(a^2+b^2)*omega_Z-(2*b^2)/(b^2+c^2)*omega_x+P_ocos(t)+ L*tau_y`

`(delomega_z)/(delt) = [(2*c^2)/(b^2+c^2)-(2*c^2)/(a^2+c^2)]*omega_x*omega_y-P_ocos(t)*(2*c^2)/(a^2+c^2)*omega_y-P_osin(t)*(2*c^2)/(b^2+c^2)*omega_x+ L*tau_z`

`L = 15/(16*pi)[[((b^2+c^2)/(b^2*c^2)),0,0],[0,((a^2+c^2)/(a^2*c^2)),0],[0,0,((b^2+a^2)/(b^2*a^2))]]`