function [vx,vy,vz]=circ3(x,y,z)
	vx = 4*(x^2+y^2-1)*x*e1;
	vy = 4*(x^2+y^2-1)*y*e2;
	vz = 2*z*e3;
