On Tue, Apr 22, 2008 at 7:59 AM, Kai <[EMAIL PROTECTED]> wrote:
>
> Ok, I forgot to mention that i'd like to plot the real and imaginary
> parts or the absolute value of the function. With your help i got this
> far:
>
>
> sage: E = pari(['1', 'I'])
>
> sage: f = lambda a,b: E.ellwp(a+b*I)
> sage: g= lambda a,b: real(E.ellwp(a+b*I))
> sage: P = plot3d(g,(0.1,0.9),(0.1,0.9), adaptive=True,
> color=rainbow(60,'rgbtuple'))
> sage: P.show(figsize=[10,10])
>
> The resulting plot is already close to what i had in mind in the first
> place. Here I chose the interval (0.1,0.9)x(0.1,0.9) because it
> contains no poles of the weierstrass p function. For larger intervals,
> containing poles the plot becomes somehow biased, since one axis
> reaches nearly to infinity. My question now is how i can "delimit" the
> resulting plot, so that for example every value larger than 10 wont be
> shown. I hope I explained that comprehensible...
On obvious way would be to refine f and g. E.g.,
def f (a,b):
z = real(E.ellwp(a+b*I))
if abs(z) <= 10: return z
else: return 10
>
> Something else I have to remark is that the function ellwp() gives
> strange results for lattice points. If I'm not mistaken,
>
>
> sage: E = pari(['1', 'I'])
> sage: E.ellwp(1+I)
> -1/2*I
>
> for example doesnt make sense, since the weierstrass p-function has a
> pole at 1+I.
>
> Thanks,
>
>
> Kai
>
>
>
> >
>
--
William Stein
Associate Professor of Mathematics
University of Washington
http://wstein.org
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
