I have this problem to solve. I want to compute the inclination of a plane $\theta(t)$ at every frame of a simulation given the following rule for its angular speed of rotation $\omega(t)$
\omega(t) = - \frac{\mathbf{V}\cdot\hat{\mathbf{n}}}{(||\mathbf{P}-\mathbf{P}'||)}\frac{1}{\tan(\theta(t))}
given that the starting angle \theta(t=0)=45°
The quantity \frac{\mathbf{V}\cdot\hat{\mathbf{n}}}{(||\mathbf{P}-\mathbf{P}'||)} is the ratio between the velocity of the observer along his line of sight ($\mathbf{V}\cdot\hat{\mathbf{n}}$) and the distance of the observer $\mathbf{P}$ from a point \mathbf{P}' and is computed externally at every frame.
I'm unable to figure out how to solve this equation in Mathematica, is it possible analytically?
If this is not possible, how can I implement a iterative method that can deal with the imprecision of computation of the $\frac{1}{\tan(\theta(t))} $ term?
In a first implementation I choose:
double omegaz= -speedAlongLineOfSight/((cyclopeanEye-projPoint).norm() *tan( toRadians(theta) )); theta = theta - (omegaz)*(deltaT);
but I think this continuosly sums errors, creating instability of the solution.