Practical NO.

11

̂ , evaluate 𝒓̅(0), 𝒓̅(𝝅), 𝒓̅(𝟐𝝅)

Q. Let (1) 𝒓̅(t) =4 cost 𝒊̂ + 3 sint 𝒋̂ and (2) 𝒓̅(t) =3tant 𝒊̂ +4sect 𝒋̂ + 5t 𝒌 𝟑

(%i2) r(t):=[4*cos(t),3*sin(t)];

(%o2) r(t):=[4*cos(t),3*sin(t)]

(%i3) r(0);

(%o3) [4,0]

(%i4) r(%pi/2);

(%o4) [0,3]

(%i5) r(2*%pi/3);

(%o5) [-2,3^(3/2)/2]

(%i1) r(t):=[3*tan(t),4*sec(t),5*t];

(%o1) r(t):=[3*tan(t),4*sec(t),5*t]

(%i2) r(0);

(%o2) [0,4,0]

(%i3) r(%pi/2);

tan: %pi/2 isn't in the domain of tan.

#0: r(t=%pi/2)

-- an error. To debug this try: debugmode(true);

(%i4) r(2*%pi/3);

(%o4) [-3^(3/2),-8,(10*%pi)/3]

̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

Q. If 𝒓̅(t) =t 𝒊̂ + 4cost𝒋̂ +3sint 𝒌
𝒕→𝟐

(%i5) r(t):=[t,4*cos(t),3*sin(t)];

(%o5) r(t):=[t,4*cos(t),3*sin(t)]

(%i6) limit(r(t),t,2);

(%o6) [2,4*cos(2),3*sin(2)]
(%i7) float(%);

(%o7) [2.0,-1.66458734618857,2.727892280477045]s

̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

Q. If 𝒓̅(t) =cost 𝒊̂ – cost 𝒋̂ + sint 𝒌
𝒕→𝟐

(%i8) r(t):=[cos(t),-cos(t),sin(t)];

(%o8) r(t):=[cos(t),-cos(t),sin(t)]

(%i9) limit(r(t),t,2);

(%o9) [cos(2),-cos(2),sin(2)]

(%i10) float(%);

(%o10) [-0.4161468365471424,0.4161468365471424,0.9092974268256817]

̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

Q. If 𝒓̅(t) =sin(t/2) 𝒊̂ + cost(2t/3) 𝒋̂ + tan(5t/4) 𝒌
𝒕→𝝅

(%i11) r(t):=[sin(t/2),cos(2*t/3),tan(5*t/4)];

(%o11) r(t):=[sin(t/2),cos((2*t)/3),tan((5*t)/4)]

(%i12) limit(r(t),t,%pi);

(%o12) [1,-1/2,1]

̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

Q. If 𝒓̅(t) =t3 𝒊̂ + sin(𝝅𝒕/𝟐) 𝒋̂ +log(t+2) 𝒌
𝒕→−𝟏

(%i13) r(t):=[t^3,sin(%pi/2),log(t+2)];

(%o13) r(t):=[t^3,sin(%pi/2),log(t+2)]

(%i14) limit(r(t),t,-1);

(%o14) [-1,1,0]

̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

Q. If 𝒓̅(t) =sin(t/2) 𝒊̂ + cos(2t/3) 𝒋̂ +tan(5t/4) 𝒌
𝒕→𝝅

(%i15) r(t):=[sin(t/2),cos(2*t/3),tan(5*t/4)];

(%o15) r(t):=[sin(t/2),cos((2*t)/3),tan((5*t)/4)]

(%i16) limit(r(t),t,%pi);

(%o16) [1,-1/2,1]

Q. If 𝒓̅(t) =(sin(t-2)2/t-2)) 𝒊̂ + (t2-4)/(t-2) 𝒋̂ , then find 𝒍𝒊𝒎 𝒓̅(𝒕)

𝒕→𝟐

(%i17) r(t):=[sin(t-2)^2/(t-2),(t^2-4)/(t-2)];

(%o17) r(t):=[sin(t-2)^2/(t-2),(t^2-4)/(t-2)]

(%i18) limit(r(t),t,2);
 (%o18) [0,4]

̂ , then find derivative of 𝒓̅(t)

