1
votes

The following equation is to be solved for M by MATLAB:

(Atemp/At)^2=1/M^2*((2/(gamma+1))*(1+(gamma-1)*M^2/2))^((gamma+1)/(gamma-1))

It is not possible to solve this equation symbolically. In Maple it is easily possible to solve such an equation implicitly; now, is there also a pre-made function in Matlab that does this for me? I could program one myself, but as my skills are limited, its performance would not fit my needs.

1
For the interested ones, this is the area-Mach number relation for supersonic flow - Ingo

1 Answers

2
votes

I would try using fzero, or if that encounters problems because of complex values/infinities, fminbnd.