Re: [sage-support] sagetex and stereographic plots

2010-03-12 Thread Dan Drake
On Fri, 12 Mar 2010 at 01:27PM +0100, G. Damm wrote: > is it possible to make stereographic 3d-plots with sagetex? > I'd want to make a beamer-presentation with these plots and want to > help my students see the 3d. Can Sage make stereographic 3-d plots? If there's a way to get a PNG image of suc

[sage-support] test if an expression is polynomial

2010-03-12 Thread ma...@mendelu.cz
Dear support the following two commands can be used to test, if a function is a polynomial or quotient of two polynomials in x variable. sage: (x^3+x+3+1/(x))._maxima_().polynomialp([x]).sage() 0 sage: ((x^3+x+3+1)/(x^4+3))._maxima_().polynomialp([x],"constantp", "integerp").sage() 1 Is there a

[sage-support] Re: convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread Yann
And I guess the answer to Paul's question is then: sage: (sinh(log(t)))._maxima_().exponentialize().sage() 1/2*t - 1/2/t sage: (cos(log(t)))._maxima_().exponentialize().sage() 1/2*e^(-I*log(t)) + 1/2*e^(I*log(t)) -- To post to this group, send email to sage-support@googlegroups.com To unsubscri

[sage-support] Re: convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread ma...@mendelu.cz
On 12 bře, 21:19, Guillaume wrote: > > No, AFAIK, nothing other than explicit substitution with .subs(). > > Hello, > > there are a few weird results. I'd like to solve this homogenous edo : > >  $tx'=x+\sqrt{x^2+y^2}$. > > using x=tu > > sage: t=var('t') > sage: x(t) = function('x',t) > sage: id

