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... 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 --~--~---------~--~----~------------~-------~--~----~ 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 -~----------~----~----~----~------~----~------~--~---