[x] There are repeated parts in my formula, to avoid this we have functions ;)
[x] The formula in http://www.edwilliams.org/avform.htm#flat is defined for two points p0=(lon0,lat0) (which is the fixed point) and p=(lon,lat). So then when you have the formula you should call it with p1,p2 first
and p1,p4 then compare.
[x] Get your parameters like radius from the input spheroid and note handwritten
r.e. use get_radius<>(spheroid)