Inverting an equation.

Hello everyone,

I've got the following formula, with a/a* and gamma known.

a/a*=(1/m)*(2/(gamma+1))*(1+((gamma-1)/2)*m^2)^((gamma+1)/(2*(gamma-1)))

What I need are the corresponding Machnumbers M. I would like to write it in a script. Can anyone help me with this problem?

Many thanks!

Answer 1

[Walter Roberson]

This does not have a simple solution. It involves finding the roots of an expression with one component of the root being taken to the power of 1/2 * (gamma+1)/(gamma-1) . Due to the form of the expression, you cannot simply use logs to transform the expression.

Effectively you will have to use fzero() to find the answer.