I thought it would be cool to see a vector field plot with the real
and imaginary parts, but for some reason it doesn't work. Anyone know
why not?
i=CDF.0
tau=i
Om=pari([1,tau])
var('x,y')
def g(a,b):
q = CDF(Om.ellwp(a+b*i))
return (float(q.real()),float(q.imag()))
plot_vector_field(g(x,y), (.1,.9), (.1,.9))
On May 20, 6:11 am, Marshall Hampton <[EMAIL PROTECTED]> wrote:
> If you force a cast to real type like RDF it works, for example:
>
> sage: def f(a,b): return RDF(abs(Om.ellwp(a+b*i)))
>
> instead of what you have. This may not be the optimal solution but it
> works.
>
> -M. Hampton
>
> On May 20, 1:19 am, Kai <[EMAIL PROTECTED]> wrote:
>
> > Hi there,
> > i'm having yet another problem with plotting the Weierstrass-p-
> > Function. This time i'd like do create a contour plot of it. I did the
> > following:
>
> > sage: i=CDF.0
> > sage: tau=i
> > sage: Om=pari([1,tau])
> > sage: def f(a,b): return abs(Om.ellwp(a+b*i))
> > sage: T=contour_plot(f, (0.1,0.9), (0.1,0.9), cmap='autumn',
> > contours=10)
> > sage: T.show()
>
> > and got the error message: "ValueError: setting an array element with
> > a sequence". It seems as if calling the contour_plot(..) function
> > works fine. The error message only appears after trying to show the
> > plot.
>
> > Thanks,
> > Kai
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
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
-~----------~----~----~----~------~----~------~--~---
