Hello,
I'm trying to recreate a plot that I did in Maple
P := (x, y) -> WeierstrassP(x + y*I, 1, 0);
PP := (x, y) -> WeierstrassPPrime(x + y*I, 1, 0);
IR := (theta, z) -> cos(theta)*Re(z) + sin(theta)*Im(z);
Gr1 := theta -> [[Re(P(x, y)), Im(P(x, y)), 0.3*IR(theta, PP(x, y))], x =
0.001 .. 3.74, y = 0.001 .. 3.74, color = [sin(2*Pi*x/3.74), 0.5,
sin(2*Pi*y/3.74)], view = [-1 .. 1, -1 .. 1, -1 .. 1], grid = [35, 35]];
plots[animate](plot3d, Gr1(t), t = 0 .. Pi - 0.1)
Here's what I have so far.
var('z')
E = EllipticCurve(QQ, [0,0,0,-1/4,0])
wp = E.weierstrass_p().laurent_polynomial()
wpp = derivative(wp,z)
f = (lambda u,v: wp(u+i*v).real(), lambda u,v: wp(u+i*v).imag(), lambda
u,v: wpp(u+i*v).imag())
parametric_plot3d(f, (0.001,3.74), (0.001,3.74))
the only problem I see is that I have no control over the view box. With
Maple I can have the parameter view = [-1 .. 1, -1 .. 1, -1 .. 1] such that
the software only shows stuff in that range.
contour_plot has the optional region parameter, so it's possible someone
already implemented something similar in another class.
Best,
RJ Acuña
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sage-support/29ddc42a-3753-436a-9c8d-a969dcac49d8n%40googlegroups.com.
