let hyp a b = sqrt (a**2.0 + b**2.0)
let gamma a b = (hyp a b) - a
let delta a b = (hyp a b) - b
let e = -(c + d) - 2.0 * sqrt (c * d)
let (e,f) = (coefft h1 h2)
(e * x**2.0) + (f * x) + h2
let AngArray = [| for i in 0 .. imax -> float(i) / float(imax) |]
|> Array.map (fun x -> dmax x 700.0 250.0)
|> Array.iter2 (fun x y -> printfn "%f\t%f" (x * 90.0) y ) AngArray