Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: A simple Solve gives me errors
Replies: 2   Last Post: Apr 24, 2014 4:03 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Steve Gray

Posts: 103
Registered: 7/6/07
A simple Solve gives me errors
Posted: Apr 23, 2014 4:53 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

I want to find the center C of a circle in 3D where the plane is defined
by points p1,p2,p3. I wrote this:

p1 = {p1x, p1y, p1z}; (* tell it that p1,p2,p3 are vectors *)
p2 = {p2x, p2y, p2z}; (* center of circle ctr = f p1+g p2+(1-f-g)p3 *)
p3 = {p3x, p3y, p3z}; (* keep the center in the plane of p1,p2,p3 *)

ctr = f p1 + g p2 + (1 - f - g) p3;
eq1 = (ctr - p1).(ctr - p1) - (ctr - p2).(ctr - p2) := 0;
eq2 = (ctr - p1).(ctr - p1) - (ctr - p3).(ctr - p3) := 0;

I get errors of a type I've seen before:

SetDelayed::write: Tag Norm in Norm[{-p1x+f p1x+g p2x+(1-f-g) p3x,-p1y+f
p1y+g p2y+(1-f-g) p3y,-p1z+f p1z+g p2z+(1-f-g) p3z}] is Protected. >>

And a similar one. I don't know what to do about this.
I want to do
Solve[{eq1, eq2}, {f, g}]




Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.