Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

answered 2011-11-14 02:53:09 -0500

DSM gravatar image

Hi! I think you must have assigned one too many of your variables (like assigning a value to RB) to come up with your last expression. Also, I don't see what predefined, lowercase q you're referring to.

In any case, here's roughly what I'd do. First, define the variables and the equation, and then solve it:

sage: var("F, Q, M, RB")
(F, Q, M, RB)
sage: deg = pi/180;
sage: 
sage: eq = -2*F*cos(45*deg)-3*Q+4*RB*cos(45*deg)-RB*sin(45*deg)-M == 0
sage: eq
-sqrt(2)*F + 3/2*sqrt(2)*RB - M - 3*Q == 0
sage: 
sage: solve(eq,RB)
[RB == 1/3*sqrt(2)*M + sqrt(2)*Q + 2/3*F]
sage: for sol in solve(eq,RB):
....:         print sol.rhs()
....: 
1/3*sqrt(2)*M + sqrt(2)*Q + 2/3*F

Then solve it. By default, solve returns a list of equations, so you can loop over them like I did or index into them using [0], [1], [2], etc. I prefer the solution_dict style:

sage: solve(eq,RB, solution_dict=True)
[{RB: 1/3*sqrt(2)*M + sqrt(2)*Q + 2/3*F}]
sage: sols = solve(eq,RB, solution_dict=True)
sage: sols[0][RB]
1/3*sqrt(2)*M + sqrt(2)*Q + 2/3*F

in which you get a list of dictionaries, with each dictionary corresponding to a possible solution, and then you simply use the variable name as the index. We can then extract a solution:

sage: rb = sols[0][RB]
sage: rbval = rb.subs(M=23, Q=6, F=1/4)
sage: rbval
41/3*sqrt(2) + 1/6

and I've substituted some basically random M, Q, and F values into it. This can then be evaluated numerically (the evalf equivalent) by using the numerical_approx method:

sage: rbval.n()
19.4942520190990
sage: rbval.numerical_approx()
19.4942520190990
sage: rbval.n(100)
19.494252019098965666956412564

There are lots of ways to turn a symbolic expression into a "number", though. See my answer to this question on a similar subject.

Does that make sense?