input_case nil$ % 3-dimensional vectors % --------------------- load_package avector; vec u,v,w; u:=avec(u1,u2,u3); v:=avec(v1,v2,v3)$ w:=avec(w1,w2,w3)$ u/2+x*v-w; u.v; u^v; u.(v^w); vmod(v); coordinates r,theta,phi$ scalefactors(1,r,r*sin(theta)); depend {f,v1,v2,v3},r,theta,phi; % f,v1,v2,v3 depend on r,theta,phi grad(f); delsq(f); div(v); curl(v); getcsystem 'cylindrical$ g:=r*z*sin(phi)$ v:=grad(g); u:=avec(r0,a*t,t)$ % line integral lineint(v,u,t); deflineint(v,u,t,0,2*pi); % volume integral volintorder:=avec(2,0,1)$ volintegral(r**n,avec(0,0,0),avec(z/h*r,h,2*pi))$ on factor; ws; sub(n=0,ws); end;