Ask Your Question

Revision history [back]

Problem solving a system of equations

Hello everyone,

Here I copy my code in sage.

var('r2 si co r12 r22 r32 d32')

eq1 = r12==r2d32(1-si(co+sqrt(3(1-co*co)))/2)/2

eq2 = r22==r2d32(1-si(co-sqrt(3(1-co*co)))/2)/2

eq3 = r32==r2d32(1+si*co)/2

eq1.show()

eq2.show()

eq3.show()

solve([eq1,eq2,eq3],r2,si,co)

And this is the answer that I get when I execute:

Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1013, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 1, in <module> File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/symbolic/relation.py", line 1050, in solve sol_list = string_to_list_of_solutions(repr(s)) File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/symbolic/relation.py", line 580, in string_to_list_of_solutions v = symbolic_expression_from_maxima_string(s, equals_sub=True) File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/calculus/calculus.py", line 2159, in symbolic_expression_from_maxima_string raise TypeError("unable to make sense of Maxima expression '%s' in Sage"%s) TypeError: unable to make sense of Maxima expression '[if((-pi/2<parg(-((3<em>_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2)))and(parg(-((3_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==-((2_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2_SAGE_VAR_r32+2_SAGE_VAR_r22+2_SAGE_VAR_r12)/(3_SAGE_VAR_d32),_SAGE_VAR_si==-(2sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union()),if((-pi/2<parg(((3< em="">_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2)))and(parg(((3_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==((2_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2_SAGE_VAR_r32+2_SAGE_VAR_r22+2_SAGE_VAR_r12)/(3_SAGE_VAR_d32),_SAGE_VAR_si==(2sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union())]' in Sage

Problem solving a system of equations

Hello everyone,

Here I copy my code in sage.

var('r2 si co r12 r22 r32 d32')

d32') eq1 = r12==r2*d32*(1-si*(co+sqrt(3*(1-co*co)))/2)/2 eq2 = r22==r2*d32*(1-si*(co-sqrt(3*(1-co*co)))/2)/2 eq3 = r32==r2*d32*(1+si*co)/2 eq1.show() eq2.show() eq3.show() solve([eq1,eq2,eq3],r2,si,co)

eq1 = r12==r2d32(1-si(co+sqrt(3(1-co*co)))/2)/2

eq2 = r22==r2d32(1-si(co-sqrt(3(1-co*co)))/2)/2

eq3 = r32==r2d32(1+si*co)/2

eq1.show()

eq2.show()

eq3.show()

solve([eq1,eq2,eq3],r2,si,co)

And this is the answer that I get when I execute:

execute:

Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1013, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 1, in <module> File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/symbolic/relation.py", line 1050, in solve sol_list = string_to_list_of_solutions(repr(s)) File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/symbolic/relation.py", line 580, in string_to_list_of_solutions v = symbolic_expression_from_maxima_string(s, equals_sub=True) File "/ext/sage/sage-8.1/local/lib/python2.7/site-packages/sage/calculus/calculus.py", line 2159, in symbolic_expression_from_maxima_string raise TypeError("unable to make sense of Maxima expression '%s' in Sage"%s) TypeError: unable to make sense of Maxima expression '[if((-pi/2<parg(-((3<em>_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2)))and(parg(-((3_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==-((2_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2_SAGE_VAR_r32+2_SAGE_VAR_r22+2_SAGE_VAR_r12)/(3_SAGE_VAR_d32),_SAGE_VAR_si==-(2sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union()),if((-pi/2<parg(((3< em="">_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2)))and(parg(((3_SAGE_VAR_r22-3_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==((2_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2_SAGE_VAR_r32^2+((-2_SAGE_VAR_r22)-2_SAGE_VAR_r12)_SAGE_VAR_r32+2_SAGE_VAR_r22^2-2_SAGE_VAR_r12_SAGE_VAR_r22+2_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2_SAGE_VAR_r32+2_SAGE_VAR_r22+2_SAGE_VAR_r12)/(3_SAGE_VAR_d32),_SAGE_VAR_si==(2sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union())]' '[if((-pi/2<parg(-((3*_SAGE_VAR_r22-3*_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2)))and(parg(-((3*_SAGE_VAR_r22-3*_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==-((2*_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2*_SAGE_VAR_r32+2*_SAGE_VAR_r22+2*_SAGE_VAR_r12)/(3*_SAGE_VAR_d32),_SAGE_VAR_si==-(2*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union()),if((-pi/2<parg(((3*_SAGE_VAR_r22-3*_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2)))and(parg(((3*_SAGE_VAR_r22-3*_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2))<==pi/2),[_SAGE_VAR_co==((2*_SAGE_VAR_r32-_SAGE_VAR_r22-_SAGE_VAR_r12)*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(2*_SAGE_VAR_r32^2+((-2*_SAGE_VAR_r22)-2*_SAGE_VAR_r12)*_SAGE_VAR_r32+2*_SAGE_VAR_r22^2-2*_SAGE_VAR_r12*_SAGE_VAR_r22+2*_SAGE_VAR_r12^2),_SAGE_VAR_r2==(2*_SAGE_VAR_r32+2*_SAGE_VAR_r22+2*_SAGE_VAR_r12)/(3*_SAGE_VAR_d32),_SAGE_VAR_si==(2*sqrt(_SAGE_VAR_r32^2+((-_SAGE_VAR_r22)-_SAGE_VAR_r12)*_SAGE_VAR_r32+_SAGE_VAR_r22^2-_SAGE_VAR_r12*_SAGE_VAR_r22+_SAGE_VAR_r12^2))/(_SAGE_VAR_r32+_SAGE_VAR_r22+_SAGE_VAR_r12)],union())]' in Sage

Sage