Q. If 𝒓̅(t) =cost 𝒊̂ – cost 𝒋̂ + sint 𝒌

(%i19) r(t):=[cos(t),-cos(t),sin(t)];

(%o19) r(t):=[cos(t),-cos(t),sin(t)]

(%i20) diff(r(t),t);

(%o20) [-sin(t),sin(t),cos(t)]

̂ , then find derivative of 𝒓̅(t), Also find derivative value at

Q. If 𝒓̅(t) =t3 𝒊̂ + sin(𝝅𝒕/𝟐) 𝒋̂ +log(t+2) 𝒌
t=0, 1, 𝝅

(%i35) r(t):=[t^3,sin(%pi*t/2),log(t+2)];

(%o35) r(t):=[t^3,sin((%pi*t)/2),log(t+2)]

(%i36) diff(r(t),t);

(%o36) [3*t^2,(%pi*cos((%pi*t)/2))/2,1/(t+2)]

(%i37) r(t):=[t^3,sin(%pi*t/2),log(t+2)];

(%o37) r(t):=[t^3,sin((%pi*t)/2),log(t+2)]

(%i38) define(rp(t),diff(r(t),t));

(%o38) rp(t):=[3*t^2,(%pi*cos((%pi*t)/2))/2,1/(t+2)]

(%i39) float(rp(1));

(%o39) [3.0,0.0,0.3333333333333333]

(%i40) float(%);

(%o40) [3.0,0.0,0.3333333333333333]

(%i41) float(rp(%pi/2));

(%o41) [7.402203300817018,-1.227124770575646,0.2800495767557787]

(%i42) float(rp(0));

(%o42) [0.0,1.570796326794897,0.5]

̂ , then find ∫ 𝒓̅(𝒕)dt

Q. If 𝒓̅(t) =cost 𝒊̂ + (1+sect tant) 𝒋̂ + sec2t 𝒌

(%i43) r(t):=[cos(t),1+sec(t)*tan(t),(sec(t))^2];

(%o43) r(t):=[cos(t),1+sec(t)*tan(t),sec(t)^2]

(%i44) integrate(r(t),t);

(%o44) [sin(t),1/cos(t)+t,tan(t)]
 Practical NO. 12

̂ , then find ∫𝟏 𝒓̅(𝒕)

Q. If 𝒓̅(t) =sect tant 𝒊̂ + tant 𝒋̂ + (2 sint cost) 𝒌 𝟎

(%i45) r(t):=[sec(t)*tan(t),tan(t),2*sin(t)*cos(t)];

(%o45) r(t):=[sec(t)*tan(t),tan(t),2*sin(t)*cos(t)]

(%i46) integrate(r(t),t,0,%pi/3);

(%o46) [1,log(2),3/4]

Q. Find T, N and the curvature 𝜿 for the plane curve r(t) = (2 cost, 2 sint)

(%i47) load(eigen);

(%o47) "C:/maxima-5.45.1/share/maxima/5.45.1/share/matrix/eigen.mac"

(%i48) r(t):=[2*cos(t),2*sin(t)];

(%o48) r(t):=[2*cos(t),2*sin(t)]

(%i49) define(rp(t),diff(r(t),t));

(%o49) rp(t):=[-2*sin(t),2*cos(t)]

(%i50) uvect(rp(t));

(%o50) [-(2*sin(t))/sqrt(4*sin(t)^2+4*cos(t)^2),(2*cos(t))/sqrt(4*sin(t)^2+4*cos(t)^2)]

(%i51) trigsimp(%);

(%o51) [-sin(t),cos(t)]

Therefore the unit tangent vector is obtain as follows

(%i52) define(T(t),%);

(%o52) T(t):=[-sin(t),cos(t)]

(%i54) define(Tp(t),diff(T(t),t));

(%o54) Tp(t):=[-cos(t),-sin(t)]

(%i57) uvect(Tp(t));

(%o57) [-cos(t)/sqrt(sin(t)^2+cos(t)^2),-sin(t)/sqrt(sin(t)^2+cos(t)^2)]

(%i58) trigsimp(%);

(%o58) [-cos(t),-sin(t)]

Therefore the unit normal vector is obtain as follows

