Project-I-sol.mw

restart; 1 

with(plots); 1 

[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
[Interactive, animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, cylinderplot, densityplot, ...
 

with(linalg); 1 

[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, cold...
 

interface(showassumed = 0); 1 

1 

assume(t::real); 1 

 

lowt := -2; 1 

-2 

hight := 2; 1 

 

r := [cos(t), sin(t), t]; 1 

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

v := diff(r, t); 1 

[-sin(t), cos(t), 1] 

s := sqrt(innerprod(v, v)); 1 

(1+sin(t)^2+cos(t)^2)^(1/2) 

s := simplify(s); 1 

2^(1/2) 

T := v/s; 1 

1/2*([-sin(t), cos(t), 1])*2^(1/2) 

N := (diff(T, t))/simplify(sqrt(innerprod(diff(T, t), diff(T, t)))); 1 

(1/2*([-cos(t), -sin(t), 0])*2^(1/2)+1/8*([0, 0, 0])*2^(1/2))*2^(1/2) 

 

 

 

 

 

rf := unapply(r, t); 1 

proc (t) options operator, arrow; [cos(t), sin(t), t] end proc 

Tf := unapply(T, t); 1 

proc (t) options operator, arrow; 1/2*([-sin(t), cos(t), 1])*2^(1/2) end proc 

Nf := unapply(N, t); 1 

proc (t) options operator, arrow; (1/2*([-cos(t), -sin(t), 0])*2^(1/2)+1/8*([0, 0, 0])*2^(1/2))*2^(1/2) end proc 

c1 := spacecurve(convert(r, list), t = lowt .. hight, axes = normal, thickness = 3); 1
c1 := spacecurve(convert(r, list), t = lowt .. hight, axes = normal, thickness = 3); 1
 

 

Plot 

 

nframes := 20; 1 

20 

deltat := (hight-lowt)/nframes; 1 

1/5 

 

 

trips := seq(display(arrow(rf(lowt+deltat*i), evalf(Tf(lowt+deltat*i)), shape = arrow, color = black), arrow(rf(lowt+deltat*i), evalf(Nf(lowt+deltat*i)), shape = arrow, color = red), thickness = 3), i...
trips := seq(display(arrow(rf(lowt+deltat*i), evalf(Tf(lowt+deltat*i)), shape = arrow, color = black), arrow(rf(lowt+deltat*i), evalf(Nf(lowt+deltat*i)), shape = arrow, color = red), thickness = 3), i...
trips := seq(display(arrow(rf(lowt+deltat*i), evalf(Tf(lowt+deltat*i)), shape = arrow, color = black), arrow(rf(lowt+deltat*i), evalf(Nf(lowt+deltat*i)), shape = arrow, color = red), thickness = 3), i...
 

INTERFACE_PLOT3D(CURVES([[-.4161468365, -.9092974268, -2.], [.2268235400, -1.203557677, -1.292893219]], [[.1124654960, -1.182967907, -1.463182089], [.2268235400, -1.203557677, -1.292893219], [0.839934...
INTERFACE_PLOT3D(CURVES([[-.4161468365, -.9092974268, -2.], [.2268235400, -1.203557677, -1.292893219]], [[.1124654960, -1.182967907, -1.463182089], [.2268235400, -1.203557677, -1.292893219], [0.839934...
INTERFACE_PLOT3D(CURVES([[-.4161468365, -.9092974268, -2.], [.2268235400, -1.203557677, -1.292893219]], [[.1124654960, -1.182967907, -1.463182089], [.2268235400, -1.203557677, -1.292893219], [0.839934...
INTERFACE_PLOT3D(CURVES([[-.4161468365, -.9092974268, -2.], [.2268235400, -1.203557677, -1.292893219]], [[.1124654960, -1.182967907, -1.463182089], [.2268235400, -1.203557677, -1.292893219], [0.839934...
INTERFACE_PLOT3D(CURVES([[-.4161468365, -.9092974268, -2.], [.2268235400, -1.203557677, -1.292893219]], [[.1124654960, -1.182967907, -1.463182089], [.2268235400, -1.203557677, -1.292893219], [0.839934...
 

 

d1 := display(trips, insequence = true); 1 

 

Plot 

 

 

display(c1, d1, scaling = constrained); 1 

 

Plot