[sage-support] Re: convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread ma...@mendelu.cz
On 12 bře, 16:48, Burcin Erocal wrote: > On Fri, 12 Mar 2010 15:23:43 +0100 > > Paul Zimmermann wrote: > > is there a way in Sage to convert expressions involving trigonometric > > or hyperbolic functions to exponentials, like the convert/exp > > function of Maple? > > > > convert(sinh(log(t)),

[sage-support] Re: convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread Guillaume
> No, AFAIK, nothing other than explicit substitution with .subs(). Hello, there are a few weird results. I'd like to solve this homogenous edo : $tx'=x+\sqrt{x^2+y^2}$. using x=tu sage: t=var('t') sage: x(t) = function('x',t) sage: id(t)=t sage: u=function('u',t) sage: d=diff(u*id,t) appar

[sage-support] Re: possible error in manual

2010-03-12 Thread John H Palmieri
On Mar 12, 9:18 am, Michael Beeson wrote: > On the pagehttp://www.sagemath.org/doc/reference/sage/symbolic/expression.html > > the first example has eqn.subs(x==5)  and I think it should be > eqn.subs(x=5) It seems to work either way, and there are examples in the documentation using both: type

[sage-support] Re: norm or absolute value

2010-03-12 Thread Yann
Hi, Is this helping? sage: var('a,b,z') (a, b, z) sage: f=a*z+i*b*z^2 sage: f.norm() b*z^2*conjugate(b)*conjugate(z)^2 - I*a*z*conjugate(b)*conjugate(z)^2 + I*b*z^2*conjugate(a)*conjugate(z) + a*z*conjugate(a)*conjugate(z) sage: f.norm().full_simplify() b^2*z^4 + a^2*z^2 sage: f.norm().factor() (b

[sage-support] norm or absolute value

2010-03-12 Thread Michael Beeson
sage: var('z'); var('a');var('b'); sage: F = a*z + i*b*z^2 sage: Now F.norm() or something should give me (a*z)^2 + (b*z^2)^2, but I can't find a command to do that. I want to do this when F is a polynomial of degree 6 with complex rational coefficients to eliminate i and produce a polynomial o

[sage-support] Re: plot3d and expression evaluation

2010-03-12 Thread stefan
I think I've got a vague idea what happens here... The same thing also happens when I define the max() function myself: sage: def my_max(x,y): sage: if(x>y): return x sage: else: return y sage: fermi2(x,y,d,L) = 1 - 1/( exp( ( my_max(abs(x),abs(y))-L) /d) + 1) sage: fermi2 (x, y, d, L)

Re: [sage-support] newbie to cython: wrap int func( int n, float x[] )

2010-03-12 Thread Robert Bradshaw
On Mar 12, 2010, at 9:05 AM, gerhard wrote: Trying to wrap an existing library. I managed to at least get started with a .spyx file as follows: cdef extern from "stdlib.h": void *malloc(size_t size) int free(void*) int sizeof() cdef extern from "func.h

[sage-support] possible error in manual

2010-03-12 Thread Michael Beeson
On the page http://www.sagemath.org/doc/reference/sage/symbolic/expression.html the first example has eqn.subs(x==5) and I think it should be eqn.subs(x=5) -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr.

[sage-support] Re: plot3d and expression evaluation

2010-03-12 Thread stefan
On Mar 12, 5:47 pm, Harald Schilly wrote: > I'm not sure but there might be a bug or problem evaluating the > expression. Anyways, going the "pure" python way works: > > sage: def fermi(x,y,d,L): return  1 - 1/( exp( ( max(abs(x),abs(y))- > L) /d)  + 1) > > sage: plot3d(lambda x,y : fermi(x,y,0.

[sage-support] newbie to cython: wrap int func( int n, float x[] )

2010-03-12 Thread gerhard
Trying to wrap an existing library. I managed to at least get started with a .spyx file as follows: cdef extern from "stdlib.h": void *malloc(size_t size) int free(void*) int sizeof() cdef extern from "func.h": int func( int n, float* x ) cdef doub

[sage-support] Re: plot3d and expression evaluation

2010-03-12 Thread Harald Schilly
On Mar 12, 5:39 pm, stefan wrote: > I have encountered a somewhat strange problem ... Ah, I got it, the max function is evaluated during the definition of the function, look: sage: fermi(x,y,d,L) = 1 - 1/( exp( ( max(abs(x),abs(y))-L) /d) + 1) sage: fermi (x, y, d, L) |--> -1/(e^(-(L - abs(x))/

[sage-support] Re: plot3d and expression evaluation

2010-03-12 Thread Harald Schilly
On Mar 12, 5:39 pm, stefan wrote: > Basically, I am trying to plot this expression: > > #sage> fermi(x,y,d,L) = 1 - 1/( exp( ( max(abs(x),abs(y))-L) /d)  + 1) I'm not sure but there might be a bug or problem evaluating the expression. Anyways, going the "pure" python way works: sage: def fermi(

[sage-support] plot3d and expression evaluation

2010-03-12 Thread stefan
Hello group, I have encountered a somewhat strange problem in plotting a simple function. It seems to be related to the issues described in the tutorial section "Some Common Issues with Functions", but the lambda- function trick does not work here. I have uploaded a worksheet with what I have bee

[sage-support] sagetex and stereographic plots

2010-03-12 Thread G. Damm
Hello all, is it possible to make stereographic 3d-plots with sagetex? I'd want to make a beamer-presentation with these plots and want to help my students see the 3d. Thanks, Georg signature.asc Description: This is a digitally signed message part.

Re: [sage-support] convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread Burcin Erocal
On Fri, 12 Mar 2010 15:23:43 +0100 Paul Zimmermann wrote: > is there a way in Sage to convert expressions involving trigonometric > or hyperbolic functions to exponentials, like the convert/exp > function of Maple? > > > convert(sinh(log(t)),exp); > 1 >

[sage-support] convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread Paul Zimmermann
Hi, is there a way in Sage to convert expressions involving trigonometric or hyperbolic functions to exponentials, like the convert/exp function of Maple? > convert(sinh(log(t)),exp); 1 t/2 - ---

Re: 回复: [sage-support] Re: 3D plot in sage

2010-03-12 Thread wxuyec
Actually I do need them to be plotted by color.it seems that sage is not designed for this way.could anyone give some information aboutwhat open source software can do this?Thanks a lot!YC - 原文 - 发件人: Jason Grout 主 题: Re: 回复: [sage-support] Re: 3D plot in sage时 间: 2010年3月12日 01:26:30On

Re: 回复: [sage-support] Re: 3D plot in sage

2010-03-12 Thread Jason Grout
On 03/12/2010 03:51 AM, Marshall Hampton wrote: > One option is several implicit plots, colored by value, i.e. something > like: > > var('x,y,z') > f=cos(x)*cos(y)+cos(y)*cos(z)+cos(z)*cos(x) > > imps = [] > plot_range = srange(0,3,.75) > color_range = [i/N(len(plot_range)) for i in range(len(plo

Re: 回复: [sage-support] Re: 3D plot in sage

2010-03-12 Thread Marshall Hampton
One option is several implicit plots, colored by value, i.e. something like: var('x,y,z') f=cos(x)*cos(y)+cos(y)*cos(z)+cos(z)*cos(x) imps = [] plot_range = srange(0,3,.75) color_range = [i/N(len(plot_range)) for i in range(len(plot_range))] for i,q in enumerate(plot_range): red = color_rang