(%i59) define(N(t),%);
 (%o59) N(t):=[-cos(t),-sin(t)]

(%i60) k:sqrt(innerproduct(Tp(t),Tp(t))/innerproduct(rp(t),rp(t)));

(%o60) sqrt(sin(t)^2+cos(t)^2)/sqrt(4*sin(t)^2+4*cos(t)^2)

(%i61) trigsimp(%);

(%o61) ½

Q. Find T, N and the curvature 𝜿 for the plane curve r(t) = (3sint, 3cost, 4t)

(%i1) load(eigen);

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/matrix/eigen.mac"

(%i2) r(t):=[3*sin(t),3*cos(t),4*t];

(%o2) r(t):=[3*sin(t),3*cos(t),4*t]

(%i3) define(rp(t),diff(r(t),t));

(%o3) rp(t):=[3*cos(t),-3*sin(t),4]

(%i4) uvect(rp(t));

(%o4) [(3*cos(t))/sqrt(9*sin(t)^2+9*cos(t)^2+16),-
(3*sin(t))/sqrt(9*sin(t)^2+9*cos(t)^2+16),4/sqrt(9*sin(t)^2+9*cos(t)^2+16)]

(%i5) trigsimp(%);

(%o5) [(3*cos(t))/5,-(3*sin(t))/5,4/5]

Therefore the unit tangent vector is obtain as follows

(%i6) define(T(t),%);

(%o6) T(t):=[(3*cos(t))/5,-(3*sin(t))/5,4/5]

(%i8) define(Tp(t),diff(T(t),t));

(%o8) Tp(t):=[-(3*sin(t))/5,-(3*cos(t))/5,0]

(%i9) uvect(Tp(t));

(%o9) [-(3*sin(t))/sqrt(9*sin(t)^2+9*cos(t)^2),-(3*cos(t))/sqrt(9*sin(t)^2+9*cos(t)^2),0]

(%i10) trigsimp(%);

(%o10) [-sin(t),-cos(t),0]

Therefore the unit normal vector is obtain as follows

(%i11) define(N(t),%);

(%o11) N(t):=[-sin(t),-cos(t),0]

(%i15) k:sqrt(innerproduct(Tp(t),Tp(t))/innerproduct(rp(t),rp(t)));
 (%o15) sqrt(9*sin(t)^2+9*cos(t)^2)/(5*sqrt(9*sin(t)^2+9*cos(t)^2+16))

(%i16) trigsimp(%);

(%o16) 3/25

Q. Find T, N and the curvature 𝜿 for the plane curve r(t) = (6 sin2t, 6 cos2t, 5t)

(%i1) load(eigen);

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/matrix/eigen.mac"

(%i17) r(t):=[6*sin(2*t),6*cos(2*t),5*t];

(%o17) r(t):=[6*sin(2*t),6*cos(2*t),5*t]

(%i18) define(rp(t),diff(r(t),t));

(%o18) rp(t):=[12*cos(2*t),-12*sin(2*t),5]

(%i19) uvect(rp(t));

(%o19) [(12*cos(2*t))/sqrt(144*sin(2*t)^2+144*cos(2*t)^2+25),-
(12*sin(2*t))/sqrt(144*sin(2*t)^2+144*cos(2*t)^2+25),5/sqrt(144*sin(2*t)^2+144*cos(2*t)^2+25)]

(%i20) trigsimp(%);

(%o20) [(12*cos(2*t))/13,-(12*sin(2*t))/13,5/13]

Therefore the unit tangent vector is obtain as follows

(%i21) define(T(t),%);

(%o21) T(t):=[(12*cos(2*t))/13,-(12*sin(2*t))/13,5/13]

(%i22) define(Tp(t),diff(T(t),t));

(%o22) Tp(t):=[-(24*sin(2*t))/13,-(24*cos(2*t))/13,0]

(%i23) uvect(Tp(t));

(%o23) [-(24*sin(2*t))/sqrt(576*sin(2*t)^2+576*cos(2*t)^2),-
(24*cos(2*t))/sqrt(576*sin(2*t)^2+576*cos(2*t)^2),0]

(%i24) trigsimp(%);

(%o24) [-sin(2*t),-cos(2*t),0]

Therefore the unit normal vector is obtain as follows

