[sage-support] Re: solving equation system with multiple solutions

2014-06-25 Thread Dominique Laurain
I have same difficulty of not finding roots, even for one simple equation... :-) R.=QQ[] f1(x) = x^6 + x^5 + x^4 + x^3 + x^2 + x + 1 eqn = solve([f1 == 0],x) print(eqn) [ 0 == x^6 + x^5 + x^4 + x^3 + x^2 + x + 1 ] while removing x^6: R.=QQ[] f1(x) = x^5 + x^4 + x^3 + x^2 + x + 1 eqn = s

[sage-support] Re: solving equation system with multiple solutions

2014-06-25 Thread Dominique Laurain
I have same difficulty, even for one simple equation... :-) 2 R.=QQ[] 3 f1(x) = x^6 + x^5 + x^4 + x^3 + x^2 + x + 1 4 eqn = solve([f1 == 0],x) 5 print(eqn) 6 [ 0 == x^6 + x^5 + x^4 + x^3 + x^2 + x + 1 ] while removing x^6: R.=QQ[] 3 f1(x) = x^5 + x^4 + x^3 + x^2 + x + 1 4 eqn = sol

[sage-support] Re: Solving equation

2014-05-28 Thread Dima Pasechnik
On 2014-05-28, George Hokke wrote: > Hi, > what I want to do is to solve an equation in which the function contains a > numerical integral in its definition. > Something like this: > > sage: d=lambda y: numerical_integral(x**2+y,0,1)[0] > sage: d(0) > 0. > > works until here. > But no

[sage-support] Re: solving equation question --- rounding error ?

2011-01-18 Thread Keshav Kini
Hello, sage: A = matrix([[1, 0.106, 1.212], [3.8759765625, 0.04801171875, : 3.972], [3.0625, 0.09325, 3.249]]) sage: A.rank() 3 sage: A.det() 0.000 Though sage computes the rank to be 3, the determinant is negligible. Mathematica says the rank of this matrix is 2, and that its det

[sage-support] Re: solving equation question --- rounding error ?

2011-01-18 Thread Ben Edwards
> Thus (0,0,0) is the unique solution of your system. Uh... not quite 'Thus'. The system in fact has an infinite number of unique solutions, as the original poster pointed out. Though I don't know why sage converges on [0,0,0]. Also just because a second sage method gives the same result as the fi