[sage-support] Re: What changed?

2010-01-02 Thread Sterling
Ah, yes, I forgot the variable y. Regardless, I still get a "math domain error" that wasn't present before. -- 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 g

[sage-support] What changed?

2010-01-02 Thread Sterling
x,0,4),(y,0,4),fill=False,contours=30) c = a + b c Thanks. This isn't imperative by any means. I was just going through some old notebooks from last semester. -Sterling -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sag

[sage-support] Re: Mathematica code conversion

2009-11-23 Thread Sterling
Is this the optimal code for what I'm trying to do? On my MacBook, it takes a good minute or so before the graph appears. Not that I'm complaining... On Nov 21, 9:56 pm, Robert Bradshaw wrote: > On Nov 21, 2009, at 7:50 PM, Jason Grout wrote: > > > > > > >

[sage-support] Mathematica code conversion

2009-11-21 Thread Sterling
f(z) = 4*log(z^3)-2*log(z^3-8) g = lambda x,y: imag(f(x+y*I)) a = contour_plot(g,(x,0,4),(y,0,4),fill=False,contours=30) Robert Bradshaw suggested I use: g = lambda x,y: imag(f(x+y*I)) if y < sqrt(3)*x else float('nan') It works, but as Robert said, it isn't really pretty. An

[sage-support] Complex Integration

2009-10-19 Thread Sterling
Does SAGE support complex integration? This doesn't seem to work: z = var('z') integrate(1/z,z,-i,i) It returns an error saying the lower limit needs to be real. No rush, Sterling --~--~-~--~~~---~--~~ To post to this group, send email

[sage-support] Numpy norm

2009-09-30 Thread Sterling
When I use numpy to calculate the norm, this works fine: linalg.norm(y, ord=2) But if I want to use: linalg.norm(y, ord=inf) NameError: name 'inf' is not defined. What am I doing wrong? --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegr

[sage-support] Re: Evaluating a Jacobian

2009-09-25 Thread Sterling
Now how do I evaluate f itself at those same points. I can't seem to figure it out. On Sep 24, 9:47 pm, Jason Grout wrote: > Sterling wrote: > > How do I evaluate a Jacobian at certain values? For example, I type: > > > x1,x2,x3 = var('x1 x2 x3') > > &g

[sage-support] Evaluating a Jacobian

2009-09-24 Thread Sterling
How do I evaluate a Jacobian at certain values? For example, I type: x1,x2,x3 = var('x1 x2 x3') f1(x1,x2,x3) = 3*x1 - cos(x2*x3) - (1/2) f2(x1,x2,x3) = x1^2 - 81*(x2 + 0.1)^2 + sin(x3) + 1.06 f3(x1,x2,x3) = e^(-x1*x2) + 20*x3 + (10*pi - 3)/3 f = (f1,f2,f3) j = jacobian(f, [x1,x2,x3]) I though