diff options
Diffstat (limited to 'src/gr/ellipsoid.maple')
| -rw-r--r-- | src/gr/ellipsoid.maple | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/src/gr/ellipsoid.maple b/src/gr/ellipsoid.maple index c9b0ae4..0c485fa 100644 --- a/src/gr/ellipsoid.maple +++ b/src/gr/ellipsoid.maple @@ -6,10 +6,11 @@ # angular coordinate system has center (U,V,W) # # direction cosines wrt angular coordinate center are (alpha,beta,gamma) -# i.e. a point has coordinates (U+alpha*r, V+beta*r, W+gamma*r) +# but Maple predefines gamma = Euler's constant, so we use (xcos,ycos,zcos) +# instead, i.e. a point has coordinates (U+xcos*r, V+ycos*r, W+zcos*r) # # then the equation of the ellipsoid is -# (U+alpha*r - A)^2 (V+beta*r - B)^2 (W+gamma*r - C)^2 +# (U+xcos*r - A)^2 (V+ycos*r - B)^2 (W+zcos*r - C)^2 # ----------------- + ---------------- + ----------------- = 1 # a^2 b^2 c^2 # @@ -18,7 +19,7 @@ # BV = B - V # CW = C - W # -eqn := (alpha*r - AU)^2/a^2 + (beta*r - BV)^2/b^2 + (gamma*r - CW)^2/c^2 = 1; +eqn := (xcos*r - AU)^2/a^2 + (ycos*r - BV)^2/b^2 + (zcos*r - CW)^2/c^2 = 1; read "../maple/util.mm"; read "../maple/codegen2.mm"; |