MathGroup Archive 1999

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Solve


steffenz at mail.fo.fh-koeln.de wrote:
> 
> Hi !
> 
> I try to use solve but mathematica neads houers eating up all memory
> but did not finish.
> Is there any mistake in my programming or is this problem to
> complicate?
> 
> bye
> 
> Steffen
> 
> Einheits[vec_]:=vec/Sqrt[vec.vec];
> Wink[vec1_,vec2_]:=ArcCos[Einheits[vec1].Einheits[vec2]];
> 
> v[alp_,bet_]:=Einheits[{Tan[bet],Tan[alp],1}];
> vlam[alp1_,bet1_,alp2_,bet2_]:=Einheits[-v[alp1,bet1]+v[alp2,bet2]];
> quell[vec_,fac_:1.3]:=Einheits[vec*{1,1,1/fac}];
> 
> Solve[{Evaluate[Wink[v[a12,b12],v[a22,b22]]==Wink[v[a11,b11],v[a21,b21
> ]]],
> Evaluate[vlam[a12,b12,a22,b22]==quell[vlam[a11,b11,a21,b21]]]},{a11,b1
> 1,a21,
>     b21}]


Hi Steffen,

Have you looked at the expressions you transfered to Solve? Probably
not: 
There are about 180 lines with three tangent-functions per line!

Try to simplify the equations first:
eq1 = Evaluate[Wink[v[a12,b12],v[a22,b22]]==Wink[v[a11,b11],v[a21,b21]]]
eq2 = Evaluate[vlam[a12,b12,a22,b22]==quell[vlam[a11,b11,a21,b21]]]

smpEq1 = Simplify[eq1]
smpEq2 = Simplify[eq2]

Solve[{smpEq1, smpEq2}, {a11,b11,a21,b21}]

Good luck!

Harald

-- 
Harald Giese
Email: giese at dkrz.de
> NEUE TELEFON-NUMMER | NEW PHONE-NUMBER <
Phone: +49 (0)40 42838 5796; Fax: +49 (0)40 5605724
Institut fuer Meereskunde der Universitaet Hamburg
(Institute of Oceanography of the University of Hamburg)
Troplowitzstrasse 7, D-22529 Hamburg


  • Prev by Date: Re: ARghh!!
  • Next by Date: Re:Exporting PRINT results from Mathematica to other applications
  • Previous by thread: Solve
  • Next by thread: ReadList