(%i25) define(N(t),%);

(%o25) N(t):=[-sin(2*t),-cos(2*t),0]

(%i26) k:sqrt(innerproduct(Tp(t),Tp(t))/innerproduct(rp(t),rp(t)));

(%o26) sqrt(576*sin(2*t)^2+576*cos(2*t)^2)/(13*sqrt(144*sin(2*t)^2+144*cos(2*t)^2+25))
(%i27) trigsimp(%);

(%o27) 24/169

̂ from t=0 to t=𝝅. Also find

Q. Find the arc length along the curve 𝒓̅(t) = 2 cost𝒊̂ + 2 sint 𝒋̂ + √𝟓t 𝒌
the unit tangenst vector of the curve.

(%i1) load(eigen);

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/matrix/eigen.mac"

(%i2) r(t):=[2*cos(t),2*sin(t),sqrt(5)*t];

(%o2) r(t):=[2*cos(t),2*sin(t),sqrt(5)*t]

(%i3) define(rp(t),diff(r(t),t));

(%o3) rp(t):=[-2*sin(t),2*cos(t),sqrt(5)]

(%i4) uvect(rp(t));

(%o4) [-
(2*sin(t))/sqrt(4*sin(t)^2+4*cos(t)^2+5),(2*cos(t))/sqrt(4*sin(t)^2+4*cos(t)^2+5),sqrt(5)/sqrt(4*sin(t
)^2+4*cos(t)^2+5)]

(%i5) trigsimp(%);

(%o5) [-(2*sin(t))/3,(2*cos(t))/3,sqrt(5)/3]

Therefore the unit tangent vector is obtain as follows

(%i6) define(T(t),%);

(%o6) T(t):=[-(2*sin(t))/3,(2*cos(t))/3,sqrt(5)/3]

(%i13) define(Tp(t),diff(T(t),t));

(%o13) Tp(t):=[-(2*cos(t))/3,-(2*sin(t))/3,0]

(%i15) k:sqrt(Tp(t).Tp(t))/sqrt(rp(t).rp(t));

(%o15) sqrt((4*sin(t)^2)/9+(4*cos(t)^2)/9)/sqrt(4*sin(t)^2+4*cos(t)^2+5)

Curvature is as follows

(%i16) trigsimp(%);

(%o16) 2/9

Arc length is calculated as follows

(%i19) L:integrate(trigsimp(sqrt(rp(t).rp(t))),t,0,%pi);

(%o19) 3*%pi

Q. Find the curvature arc length along the curve 𝑟̅ (t) = 3 cost𝑖̂ + 3 sint 𝑗̂ + t 𝑘̂ from t=0 to t=1.
(%i30) load(eigen);

(%o30) "C:/maxima-5.45.1/share/maxima/5.45.1/share/matrix/eigen.mac"

(%i31) r(t):=[3*cos(t),3*sin(t),t];

(%o31) r(t):=[3*cos(t),3*sin(t),t]

(%i32) define(rp(t),diff(r(t),t));

(%o32) rp(t):=[-3*sin(t),3*cos(t),1]

(%i33) uvect(rp(t));

(%o33) [-
(3*sin(t))/sqrt(9*sin(t)^2+9*cos(t)^2+1),(3*cos(t))/sqrt(9*sin(t)^2+9*cos(t)^2+1),1/sqrt(9*sin(t)^2+9
*cos(t)^2+1)]

(%i34) trigsimp(%);

(%o34) [-(3*sin(t))/sqrt(10),(3*cos(t))/sqrt(10),1/sqrt(10)]

The unit tangent vector of the curve is

(%i35) define(T(t),%);

(%o35) T(t):=[-(3*sin(t))/sqrt(10),(3*cos(t))/sqrt(10),1/sqrt(10)]

(%i36) define(Tp(t),diff(T(t),t));

(%o36) Tp(t):=[-(3*cos(t))/sqrt(10),-(3*sin(t))/sqrt(10),0]

(%i37) k:sqrt(Tp(t).Tp(t))/sqrt(rp(t).rp(t));

(%o37) sqrt((9*sin(t)^2)/10+(9*cos(t)^2)/10)/sqrt(9*sin(t)^2+9*cos(t)^2+1)

