Groups > comp.soft-sys.math.mathematica > #16941

Re: Cirlce in 3D?

Path	csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!news.glorb.com!news-out.readnews.com!news-xxxfer.readnews.com!newspump.sol.net!post2.nntp.sol.net!posts.news.twtelecom.net!nnrp3.twtelecom.net!not-for-mail
From	"djmpark" <djmpark@comcast.net>
Newsgroups	comp.soft-sys.math.mathematica
Subject	Re: Cirlce in 3D?
Date	Tue, 13 May 2014 02:26:11 +0000 (UTC)
Sender	steve@smc.vnet.net
Approved	Steven M. Christensen <steve@smc.vnet.net>, Moderator
Message-ID	<lkrvs3$1t$1@smc.vnet.net> (permalink)
References	<10319860.232311.1399870151373.JavaMail.root@m03>
Lines	100
Organization	Time-Warner Telecom
NNTP-Posting-Date	13 May 2014 02:32:33 GMT
NNTP-Posting-Host	18471be1.news.twtelecom.net
X-Trace	DXC=O0aS1Y9oQMCk^S9D<aY5lLC_A=>8kQj6M;[h;PUXBgbD?SWP;HjUf@FEFiONJ7[GoFHPc4aAJMdYD
X-Complaints-To	abuse@twtelecom.net
Xref	csiph.com comp.soft-sys.math.mathematica:16941

Show key headers only | View raw

Using John Browne's Grassmann algebra application and with the three points
at Cartesian coordinates {ax, ay, az}, {bx, by, bz}, and {cx, cy, cz} I
calculated the center of the circle at:

{x -> -(4 (az (by - cy) + bz cy - by cz + 
         ay (-bz + 
            cz)) (-(ax^2 + ay^2 + az^2 - bx^2 - by^2 - bz^2) (-by cx +
             ay (-bx + cx) + ax (by - cy) + bx cy) - 
         2 (az - bz) (az by cx - ay bz cx - az bx cy + ax bz cy + 
            ay bx cz - ax by cz)) - (2 (ax^2 + ay^2 + az^2 - bx^2 - 
            by^2 - bz^2) (bz - cz) - 
         2 (az - bz) (bx^2 + by^2 + bz^2 - cx^2 - cy^2 - 
            cz^2)) (-2 (ay - by) (ay (bx - cx) + by cx - bx cy + 
            ax (-by + cy)) - 
         2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
            ax (-bz + cz))))/(-4 (az (bx - cx) + bz cx - bx cz + 
         ax (-bz + cz)) (-2 (ay - by) (ay (bx - cx) + by cx - bx cy + 
            ax (-by + cy)) - 
         2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
            ax (-bz + cz))) + 
      4 (az (by - cy) + bz cy - by cz + 
         ay (-bz + cz)) (2 (ax - bx) (-by cx + ay (-bx + cx) + 
            ax (by - cy) + bx cy) + 
         2 (az - bz) (az (by - cy) + bz cy - by cz + ay (-bz + cz)))),
  y -> (-2 (ax^2 + ay^2 + az^2 - bx^2 - by^2 - bz^2) (bz - cz) + 
     2 (az - bz) (bx^2 + by^2 + bz^2 - cx^2 - cy^2 - 
        cz^2) + (4 (az (bx - cx) + bz cx - bx cz + 
          ax (-bz + cz)) (4 (az (by - cy) + bz cy - by cz + 
             ay (-bz + 
                cz)) (-(ax^2 + ay^2 + az^2 - bx^2 - by^2 - 
                 bz^2) (-by cx + ay (-bx + cx) + ax (by - cy) + 
                bx cy) - 
             2 (az - bz) (az by cx - ay bz cx - az bx cy + ax bz cy + 
                ay bx cz - ax by cz)) - (2 (ax^2 + ay^2 + az^2 - 
                bx^2 - by^2 - bz^2) (bz - cz) - 
             2 (az - bz) (bx^2 + by^2 + bz^2 - cx^2 - cy^2 - 
                cz^2)) (-2 (ay - by) (ay (bx - cx) + by cx - bx cy + 
                ax (-by + cy)) - 
             2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
                ax (-bz + cz)))))/(-4 (az (bx - cx) + bz cx - bx cz + 
           ax (-bz + cz)) (-2 (ay - by) (ay (bx - cx) + by cx - 
              bx cy + ax (-by + cy)) - 
           2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
              ax (-bz + cz))) + 
        4 (az (by - cy) + bz cy - by cz + 
           ay (-bz + cz)) (2 (ax - bx) (-by cx + ay (-bx + cx) + 
              ax (by - cy) + bx cy) + 
           2 (az - bz) (az (by - cy) + bz cy - by cz + 
              ay (-bz + cz)))))/(4 (az (by - cy) + bz cy - by cz + 
       ay (-bz + cz))), 
 z -> (-az^2 by + by cx^2 + az^2 cy - bx^2 cy - by^2 cy - bz^2 cy + 
      by cy^2 + ax^2 (-by + cy) + ay^2 (-by + cy) + by cz^2 + 
      ay (bx^2 + by^2 + bz^2 - cx^2 - cy^2 - cz^2))/(2 (-bz cy + 
        az (-by + cy) + ay (bz - cz) + by cz)) + ((-by cx + 
        ay (-bx + cx) + ax (by - cy) + 
        bx cy) (4 (az (by - cy) + bz cy - by cz + 
           ay (-bz + 
              cz)) (-(ax^2 + ay^2 + az^2 - bx^2 - by^2 - 
               bz^2) (-by cx + ay (-bx + cx) + ax (by - cy) + 
              bx cy) - 
           2 (az - bz) (az by cx - ay bz cx - az bx cy + ax bz cy + 
              ay bx cz - ax by cz)) - (2 (ax^2 + ay^2 + az^2 - bx^2 - 
              by^2 - bz^2) (bz - cz) - 
           2 (az - bz) (bx^2 + by^2 + bz^2 - cx^2 - cy^2 - 
              cz^2)) (-2 (ay - by) (ay (bx - cx) + by cx - bx cy + 
              ax (-by + cy)) - 
           2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
              ax (-bz + cz)))))/((az (by - cy) + bz cy - by cz + 
        ay (-bz + cz)) (-4 (az (bx - cx) + bz cx - bx cz + 
           ax (-bz + cz)) (-2 (ay - by) (ay (bx - cx) + by cx - 
              bx cy + ax (-by + cy)) - 
           
           2 (az - bz) (az (bx - cx) + bz cx - bx cz + 
              ax (-bz + cz))) + 
        4 (az (by - cy) + bz cy - by cz + 
           ay (-bz + cz)) (2 (ax - bx) (-by cx + ay (-bx + cx) + 
              ax (by - cy) + bx cy) + 
           2 (az - bz) (az (by - cy) + bz cy - by cz + 
              ay (-bz + cz)))))}

I have a notebook on it but it uses the GrassmannCalculus and Presentations
Applications. Presentations has Circle3D and Disk3D primitives. The notebook
also calculates and displays random cases.


David Park
djmpark@comcast.net 
http://home.comcast.net/~djmpark/index.html 



From: Ste[hen Gray [mailto:stevebg@roadrunner.com] 

I'm looking for a neat formula to find the center of a circle in 3D through
3 points. I also need a good way to display it, preferably thickened so I
can show several and see whether they are linked, etc. To my surprise I did
not find anything on the Wolfram sites about these problems. (I have
Mathematica 7, if that matters.)

Back to comp.soft-sys.math.mathematica | Previous | Next | Find similar

Thread

Re: Cirlce in 3D? "djmpark" <djmpark@comcast.net> - 2014-05-13 02:26 +0000

csiph-web