nda activate sagetest
> sage
>
> On Sun, Aug 21, 2022 at 6:24 PM Robert Parini wrote:
>
>> Using conda on macOS 12.4 (with Apple silicon) I get the attached error
>> after installing sage with:
>>
>> conda create -n sagetest sage
>> conda activate sagetest
:
> Can you share the output of the following command?
>
> conda list -n sagetest
>
> On Sun, Aug 21, 2022 at 6:24 PM Robert Parini wrote:
>
>> Using conda on macOS 12.4 (with Apple silicon) I get the attached error
>> after installing sage with:
>>
>&
Using conda on macOS 12.4 (with Apple silicon) I get the attached error after
installing sage with:
conda create -n sagetest sage
conda activate sagetest
sage
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and s
r 5, 2019 at 1:43:23 AM UTC-8, Robert Samal wrote:
>
> I noticed the following strange behavior of
> graph6_string()/sparse6_string() functions of graphs:
>
> sage: K2=graphs.CompleteGraph(2)
> sage: P=K2.cartesian_product(K2)
>
> sage: print(P.sparse6_string())
> sage
I noticed the following strange behavior of
graph6_string()/sparse6_string() functions of graphs:
sage: K2=graphs.CompleteGraph(2)
sage: P=K2.cartesian_product(K2)
sage: print(P.sparse6_string())
sage: print(Graph(P.graph6_string()).sparse6_string())
:CoKN
:Cci
To explain: I understand, that
I observed the following weird behavior of the symbolic engine.
sage: x/x
1
sage: x^2/x
x
sage: (x^2+x)/x
(x^2 + x)/x
sage: assume(x>0)
sage: assume(x,'real')
sage: assumptions()
[x > 0, x is real]
sage: (x^2+x)/x
(x^2 + x)/x
To clarify: first, I consider the first two simplifications slightly
Indeed it works in Sage 8.4.
Thanks!
On Wednesday, October 9, 2019 at 8:34:41 AM UTC-7, Dima Pasechnik wrote:
>
> This got broken in Sage 8.5.
> (still works in 8.4)
>
>
>
> On Wed, Oct 9, 2019 at 6:09 AM David Joyner > wrote:
>
>>
>>
>> On Wed,
Sorry, F=GF(3), I made my original example shorter and didn't read it
properly.
So the full problematic code is
B=matrix(GF(3), 2,2,[1,0,1,0], sparse=True)
v=vector(GF(3), [1,1])
B.solve_right(v)
Thanks,
Robert
On Tuesday, October 8, 2019 at 5:17:59 PM UTC-7, Robert Samal wrote:
>
I am trying to solve a rather large linear systems of equations of GF(3).
As the matrices are sparse, I thought that adding "sparse=True" to the
constructor of the matrix could be of help. However, I ran to a strange
error message.
B=matrix(GF(3), 2,2,[1,0,1,0], sparse=True)
v=vector(F, [1,1])
Hi, Dima; thanks for taking the time to help out a noob -- again.
This is what I did:
sage: a,b,y,z=var('a,b,y,z')
sage: p=z*x^2+x^2-(x^2+y^2)*(a*x-2*b*y)+z*y^2+y^2
sage: type(p)
The book says that I should be able to do this:
sage: p.collect(x).collect(y)
And get this in return:
sage: 2*b*y*
I must have set some other strange default.--Rob
On 5/10/18, Jeroen Demeyer wrote:
> On 2018-05-10 19:38, Robert Gross wrote:
>> I have aliased cp to "cp -i" and "mv" to "mv -i".
>
> What do you mean this *exactly*?
>
> --
> You recei
On 5/9/18, I wrote:
> Hi,
>
> Mac OS 10.13.4. Upgrade from 8.1 to 8.2. Hangs at
>
> [ncurses-6.0.p0] config.status: creating include/ncurses_cfg.h
>
> This happened when I used "sage --upgrade" and it occurs again with "sage
> -i ncurses".
Same happened building from scratch. The relevant par
amiliar with, has its own plotting code, dunno about
Sage, in any event just reusing the splitting algorithm isn't any big
deal.
best,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and
Never use inexact values in integral terms.
Well, some experimentation shows that with keepfloat = false (the
default), it works OK. But keepfloat = true (as in Sage) it bumps into
this error.
keepfloat causes trouble in other places too ... every now and then I
think we (Maxima project) should du
ence of 'false' is a bug.
If you can make a bug report in the Maxima bug tracker, that would very
helpful. https://sourceforge.net/p/maxima/bugs
By the way I am working with Maxima 5.40+ (almost 5.41).
best,
Robert Dodier
--
You received this message because you are subscribed to the G
prevent
that conversion. That causes trouble, given the widespread implicit
assumption about exact numbers.
It's a bug of course -- perhaps you can submit a bug report to:
http://sourceforge.net/p/maxima/bugs
but you can work around it by setting keepfloat to false, or writing 1/5
in
uces abs(z) when I ask for the real part of an expression. Consider:
sage: y = var('y')
sage: assume(x, 'real')
sage: assume(y, 'real')
sage: log(x+I*y).real_part()
log(abs(x + I*y))
Best,
Robert
--
You received this message because you are subscribed to the Google Gr
eplace abs(z) with sqrt(norm(z))?
Best,
Robert
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post t
/bugs/3280 for a related bug.
HTH
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To po
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
V
ed.
Well, Sage is punting 'solve' to Maxima, and it looks like Maxima isn't
loaded correctly or something. SET-LOCALE-SUBDIR isn't anything
important, it's only setting the subdirectory for the on-line
documentation. One could delete that call, but I'm guessing one would
then
ng to classify the
type of problem, and using asksign for that. Maxima can actually get the
result without resorting to asksign, but it doesn't try that first. The
relevant function, if anybody cares, is METHOD-BY-LIMITS in
src/defint.lisp.
I dunno if this is a bug. It is needlessly clumsy bu
uming that Sage punts to Maxima for real() here. But
integral_numerical is probably not calling Maxima, right?
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving email
an experiment) introduces its own problems though.
That exposes some bugs in gruntz, and also some results which are
different, so it would be necessary to trawl through them and verify
that they're correct.
best
Robert Dodier
--
You received this message because you are subscribed to the Go
l and
domain:complex).
Reported as: https://sourceforge.net/p/maxima/bugs/3126/
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage
/x86_64-linux-gnu/libgpg-error.so.0
(0x7f627fec5000)
On Thursday, January 7, 2016 at 6:38:53 PM UTC+1, Volker Braun wrote:
>
> Can you tell us more about what eog links to?
>
> $ sage -sh command -v eog
> /usr/bin/eog
> $ sage -sh ldd /usr/bin/eog # adjust p
e^(-1)]
For the record, when I try this directly in Maxima, the equation is
solved in both cases. (I assume that Sage calls Maxima to solve
equations; is that right?) Heaven knows Maxima has > 1 bug in which
results depend on the names of variables, but I guess this isn't
one of them.
best,
I assume you meant
sage: v = P(5)
sage: v(oo)
A positive finite number
This is because the elements of QQ coerce to the parent of oo, which
is the "signed infinity ring." This is so we have
sage: P. = PolynomialRing(QQ)
sage: w = x + 5
sage: v = w - x
w(1.0)
6.00
sage: v(1.0)
5.0
Fixes to the source files?
On Tuesday, January 6, 2015 3:23:10 AM UTC-5, Volker Braun wrote:
>
> I built a new binary. Harald, can you replace the F21 binary with the new
> one on the mirrors?
>
> buildbot@build:~$ md5sum
> binaries/sage-6.4.1-x86_64-Linux-Fedora_21_x86_64.tar.gz
> cbdf8a8a1e4ee
sage -upgrade from 6.2 error on Fedora 19. Running in a root window
upgrade terminates with a permission error
Target: x86_64-redhat-linux
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man
--infodir=/usr/share/info --with-bugurl=http://bugzilla.redhat.com/bugzilla
--enable-b
I think the Mathematica interface is still broken. I'm looking into what it
would take to fix it.
What kind of sage object are we talking about?
On Sunday, 21 December 2014 17:01:47 UTC-5, Shane Scott wrote:
>
> Am I correct in thinking the sage-mathematica interface is still broken?
> If tha
and has been broken for a couple of years.
Do I understand correctly that the Mathematica interface is still not
working?
Best,
Robert
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiv
m working on it, for the record.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post
isn't entirely
consistent -- when p isn't decidable to Maxima, you can get a
partially-evaluated conditional, but not a partially-evaluated
loop (triggers an error), and various programming functions (e.g.
length, first, integerp) might act in an unexpected way. This, too,
hasn't
e(-1,-log(2))
(%i7) %, numer;
(%o7) 1.273097216447114
(%i8) quad_qags (foo, x, 2, 3);
(%o8) [1.273097216447114,1.413421842285782E-14,21,0]
Looks like Maxima handles the definite and indefinite integrals as
expected. Or perhaps I have misunderstood the problem?
Hope this helps,
Robert Dodier
-
edundatnt as I have already made this clear above. Is this a bug, or am I
> missing something?
This is a bug. If you have time, can you please report it to the
Maxima bug tracker: http://sourceforge.net/p/maxima/bugs
best,
Robert Dodier
--
You received this message because you are subscri
I do find this behavior quite surprising--diameter should be an alias
for either relative or absolute diameter, not depending on the
interval.
On Wed, Oct 22, 2014 at 2:44 AM, John Cremona wrote:
> I am trying to use the Real Interval Field (RIF), which in principle
> does exactly what I want. B
mat_function (log, mymatrix);
For large, numerical matrices, I'm sure linalg.logm is much faster.
But for small or nonnumerical matrices, maybe Maxima is useful.
Hope this helps,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support&
an solve a
relatively narrow range of problems). But I find that Maxima's
'to_poly_solve' can solve it. Maybe someone here can say how to
call it from Sage.
> In fact, the solution is: w=t+t^2
Are you sure? Assuming some value for t, plotting the expression
doesn't seem to s
would file a bug report.
http://sourceforge.net/p/maxima/bugs
best,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@
Am 08.10.2014 um 11:54 schrieb fbotana:
> Nevertheless, it seems that QEPCAD is not anymore installed/working. Try
> http://sagecell.sagemath.org/?q=xlgebu
That's a pity. I was using it [1], but I didn't find the time to learn
how to fix (and update) the broken [2] qepcad spkg.
R
Installed pre-compiled Sage-6.3.app on Mac OS X 10.9.5
At startup of Terminal Session, there is this
~$ /Applications/Sage-6.3.app/Contents/Resources/sage/sage; exit
sys:1: RuntimeWarning: not adding directory '' to sys.path since everybody
can write to it.
Untrusted users could put files in this
On 2014-09-20, Kristoffer Ryhl-Johansen wrote:
> f(x)=log(1-x)*log(1+x)/(1+x)
>>
>> f.integrate(x,0,1)
>>
> Produces a segfault when I run it on my ubuntu 14.04 computer
Fixed by Maxima commit f7921c5265 (bug in Risch code).
best
Robert Dodier
--
You received this
ed by one of the other packages loaded by
abs_integrate, but that doesn't seem to be the case; so I guess the
problem is triggered by abs_integrate itself.
I will try to investigate some more. If someone files a bug report,
that will help us track it. http://sourceforge.net/p/maxima/bugs
b
> Harald, thanks. I think now we have six or seven ways to do quadrature in
> Sage!
Yes, but most of them are QUADPACK, right?
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and s
functions in Maxima --
I find that quad_qags(f2, x, 1, 10^11) fails (with error=5, "integral
is probably divergent or slowly convergent") but
quad_qag(f2, x, 1, 10^11, 4) succeeds, likewise quad_qagi(f2, x, 1, inf)
succeeds. If Sage is indeed calling QUADPACK, perhaps at least the
error num
http://sagemath.blogspot.com/2009/12/mathematical-software-and-me-very.html
On Thu, Aug 14, 2014 at 1:14 AM, John Cremona wrote:
> When William Stein first started the project it was an acronym SAGE
> for (I think) System for Algebra and Geometry Experimentation. But
> soon it became a much wide
On Wednesday, July 16, 2014 10:25:03 AM UTC+2, I wrote:
>
> I see the following wrong results:
>
> sage: x<2 and x<1
> x < 2
> sage: x<2 or x<1
> x < 1
>
I found a way to compute these with Sage:
sage: qepcad(qepcad_formula.and_(x < 2, x < 1), vars='(x)')
x - 1 < 0
sage: qepcad(qepcad_formul
https://github.com/sagemath/sage/pull/21 aka
http://trac.sagemath.org/ticket/16672
On Thu, Jul 17, 2014 at 9:45 AM, Mahrud Sayrafi wrote:
> Hi,
>
> In this page:
> http://www.sagemath.org/doc/constructions/linear_algebra.html#eigenvectors-and-eigenvalues
> in the eigenvectors and eigenvalues sect
Also, why do I need two steps here?:
sage: solve([x==0, x!=1], x)
[[x == 0, -1 != 0]]
solve([x == 0, -1 != 0], x)
[x == 0]
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send a
Am 17.07.2014 11:32, schrieb Robert Pollak:
> In fact, I do not even know how to create the "and" situation.
> The following should be "x<=2 and x>=0":
>
> sage: Polyhedron(ieqs=[(2,-1), (0,0)]).Hrepresentation()
> (An inequality (-1) x + 2 >= 0,)
Oop
Am 16.07.2014 21:06, schrieb slelievre:
> Robert Pollak wrote:
> I see the following wrong results:
>
> sage: x<2 and x<1
> x < 2
> sage: x<2 or x<1
> x < 1
>
> The best way to manipulate logical combination of inequalities
Am 16.07.2014 20:41, schrieb Nils Bruin:
> On Wednesday, July 16, 2014 1:25:03 AM UTC-7, robert.pollak wrote:
> sage: x<2 and x<1
> x < 2
> sage: x<2 or x<1
> x < 1
>
> That's because "and" and "or" are program flow constructs in python, as
> they are in C (they have "shortcut eval
Hello!
I see the following wrong results:
sage: x<2 and x<1
x < 2
sage: x<2 or x<1
x < 1
Is this just a syntax problem? How would I enter this correctly?
Robert
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To
The very short answer is to type "make" and wait an hour or three.
On Jul 14, 2014 4:54 AM, "Oscar Alberto Castillo Felisola" <
o.castillo.felis...@gmail.com> wrote:
> Checking it out! Thank you John.
>
> --
> You received this message because you are subscribed to the Google Groups
> "sage-suppor
> After reading qepcad.py I tried to translate this to a qepcad call:
>
> qepcad(qepcad_formula.or_(qepcad_formula.and_(x-5 > 0, 3(x-1) <=
> x-5), qepcad_formula.and_(x-5 < 0, 3(x-1) >= x-5)), vars='(x)')
>
>
> 3(x-1) tries to call 3 with the argument x-1. That explains the error
> y
Integer(0), Integer(3)(x-Integer(1)) >= x-Integer(5))), vars='(x)')
TypeError: 'sage.rings.integer.Integer' object is not callable
Can you please help me to compose the correct call?
(Because of http://trac.sagemath.org/ticket/16642 I am currently using
http://sagecell.sagemath.org/, which has
rtainly a drawback.
I don't know if e.g. SymPy could solve it; I didn't try.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, s
s n increases without
bound. So I'm guessing the limit is zero.
If you need a proof, maybe you can show the integral is bounded,
therefore the limit is zero.
Sorry I can't be more helpful,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
&
ace 6.02e-23 with 602*10^-25.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post t
host: Windows 8.1
VirtualBox 4.3.12
guest: Ubuntu 14.04 LTS
Sage 6.2 Release Date 2014-05-06
Statements that gave rise to the error:
A5 = AlternatingGroup(5)
A5_sgs = A5.subgroups()
len(A5_sgs)
=>
...
RuntimeError: Gap produced error output
Error, sorry, the GAP Tables of Marks Library is not ins
MathML and presentation MathML. I tinkered with maximaMathML
a couple of months ago and it seemed to work OK (after fixing some
bugs). Write me off-list if you're interested.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group
What exactly do you mean by "simplify a real number?"
On Thu, May 29, 2014 at 8:32 AM, SiL588 . wrote:
> Unfortunately I don't know the rules of Phyton language, i just started
> using Sage notebook to do linear algebra computation.
> I think I did what you said, I assinged m a value that was the
On Tue, Apr 29, 2014 at 10:57 PM, Robert Bradshaw
wrote:
> On Tue, Apr 29, 2014 at 9:07 AM, Volker Braun wrote:
>> On Tuesday, April 29, 2014 3:58:14 PM UTC+1, Simon King wrote:
>>>
>>> Yes there is! The hook is the hash function.
>>
>>
>> CPytho
and you'd probably want to cache it.
>> ... in some cases only a trivial hash function (such as: hash of the
>> parent) should be used.
>
>
> or, better, just 1: set([ZZ(1), QQ(1)])
This is probably a bad idea--it'll lead to very poor and
hard-to-diagnose performance
simplify_sum ('sum(x^(3*k)/factorial(2*k),k,0,inf));
=> sqrt(%pi)*bessel_i(-1/2,x^(3/2))*x^(3/4)/sqrt(2)
Hope this helps,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop
19) log((m+1)/(2*m))
(%i20) limit (%, m, inf);
(%o20) -log(2)
best,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@
s and maxima. I
> can't recall the precise details right now, but after that I "learned" to
> defer substitution to the very end.
A different way to handle that is to change the floats to rational numbers,
via ratsimp. Then Maxima will be happy with expressions which contain
those
On Fri, Mar 14, 2014 at 2:03 PM, Georgios Tzanakis wrote:
>
> On Fri, Mar 14, 2014 at 4:49 PM, Robert Bradshaw
> wrote:
>>
>> Note that
>>
>> L[i][rows[i]] + j %w == 0:
>>
>> would probably be just (or nearly) as fast as
>>
>>
Note that
L[i][rows[i]] + j %w == 0:
would probably be just (or nearly) as fast as
(((L[i])[(rows[i])])+j %w)==0
If you're going to be dealing with arrays of ints you might want to
look into NumPy and/or memory views for even more speed.
On Thu, Mar 13, 2014 at 7:58 PM, Georgios Tzanakis w
to, I see that the patch is 7 months old. I'm
curious, is it in active development? I sure would like to see a solid
piecewise implementation in an upcoming Sage release.
Thanks again.
--Robert
On Tuesday, 28 January 2014 22:14:39 UTC-5, kcrisman wrote:
>
>
>
> On Tuesday, Janu
, 0), 1+x], [(0, 1), 1-x], [(1, 10),
0*x^0]], x)
This seems to evaluate correctly within its domain, but I can't plot it,
say, with this (or variations thereof):
plot(t(x), x, -4, 4)
What am I doing wrong?
--Robert
--
You received this message because you are subscribed to the Google G
uld not be enabled in the top-level make, in my opinion. Typically,
> "make -jN" makes parallel compiles within the same package (if the package
> supports it). The parallel build in sage compiles different packages in
> parallel, but each package still compiles as -j1.
>
>
How hard would it be to let "make -jN" actually work from the top-level make?
On Tue, Dec 31, 2013 at 4:57 PM, Joseph P. Skudlarek wrote:
> This is a request to update the README.txt file used when building from
> sources -- the README.txt buries the fact that "-jN" in "make -jN" is
> effectively
(%i12) asin(sqrt(3)*sqrt(2)/3) + asin(sqrt(3)/3) - asin(x) + u;
(%o12) -asin(x)+u+%pi/2
(%i13) asin(sqrt(3)*sqrt(2)/3) + asin(sqrt(3)/3) + asin(sqrt(1 - u^2)) +
asin (u);
(%o13) %pi
This is just what I got from some half-hearted hacking; I'm sure there
are serious limitations.
Good
sage: Integers(45)['t']
Univariate Polynomial Ring in t over Ring of integers modulo 45
I don't think we have linear algebra over general non-integral-domains, but
sage: R = GF(5)['x']
sage: M = random_matrix(R, 4, 4); b = random_vector(R, 4); x = M \ b
sage: M*x
(4*x^2 + x + 4, x^2 + 2*x + 4, 4*
Using a Python list is probably the fastest way to iterate over an
array of Python objects--it's a PyObject** under the hood and Cython
uses the C API calls to get at it. Your "check" might be the
bottleneck, especially if it's a Python call.
Also, no need to write this as a while loop; just use "
On Mon, Aug 19, 2013 at 6:30 PM, Dima Pasechnik wrote:
> On 2013-08-19, Vincent Knight wrote:
>> --001a1133aa8653f2ed04e4510b09
>> Content-Type: text/plain; charset=ISO-8859-1
>>
>> Thanks for the answer kcrisman but I'm afraid I'm still not sure I
>> understand.
>>
>> If by 'unsigned infinity' y
...)).
'integrate is a formal integral -- it doesn't invoke the code to solve
definite integrals -- so it won't bump into that error.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe
On 2013-07-18, Ed Scheinerman wrote:
> sage: sum(1/binomial(n,k),k,0,n)
> (n + 1)*2^(-n)
>
> and that answer is wrong.
That's a bug in Maxima's simplify_sum -- reported as bug # 2614.
https://sourceforge.net/p/maxima/bugs/2614/
best
Robert Dodier
--
You received this m
oat called a bigfloat,
which is not a CL type.
HTH
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups
I've recently had issues trying to log in to Sage via sagenb.org
I had set up a Sage account using my Google account. Now when I attempt to log
in it takes me to a page where I choose which Google account I want to use. So
far so good. But when I click on an account, it takes me to the new acc
nno what might be available in Sage
proper.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.c
sage: A = random_matrix(GF(2), 1, 1)
sage: A.det()
1
sage: b = random_vector(GF(2), 1)
sage: %time x = A \ b
CPU times: user 1.61 s, sys: 0.06 s, total: 1.67 s
Wall time: 1.67 s
sage: A * x == b
True
On Wed, Apr 17, 2013 at 1:45 PM, Juan Grados wrote:
> I have the equation Ax=b where
The syntax "R. = QQ[]" creates a polynomial ring in two
variables, with generators A and d (bound to the current session). A^d
is not a polynomial in A and d over QQ.
sage: R.=QQ[]
sage: R
Multivariate Polynomial Ring in A, d over Rational Field
On Sun, Mar 17, 2013 at 12:51 AM, Rolandb wrote:
>
ists makes more sense, I know a
lot of people who have switched from MATLAB to those tools (of which
Sage is a superset). It also depends on how much freedom you'll
have--if you get a job where everyone uses MATLAB you might not have
much choice but if you're doing research than you can use
ma source on Linux.
(%i1) load (to_poly_solve);
Loading maxima-grobner $Revision: 1.6 $ $Date: 2009-06-02 07:49:49 $
(%o1)
/home/robert/maxima/maxima-git/maxima-code/share/to_poly_solve/to_poly_s\
olve.mac
(%i2) to_poly_solve ([(a*x+b*y)*x*y/c=1,3*log(a + b + c) -
log(27*a*b*x*y)],[x,y]);
Maxima
Hi,
I upgraded to 5.7, and I get an error from zn_poly-0.9.p9 when running the
quick self-test: nuss_mul()... FAIL!.
This is on a Mac OX 10.6. I did manage to install the rest of sage-5.7,
and I can start sage successfully, so I can try to test anything that
anyone can suggest to pinpoint the p
It's cPickle with a capital P.
On Wed, Feb 20, 2013 at 2:30 AM, akhil wrote:
> Hello,
>
>
> I want to use cpickle to store a matrix object in a text file. Sample code
> is as follows:
>
> A = matrix(GF(2),2,3) #creating a 2 * 3 matrix having all entries
> zero
>
> import cpickle as pickl
). So that
would work if the exponent is an integer. However, I doubt if Maxima can
do anything unless the exponent is a literal integer, so the question
seems pointless. I'd have to look at it again before figuring out if the
question could be skipped.
best
Robert Dodier
--
You received thi
1000 bits, not just
53 bits (as it would be if 1.2 was parsed to 53 bits then passed in to
RealField(1000)).
For all of the above, see RR?, RDF?, etc. for (lots!) more documentation.
- Robert
On Wed, Jan 23, 2013 at 8:59 AM, LFS wrote:
> Hi - I would appreciate if someone could point
Setting xmin/xmax for parametric_plot doesn't seem to do anything, but
ymin/ymax work as expected. What am I doing wrong?
t = var('t')
parametric_plot( (cos(t), sin(t)), (t, 0, 2*pi), xmin=-2, xmax=2, ymin=-2,
ymax=2)
--
You received this message because you are subscribed to the Google Groups
entations of these approaches in
Sage or any upstream project? (e.g. PARI/GP, Singular, I don't know.)
best,
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To post to this group, send email to sage-support@googlegrou
the basis for an implementation? Pointers to any other
resources would be interesting.
best
Robert Dodier
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To post to this group, send email to sage-support@googlegroups.com.
To unsubs
%if(?%and(-%pi/2 < parg(sqrt(4*c+b^2)+b),
parg(sqrt(4*c+b^2)+b) <= %pi/2),
[d = (b*sqrt(4*c+b^2)+2*c+b^2)/2],%union()))
I didn't check the result; sorry about that.
best
Robert Dodier
--
You received this message because you are su
On Wed, Oct 17, 2012 at 9:37 PM, Eric Kangas wrote:
> code:
>
> b = 11^2
>
> a = b^2
>
> pri = [int(is_prime(i)) for i in range(a)]
>
> j = [i for i in range(a)][b+1:a:b]
>
> k = [i for i in range(a)][(b*2)+1:a:b]
>
> j.insert(0,0)
>
> k.insert(0,b)
>
> m = [matrix(QQ,sqrt(a)/(b/sqrt(b)),pri[j[i]:
What is eval? It looks like a list of eigenvalues, rather than a single
one. Could you show how it is computed?
On Wed, Oct 17, 2012 at 4:55 PM, Eric Kangas wrote:
> Of course complex.
>
> Here is what happens when I try using the embeddings() function:
>
> sage: eval[0].embeddings(CC)
>
>
When you say "plot these values", do you mean as real or complex
values? To do so you need to choose an embedding, e.g.
sage: K. = QQ[sqrt(5)]; K
Number Field in sqrt5 with defining polynomial x^2 - 5
sage: K.embeddings(CC)
[
Ring morphism:
From: Number Field in sqrt5 with defining polynomial
On Thu, Sep 27, 2012 at 1:41 PM, Tom Boothby wrote:
> On Thu, Sep 27, 2012 at 12:55 AM, Robert Bradshaw
> wrote:
>
>> I would be extremely
>> surprised if any Sage developer morally objects to you licensing this
>> output as you wish (though opinions may vary w
On Wed, Sep 26, 2012 at 10:56 PM, Geoffrey Irving wrote:
> On Wed, Sep 26, 2012 at 10:42 PM, Robert Bradshaw
> wrote:
>> On Wed, Sep 26, 2012 at 8:54 PM, Geoffrey Irving wrote:
>>> On Wed, Sep 26, 2012 at 6:03 PM, Robert Bradshaw
>>> wrote:
>>>> On Wed
1 - 100 of 1029 matches
Mail list logo