Campo vectoria en 3D a traves de un cilindro con SAGE.
var ('r ph z')
S=Cylindrical('radius', ['azimuth', 'height'])
S.transform(r, 0, z)
p1=plot_vector_field3d(S.transform(r, 0, z), (r,-4,4), (ph,0,2*pi), (z,0,4))
p2=cylindrical_plot3d(2,(ph,0,2*pi),(z,0,4), opacity=0.2)
show(p1+p2)
Añadir nuevo comentario