Surface Exploration

xxxxxxxxxx

 

1

def c(r,t): return abs(r*sin(t))

2

s1='<left style=%s>'%'color:#3636ff; font-size:150%;'

3

s2='$px=%s$<br/>$py=%s$<br/>$pz=%s$</left>'

4

@interact

5

def _(n=[2..8],pp=[20,30,..,80],

6

      cm=['cool','hsv','bwr']):

7

    var('u,v,r,t') 

8

    x=real_part((u+I*v)^n)

9

    y=imag_part((u+I*v)^n)

10

    uv={real_part(u):u,imag_part(u):0,

11

    real_part(v):v,imag_part(v):0}

12

    xx,yy=x.subs(uv),y.subs(uv)

13

    rt={u:r*cos(t),v:r*sin(t)}

14

    px,py=xx.subs(rt),yy.subs(rt)

15

    pz=r*cos(t)

16

    pretty_print(html(s1+s2%(px,py,pz)))

17

    p=parametric_plot3d(

18

        (px,py,pz),(r,0,1),(t,0,4*pi),

19

        plot_points=pp,color=(c,colormaps[cm]))

20

    p.show(frame=False)

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages