Program: MKCIRCUMCIRCLESPH
mkcircumcirclesph......determine circumcircle parameters of a spherical triangle
call: [ccsph,rsph]=mkcircumcirclesph(a,b,c);
a,b,c: (numeric doublets) [deg]
sphercial coordinates of points A, B, C that define a
spherical triangle. The radius of the sphere is assumed to
be 1.
a(1) is longitude of point A
a(2) is latitude of point A
b,c analoguous.
result: cc: (numeric doublet) [deg]
spherical coordinates of circumcenter point
c(1) is longitude
c(2) is latitude
radius is implied to be 1.
rsph: (number) [deg]
radius of the circumcircle, given as angle on the unit
sphere. Thsi radius is computed as distance angle between CC
and point A.
Computation uses the cartesian computation in MKCIRCUMCIRCLE: The circle
as such is the same as in cartesian coordinates, only the center point
and radius are different since they must define a small circle on a
sphere.
Martin Knapmeyer, 03.12.2004
Read M-File Source Code
eof.