(%i38) trigsimp(%);

(%o38) 3/10

Arc length is calculated as follows

(%i40) L:integrate(trigsimp(sqrt(rp(t).rp(t))),t,0,1);

(%o40) sqrt(10)

Q. Draw the graph of the function f(x,y)=(x2-y2)2

(%i1) f(x,y):=(x^2-y^2)^2;

(%o1) f(x,y):=(x^2-y^2)^2

(%i2) load(draw);

(%o2) "C:/maxima-5.45.1/share/maxima/5.45.1/share/draw/draw.lisp"

(%i3) draw3d(explicit(f(x,y),x,-3,3,y,-3,3));
(%i3) draw3d(enhanced3d=true,explicit(f(x,y),x,-3,3,y,-3,3));

Q. Find the gradient of the scalar function f(x,y)=(x2-y2)

(%i5) load("vect");

(%o5) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i6) gdf(x,y):=grad(f(x,y));

(%o6) gdf(x,y):=grad(f(x,y))

(%i10) f(x,y):=x^2-y^2;

(%o10) f(x,y):=x^2-y^2

(%i11) ev(express(gdf(x,y)),diff);

(%o11) [2*x,-2*y,0]

(%i12) define(gdf(x,y),%);

(%o12) gdf(x,y):=[2*x,-2*y,0]

Q. Find the gradient vector field of the function f(x,y)=x2sin5y

(%i24) load(vect);

"vect: warning: removing existing rule or rules for "."."" "

(%o24) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i25) gdf(x,y):=grad(f(x,y));

(%o25) gdf(x,y):=grad(f(x,y))
(%i26) f(x,y):=x^2*sin(5*y);

(%o26) f(x,y):=x^2*sin(5*y)

(%i27) ev(express(gdf(x,y)),diff);

(%o27) [2*x*sin(5*y),5*x^2*cos(5*y),0]

(%i28) define(gdf(x,y),%);

(%o28) gdf(x,y):=[2*x*sin(5*y),5*x^2*cos(5*y),0]

Q. Find the line integral of f(x,y)=yex2 along the curve 𝒓̅(t) = 4 t𝒊̂ - 3t 𝒋̂, -1≤ t≤ 2

(%i1) f(x,y):=y*exp(x^2);

(%o1) f(x,y):=y*exp(x^2)

(%i2) [x,y]:[4*t,-3*t];

(%o2) [4*t,-3*t]

(%i3) rp:diff([x,y],t);

(%o3) [4,-3]

(%i5) integrate(f(x,y)*sqrt(rp.rp),t,-1,2);

(%o5) -15*(%e^64/32-%e^16/32)

(%i6) float(%);

(%o6) -2.922726131630445*10^27

Q. Evaluate ∫𝒄 𝑿ds , where C is the straight line segment x=t, y=t/2, from (0,0) to (4,2)

(%i7) f(x,y):=x;

(%o7) f(x,y):=x

(%i8) [x,y]:[t,t/2];

(%o8) [t,t/2]

(%i9) rp:diff([x,y],t);

(%o9) [1,1/2]

(%i10) integrate(f(x,y)*sqrt(rp.rp),t,0,4);

(%o10) 4*sqrt(5)

(%i11) float(%);

(%o11) 8.94427190999916
 ̂
Q. Evaluate ∫𝒄 𝑭d, where F(x,y,z)= z 𝒊̂ + xy 𝒋̂ -y2 𝒌 along the curve C given by 𝒓̅(t) = t2 𝒊̂ + t 𝒋̂ -
̂
√𝒕 𝒌

, where 0 ≤ t ≤ 1

(%i12) F(x,y,z):=[z,x*y,-y^2];

(%o12) F(x,y,z):=[z,x*y,-y^2]

(%i13) [x,y,z]:[t^2,t,sqrt(t)];

(%o13) [t^2,t,sqrt(t)]

(%i15) romberg(F(x,y,z).diff([x,y,z],t),t,0,1);

(%o15) 0.8500022849339077

̂
Q. Evaluate ∫𝒄 −𝒚𝒅𝒙 + 𝒛𝒅𝒚 + 𝟐𝒙𝒅𝒛, where C is the helix 𝒓̅(t) = cost𝒊̂ + sint 𝒋̂ + t 𝒌

