2704.sagews

#ezek a rajzolóeljárások

def feluletet_rajzol(phi):
    var('u,v')     
    ps=[parametric_plot3d(phi(*k),(u,0,1),thickness=4,color="red") 
     for k in [[0,u],[1,u],[u,0],[u,1]] ]
    p=parametric_plot3d(phi(u,v),(u,0,1),(v,0,1),mesh=true,opacity=0.6)
    
    #legördülő listából választhatunk, hogy a határ hányadik elemét szeretnénk megjeleniteni
    #illetve egyszerre is kirajzolhatjuk az osszeset
    @interact
    def _(n=selector([0,1,2,3]),osszeset=checkbox(True)):
        if osszeset:
            (p+sum(ps)).show(aspect_ratio=1)#,frame=False)
        else:
            (p+ps[n]).show(aspect_ratio=1)#,frame=False)
            
def tomortestet_rajzol(phi,hatarfelulet):
    var('u,v')
    ps=[parametric_plot3d(phi_gomb(*par),(u,0,1),(v,0,1),color="red") 
     for par in [[0,u,v],[1,u,v],[u,0,v],[u,1,v],[u,v,0],[u,v,1]] ]
    
    p=[parametric_plot3d(phi_gomb(*[[0,u,v],[1,u,v],[u,0,v],[u,1,v],[u,v,0],[u,v,1]][k]),
     (u,0,1),(v,0,1),color="blue",mesh=true,opacity=0.6) for k in hatarfelulet] 
    
    @interact
    def _(n=selector([0,1,2,3,4,5]),osszeset=checkbox(False)):
        if osszeset:
            (sum(p)+sum(ps)).show(aspect_ratio=1)
        else:
            (sum(p)+ps[n]).show(aspect_ratio=1)
            
            
            
#ezek pedig példák cellákra

#körgyűrű
def phi(u,v):
    a=3
    b=4
    return ((a+(b-a)*u)*cos(2*pi*v),(a+(b-a)*u)*sin(2*pi*v),0)

#félgömbfelület
def phi2(u,v):
    R=3
    return (R*sin(pi/2*u)*cos(2*pi*v),R*sin(pi/2*u)*sin(2*pi*v),R*cos(pi/2*u))

#tömör félgömb
def phi_felgomb(u,v,rho):
    R=1
    return (R*rho*sin(pi/2*u)*cos(2*pi*v),R*rho*sin(pi/2*u)*sin(2*pi*v),R*rho*cos(pi/2*u))

tomortestet_rajzol(phi_felgomb,[1,5])

n  osszeset  [removed] [removed] [removed] [removed]

feluletet_rajzol(phi2)

n  osszeset  [removed] [removed] [removed] [removed]