Hi,(adsbygoogle = window.adsbygoogle || []).push({});

We assume

Max hBy(y) = α1 for y ε (a1, a2) and Min hBy(y) = α2 for y ε (a1, a2)

Where -∞ < α2 ≤ α1 < ∞, and the complex phase speed must lie in the region defined by

(CR + α1 )^2 + CI^2 ≤ µγ^2/k^2, if CR <- α1

CI^2 ≤ µγ^2/k^2, if -α1≤ CR ≤-α2,

(CR + α2 )^2 + CI^2 ≤ µγ^2/k^2, if CR >- α2

Where CI ≥ 0 and γ^2= max [hBy(y) h0y(y)] >0, h0y(y)= -8(y-a)/L^2

The region represents a rectangle of length α1 - α2 with a quarter circle on each end, with the height of the rectangle and the radius of the circles given by µ^(1/2)γ/k.

Could you please help me to write the Fortran code for this equations?

Thanks

Logi

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Help me to write fortran code

**Physics Forums | Science Articles, Homework Help, Discussion**