계산해봤는데, 왜 해가 안 나오냐? 함수 아래같이 써서, 함수 출력값이 I0*omega0=1.08e41 이랑 똑같이 나오는 임의의 omega 값 구하려고 테스트해보는데, 대충 omega=0.0005 넣으니까, 함수값 7.45e41 정도로 나오고는 그보다 작은 omega 값 넣으니, NaN 자꾸 터지네. 뭐지.... 함수 잘못 썼다 싶으면 지적하라고 원본 첨부한다. 대충 omega = 7.27e-5 rad/s 가 조건 부합하는 것 같다.
Function fuck(c, D, I0, m, omega)
variable c, D, I0, m, omega
return I0*omega + c^3*asin(mod(omega*D/c, pi/2))/(25*omega^2) - c^2*D*cos(asin(mod(omega*D/c, pi/2)))/(25*omega) + omega*D^2*m*sec(asin(mod(omega*D/c, pi/2)))
End
