[sage-support] Problem with Solve

2011-05-27 Thread rpoirier
I want to solve a relativity problem using SOLVE. I used to be able to, but now I run into a problem. lp=60e9*1000 dtp=1.25*24*3600 beta=var('beta') gamma=1/(sqrt(1-beta^2)) solve(dtp*gamma==lp/(beta*c),beta) the result I expect is a numerical value, but I get the following: [beta == 733841/39600

[sage-support] Problem with solve and complex numbers: incorrect answer

2011-03-29 Thread Kunonin
Hello, i was working with real functions in real coefficients. I was doing some test to be sure it was all always real (without imaginary part), and i ask Sage to solve this: var('x',domain=RR); solve([imag(x)==0,imag(sqrt(-x^2-1))==0],x); and it answered me: [[x == -I], [x == I]] I would under

Re: [sage-support] Problem with solve() - n() -- not on maxima side

2011-02-21 Thread Francois Maltey
achrzesz wrote : SAGE: sage: s=solve(3*x^3-9*x+10==0,x,solution_dict=True) sage: [n(t[x]) for t in s] [1.06780542232902 - 1.84949324407141*I, 0.0277635108030695 + 1.24902476648341*I, -1.09556893313209 + 0.600468477588001*I] # WRONG! sage: s=solve(3*x^3-9*x+10==0,x, solution_dict=True,to_poly_

[sage-support] Problem with solve() - n() -- not on maxima side

2011-02-21 Thread achrzesz
SAGE: sage: s=solve(3*x^3-9*x+10==0,x,solution_dict=True) sage: [n(t[x]) for t in s] [1.06780542232902 - 1.84949324407141*I, 0.0277635108030695 + 1.24902476648341*I, -1.09556893313209 + 0.600468477588001*I] # WRONG! sage: s=solve(3*x^3-9*x+10==0,x, solution_dict=True,to_poly_solve='force') sage

[sage-support] problem with solve()?

2011-02-17 Thread Dmitry Shkirmanov
Hello, list, i am solving the simple system of the equations: # var("pa,pb,costheta,Ea,Eb,ma,mb,mc") solve([pa*pb*costheta==(-ma^2-mb^2+mc^2)/2 + Ea*Eb, pa==sqrt(Ea^2- ma^2), pb==sqrt(Eb^2-mb^2)],pa,pb,Eb) # Sage replies: # TypeError: unable to make sense of Maxima expression # what is wrong? --

[sage-support] Problem with "solve"

2009-08-13 Thread Mani chandra
Hi, I have a set of variables which I defined in an array and I'm trying to use "solve" to solve for these variables, but I get the following error: x = [] N = 3 for i in range(N): string = 'x' + str(i) temp_var = var(string) x.append(temp_var) Eqn = [] Eqn.append(x[0] + x[1] +