Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

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|improve this question
add comment

Your Answer

 
discard

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

Browse other questions tagged or ask your own question.