And 0 ≤ t ≤ 2𝝅

(%i16) F(x,y,z):=[-y,z,2*x];

(%o16) F(x,y,z):=[-y,z,2*x]

(%i17) [x,y,z]:[cos(t),sin(t),t];

(%o17) [cos(t),sin(t),t]

(%i19) romberg(F(x,y,z).diff([x,y,z],t),t,0,2*%pi);

(%o19) 3.141592657837756

Q. Find whether the vector field F(X,Y) =(x3-4y2)𝒊̂ + (4y3 -3x) 𝒋̂ is conservative or not.

(%i1) load("vect");

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i2) F(x,y):=[x^3-4*y^2,4*y^3-3*x];

(%o2) F(x,y):=[x^3-4*y^2,4*y^3-3*x]

(%i3) scalefactors([x,y]);

(%o3) done

(%i4) curl(F(x,y));

(%o4) curl([x^3-4*y^2,4*y^3-3*x])

(%i5) express(%);

(%o5) 'diff((4*y^3-3*x),x,1)-'diff((x^3-4*y^2),y,1)

(%i6) ev(%,diff);
 (%o6) 8*y-3

The curl is not zero , hence field is not conservative.

̂ is conservative or not.
Q. Find whether the vector field F(X,Y) =(z + y) 𝒊̂ + (z) 𝒋̂ +(y + x) 𝒌

(%i1) load("vect");

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i25) F(x,y,z):=[z+y,z,y+x];

(%o25) F(x,y,z):=[z+y,z,y+x]

(%i26) scalefactors([x,y,z]);

(%o26) done

(%i27) curl(F(x,y,z));

(%o27) curl([z+y,z,y+x])

(%i28) express(%);

(%o28) ['diff((y+x),y,1)-1,'diff((z+y),z,1)-'diff((y+x),x,1),'diff(z,x,1)-'diff((z+y),y,1)]

(%i29) ev(%,diff);

(%o29) [0,0,-1]

The curl is not zero , hence field is not conservative.

̂ is conservative or not.
Q. Find whether the vector field F(X,Y) =yz 𝒊̂ + xz 𝒋̂ + xy 𝒌

(%i1) load("vect");

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i20) F(x,y,z):=[y*z,x*z,x*y];

(%o20) F(x,y,z):=[y*z,x*z,x*y]

(%i21) scalefactors([x,y,z]);

(%o21) done

(%i22) curl(F(x,y,z));

(%o22) curl([y*z,x*z,x*y])

(%i23) express(%);

(%o23) ['diff((x*y),y,1)-'diff((x*z),z,1),'diff((y*z),z,1)-'diff((x*y),x,1),'diff((x*z),x,1)-'diff((y*z),y,1)]

(%i24) ev(%,diff);

(%o24) [0,0,0]
The curl is zero , hence field is conservative. Further we can find the scalar potential f in Maxima
with the potential function.

(%i1) load("vect");

(%o1) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i30) F(u,v,w):=[v*w,u*w,u*v];

(%o30) F(u,v,w):=[v*w,u*w,u*v]

(%i32) scalefactors([u,v,w]);

(%o32) done

(%i33) potential(F(u,v,w));

(%o33) u*v*w

̂
Q. Find a potential function f for the field F = 2x 𝒊̂ + 3y 𝒋̂ + 4z 𝒌

(%i2) load("vect");

(%o2) "C:/maxima-5.45.1/share/maxima/5.45.1/share/vector/vect.mac"

(%i3) F(x,y,z):=[2*x,3*y,4*z];

(%o3) F(x,y,z):=[2*x,3*y,4*z]

(%i4) scalefactors([x,y,z]);

(%o4) done

(%i5) curl(F(x,y,z));

(%o5) curl([2*x,3*y,4*z])

(%i6) express(%);

(%o6) ['diff((4*z),y,1)-'diff((3*y),z,1),'diff((2*x),z,1)-'diff((4*z),x,1),'diff((3*y),x,1)-'diff((2*x),y,1)]

(%i7) ev(%,diff);

(%o7) [0,0,0]

(%i8) potential(F(x,y,z));

(%o8) (4*z^2+3*y^2+2*x^2)/2