Hi,
I think what Minh was trying to say is that these lines:
Stochastix wrote:
> sage: a=lambda x: alpha*x-mu1+mu2
> sage: f=lambda x: (a(x)*b-c+sqrt((a(x)*b-c)^2+4.0*a(x)*b*r*mu1))/(2*a
> (x)*mu1)
> sage: g=lambda x: (r+l-mu1*f(x))/mu2
> sage: prev=lambda x: f(x)/(f(x)+g(x))
> sage: k=lambda x:
Leo Maloney wrote:
> I'm trying to compute the inverse of a 5000 x 5000 sparse matrix. I'm
> getting an EOF error after it runs for about 5 hours, and then it
> states that sage is trying to access unallocated memory. Is there a
> way I can increase the memory for this computation? Every time I
Hello,
I think this is normal. Perhaps you meant the following (note the *
which expands the single argument into its components):
sage: for x in itertools.chain(*itertools.imap(Permutations,range(4))):
print x
[]
[1]
[1, 2]
[2, 1]
[1, 2, 3]
[1, 3, 2]
[2, 1, 3]
[2, 3, 1]
[3, 1, 2]
[3, 2, 1]
The following works (at least in sage-4.4.4):
sage: R. = QQ[]
sage: SF = SymmetricFunctions(QQ)
sage: SF.from_polynomial(x^2 + y^2 + z^2)
m[2]
sage: SF.from_polynomial(x^2 + y^2)
...
ValueError: x^2 + y^2 is not a symmetric polynomial
Cheers,
Jason
On 07/20/2010 09:23 PM, kcrisman wrote:
> It do
Hi,
On 08/23/2010 03:42 PM, robin hankin wrote:
> I tried this:
>
> roots = solve(x^3+10*x^2+11*x+8==0,x)
>
> The best I can do is
>
> N(roots[1].rhs())
>
> but this is just one at a time. How do I make N() operate on all of roots?
You may like
for r in roots:
print N(r.rhs())
or
[N
Is this what you want?
sage: v = vector(2^i for i in xrange(1,10))
sage: print v
(2, 4, 8, 16, 32, 64, 128, 256, 512)
-Jason
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more op
I'm not sure I completely understand your question, but maybe this
response will contain all the ingredients you need.
First, since you are working with polynomials, I suggest creating a
polynomial ring in the variables x and y, and then your polynomial f.
sage: R. = QQ['x,y']
sage: f = 3*x^2 + 2
Hello,
For polynomial equations, the following should work in general.
sage: R. = CC['x','y']
sage: f = x-y
sage: g = x^2 - y^2
sage: I = R.ideal([f]).radical()
sage: g in I
True
In general, to see if the equation g == 0 is implied by the equations
f1==0, f2==0, ..., fn==0 you can do
I = R.idea
On 09/29/2010 01:47 PM, Johannes wrote:
> Hi list,
> is it possible to create in a given environment in sage a new one, which
> know all in the parent defined variables and values, and i i leave it
> again, ijust get back to my old values?
Is this what you're looking for?
http://www.sagemath.org/
And two fewer than that...
sage: Zmod(5)
Ring of integers modulo 5
(Sorry, couldn't resist.)
On 01/04/2011 11:33 AM, Volker Braun wrote:
> You can even shave off two more key strokes and type:
>
> sage: ZZ.quo(5)
> Ring of integers modulo 5
>
> :-)
>
> --
> To post to this group, send email
I'm forwarding to the sage-combinat group, where the people who know
this code well have a better chance of seeing it.
On 05/12/2011 06:19 PM, Matthias wrote:
> Hi,
>
> i don't really know whether this is a bug, but at least it seems a bit
> odd to me (are alphacheck supposed to have different me
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hi Stan,
I think I saw one question you asked that hasn't been answered yet:
Stan Schymanski wrote:
>
> If I construct a list out of two other lists, I usually don't
> expect the original lists to change if I manipulate the resulting
> list. How c
Hi Matthew,
Matthew J wrote:
> Sage is great software that I rave about in pretty much all of my
> classes except for probability theory. I’d like to get some info on a
> few topics to clear some things up so that I can use these for classes
> and to post to an examples worksheet. Thanks in advan
If not, can you be more
specific about what 'save' does that does not work for you?
Best,
Jason Bandlow
--~--~-~--~~~---~--~~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support-unsu
wing works:
sage: a = float(1.0)
sage: QQ(RR(a))
1
I'm happy to open a ticket if that's the right thing to do here.
Thanks,
Jason Bandlow
--~--~-~--~~~---~--~~
To post to this group, send email to sage-support@googlegroups.com
To u
d one file at the start of my session.
Any advice is much appreciated!
Thanks,
Jason Bandlow
--~--~-~--~~~---~--~~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support-unsubscr...@googlegroup
Robert Bradshaw wrote:
> On Mar 21, 2009, at 9:06 AM, Jason Bandlow wrote:
>
>> Hi all,
>>
>> Is the following missing coercion known? I couldn't find anything on
>> trac, but there's a lot there related to coercion, so I may have
>> missed it.
&g
Hi all,
When I start up a clean version of sage 3.4 on my local machine and
enter the following into a notebook cell:
M=load('http://www.math.upenn.edu/~jbandlow/sage_data/dic_of_kst_to_G_cob_mats.sobj')
# This object is a dictionary
key = (1, Partition([1]),Partition([2]))
print key in M.keys()
William Stein wrote:
> On Thu, Mar 26, 2009 at 8:50 AM, Jason Bandlow wrote:
>> Hi all,
>>
>> When I start up a clean version of sage 3.4 on my local machine and
>> enter the following into a notebook cell:
>>
>> M=load('http://www.math.upenn.edu/~jban
David Joyner wrote:
> On Tue, May 12, 2009 at 6:43 PM, amps wrote:
>> I see that there is a function to compute the character table of the
>> symmetric group, but is there one where you input two partitions and
>> it outputs the value of the character indexed by the first partition
>> evaluated a
...or is it just me?
TIA,
Jason
--~--~-~--~~~---~--~~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support-unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.
William Stein wrote:
> On Mon, May 25, 2009 at 4:58 PM, Jason Bandlow wrote:
>> ...or is it just me?
>>
>
> it's indeed down. I'll restart all the sagenb servers right now.
> They'll be back in 5 minutes.
Thanks!
-Jason
--~--~-~--~~--
Alasdair wrote:
> This is more of a python question than a Sage question, but
> anyway...I'm trying to iterate over a list, producing a sequence of
> new lists, each of which is obtained from the original list by
> deleting one object. I've tried:
>
> for x in lst:
>lstc=copy(lst)
>print
Hello all,
Regarding doctesting, I'd like to work with the following setup:
1. Create a file work.sage (or work.py) somewhere in my home directory.
2. Start a notebook session, and attach work.sage.
3. Use the notebook for generating and staring at data, while using a
text editor to modify my co
William Stein wrote:
> On Tue, Apr 15, 2008 at 9:11 AM, Jason Bandlow <[EMAIL PROTECTED]> wrote:
>
>> Hello all,
>>
>> Regarding doctesting, I'd like to work with the following setup:
>> 1. Create a file work.sage (or work.py) somewhere in my home d
Hello,
pong wrote:
> Sorry for asking such a basic question but I still have trouble with
> it.
>
> I have created example.sage with the following content
>
> print "Hello World"
> print 2^3
>
> and load it into SAGE by load "example.sage". My understanding of the
> tutorial is that SAGE will con
32-bit binaries for 3.1.1
from sagemath.org and still had the same problem. I'm running Ubuntu
Hardy on an AMD laptop.
If I try any command that calls Maxima, I have to wait for a timeout,
and then get an error. Specific system information and a traceback are
below.
Thanks very much for
Thanks Michael,
> Check your home
> directory for any file with an accent or Umlaut and you likely found
> the culprit you need to rename.
This fixed the problem. Good luck with ecls.
Cheers,
Jason
mabshoff wrote:
>
> On Sep 15, 6:34 pm, Jason Bandlow <[EMAIL
Hi all,
A student of mine noticed the following and it looks like a bug to me
(at least with the documentation).
>From the notebook with 3.1.2 (sage prompts added for readability):
sage: plot?
File:
/home/jason/sage/local/lib/python2.5/site-packages/sage/plot/plot.py
Type:
Defi
29 matches
Mail list logo