Problems solving using Solve
- To: mathgroup at smc.vnet.net
- Subject: [mg68809] Problems solving using Solve
- From: "akil" <akomur at wanadoo.nl>
- Date: Sat, 19 Aug 2006 00:41:17 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
I have a problem suing the following code:
Cone[beta_, v_,rcAB_,rcCD_,ax_,ay_,cx_,cy_,centerOfMassx_,centerOfMassy_] :=
Module[
{betaUsed, vtemp, fx, fy, ex, ey,
waarde,criticalRC,criticalHeight,xVerplaatsing,ePx2Py2},
betaUsed = beta ;
vtemp = v;
fline[x_] := Tan[betaUsed]x - Tan[betaUsed]vtemp;
fx = x /. Solve[fline[x] == cd[x,rcCD,cx,cy], x][[1]];
fy = cd[fx,rcCD,cx,cy];
criticalRC = (centerOfMassy - fy)/(centerOfMassx - fx);
criticalHeight = -criticalRC*fx + fy;
ex = (-ay + rcAB*ax + criticalHeight)/(rcAB + (-criticalRC));
ey = ab[ex,rcAB,ax,ay];
xVerplaatsing = ex - vtemp;
ePx2Py2 = Sqrt[xVerplaatsing^2 + ey^2];
waarde = xVerplaatsing/ePx2Py2;
ArcCos[-waarde]
];
afg =
Simplify[D[Cone[beta,v,rcAB,rcCD,ax,ay,cx,cy,centerOfMassx,centerOfMassy],v],
TimeConstraint -> 1 ];
Solve[afg== 0, v];
eg.
{ax,ay}= {-44.9376, 36.1678} and {cx,cy}= {17.8329, 122.477} and
{centerOfMassx,centerOfMassy}={-18.6338,72.3356} and rcAB=-0.166667 and
rcCD=-0.37931
Doing the Solve[afg== 0, v] takes ages (after an hour I quit it). Is there
any way I can make it work faster? Is any code in Cone such that it
unnecessarily makes it run faster? So are there any ways I can try to make
everything run faster? Using integers it runs, but it takes too long when
using real cones (made up of four points A B C D and a centerOfMass).
Akil