[adsense:responsibe:9545213979]
var('x y z') def rotacional(F): assert(len(F) == 3) return vector([diff(F[2],y)-diff(F[1],z), diff(F[0],z)-diff(F[2],x), diff(F[1],x)-diff(F[0],y)]) F1=vector([(2*x^2* y+y*z),( x* y^2-x* z^3),-( 6* x* y* z-2* x^2* y^2)]) def divergencia(F): assert(len(F) == 3) return diff(F[0],x) + diff(F[1],y) + diff(F[2],z) F2=rotacional(F1) divergencia(F2)
Añadir nuevo comentario