Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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.

share|cite|improve this question

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.