Vector Valued Functions

Transcription

1 Vector Valued Functions Introduction and Goals: The main goal of this lab will help you visualie the graphs of vector-valued functions and their tangent vectors. Moreover, we want you to begin to view tangent vectors as the three-dimensional equivalent to the slope of the tangent line from Calculus I. Recall that the slope of a tangent line told us the rate of change of one variable with respect to another. Specifically, if we were eamining a distance verses time graph then the slope of the tangent line gives us the instantaneous rate of change of distance with respect to time, that is, velocity. In three dimensions, if we are a particle moving along a curve then the tangent vector gives us the direction we are moving in and the length of that tangent vector gives us our speed in which we are moving along the curve. This lab also investigates the Maple commands needed to perform some of the calculations that we would do by hand, in class. Furthermore, we continue to build an arsenal of special commands to make using Maple a little easier. Before You Start: Make sure that you read and understand the mathematics from the corresponding sections in your tetbook. It would also help if you looked over the lab on Lines and Planes. Tetbook Correspondence: Stewart 5 th Edition: Stewart 5 th Edition Early Transcendentals: Thomas Calculus 0 th Edition Early Transcendentals: 0.5. Johnston & Mathews Calculus: 8.7. Maple Commands and Packages Used: Packages: plots, linalg, VectorCalculus. Commands: spacecurve, plot3d, display, solve, limit, diff, unapply, arrow, int, solve, fsolve, evalf, seq. User Defined Commands: plotvvf, plotvvf, plotvvftv, normplot, normplotresy, plotvvftv, curveint. History & Biographies: Maple Commands:

2 Defining and Plotting Space Curves: As usual, we will need the plots package and there are a few commands we are using that require the linalg package, so load them in with the with command. > with(plots): > with(linalg): We are also going to take a more sophisticated look at vectors. In previous labs we were using matrices to denote vectors. We mainly used single row matrices to denote the vector but we also looked at a single column matri representation of the vector. Using matrices for vectors is very natural since this is the way we represent them mathematically and it is a good way to get your feet wet with vectors and Maple but in general we want to do many things with vectors that we would not necessarily want to do with a matri. For eample, differentiate and integrate them componentwise as well as create vector-valued functions from them. Although you can do most vector manipulations with the matri representation some of the Maple commands we would need are a bit cumbersome. To make our Maple commands a little easier and load in several interesting function we will add the VectorCalculus package. > with(vectorcalculus): Warning, the assigned names <,> and < > now have a global binding Warning, these protected names have been redefined and unprotected: *, +,., Vector, diff, int, limit, series As you know from the tet, a space curve is given either parametrically as where f ( t), g ( t) and ( t) r ( t ) = f ( t), g( t), h( t) = y = = f ( t) g( t) h( t) h are functions of a single variable t, or in vector form as. When the VectorCalculus package is loaded we define a vector as a component list, as usual, but this time we use angled brackets. For eample, say we wanted to input the vector sin ( t ),cos( t), t, the Maple synta would be, > <sin(t),cos(t),sqrt(t)>; sin( t ) e + cos( t) e + y t e Note that the output is different than the angled bracket output from previous labs. When the VectorCalculus package is not loaded the angled brackets gave us a column vector, actually a single column matri. With the VectorCalculus package loaded the angled

3 brackets automatically see the epression as a vector, i.e. a vector data type in Maple. The e, e y and e stand for the vectors i, j and k respectively, that is,0, 0, 0,, 0 and 0,0,. We can define this as the vector valued function r ( t ) = sin( t),cos( t), t simply by placing the function definition synta at the beginning. That is, > r:=t-><sin(t),cos(t),sqrt(t)>; r := t VectorCalculus:-`<,>` ( sin( t ), cos( t ), t ) Note that the output is a bit strange. Maple seems to be calling a special function from the VectorCalculus package. If we evaluate r(t) we see that all is well. > r(t); sin( t ) e + cos( t) e + y t e Similarly, we can evaluate the vector-valued function using standard function notation. > r(pi); e + y π e > r(pi/); e + π e > r(0); e y Let s define three more vector functions so that we have something to work with for the s ( t) = cos t,sin t,sin t cos t, remainder of the lab. Define ( ) ( ) ( ) ( ) 3 q ( t ) = t cos( t), t sin( t), t and p ( t ) = t, t, t, > s:=t-><cos(t),sin(t),sin(t)^-cos(t)^>; s := t VectorCalculus:-`<,>` ( cos( t ), sin( t ), VectorCalculus:-`+` ( sin( t ), VectorCalculus:-`*` ( cos( t ), -) )) > s(t); cos( t ) e + sin( t) e + ( sin( t) cos( t) ) e y > q:=t-><t*cos(t),t*sin(t),t>;

4 q := t VectorCalculus:-`<,>` ( VectorCalculus:-`*` ( t, cos( t ) ), VectorCalculus:-`*` ( t, sin( t ) ), t) > q(t); t cos( t ) e + t sin( t ) e + t e y > p:=t-><t,t^,t^3>; p := t VectorCalculus:-`<,>` ( t, t, t 3 ) > p(t); t e + t e + t 3 e y The standard way to plot a vector-valued function is with the spacecurve command. Unfortunately, using the spacecurve command directly does not work because the spacecurve command is epecting a different form for the function. > spacecurve(r(t),t=0..0*pi,numpoints=500); Error, (in spacecurve) first argument must be array, list or set To appease the spacecurve command we need to strip the components from r and place them in matri form. This can be done with the seq command, as follows, > spacecurve([seq(r(t)[n],n=..3)],t=0..0*pi,numpoints=500, aes=boxed); Note that one option we included in this command was numpoints, recall that numpoints is the number of points plotted to draw the graph. It is common with space curves that the default number of points is too small and the curve looks polygonal. Now this command is quite long and we really don t want to type it in repeatedly, of even copy and past it repeatedly, so we are going to define a function that will make graphing space curves quick and easy. In the command above there are four things of interest, the

5 function r (t), the lower bound for t, the upper bound for t and the number of points. When we create the new function these four items will be represented by variables. > plotvvf:=(r,lb,ub,pts)->spacecurve([seq(r(t)[n],n=..3)], t=lb..ub,numpoints=pts); plotvvf := ( r, lb, ub, pts ) spacecurve ([ seq ( r( t ), n =.. 3 )], t = lb.. ub, numpoints = pts ) n > plotvvf(r,0,0*pi,500); which is the same image as above. The nifty thing with the command we just defined is that we can substitute any vector-valued function into the command and get a graph. > plotvvf(s,0,*pi,500); > plotvvf(q,-0*pi,0*pi,500);

6 > plotvvf(p,0,,500); We can also combine the space curve images together using the display command. > display(plotvvf(s,0,*pi,500),plotvvf(q,-pi,pi,500));

7 Note that in the above image we selected : and the boed aes options from the toolbar after the graph was drawn. If you want to add more options to your plotting command simply place a variable name in the list of variables and place the corresponding options in the spacecurve command. For eample, the following adds color and aes type to the list of parameters. > plotvvf:=(r,lb,ub,pts,col,as)->spacecurve([seq(r(t)[n], n=..3)],t=lb..ub,numpoints=pts,color=col,aes=as); plotvvf := ( r, lb, ub, pts, col, as ) spacecurve ( [ seq ( r( t), n =.. 3 ) ], t = lb.. ub, n numpoints = pts, color = col, aes = as ) > plotvvf(s,-pi,pi,500,black,boxed); > plotvvf(s,-pi,pi,00,red,normal);

8 Limits of Vector Valued Functions: Finding limits of vector-valued functions in Maple is as easy as finding limits of real valued function in Maple. Use the limit command and give the command the vectorvalued function, the limit point and an optional direction. To see some of the interesting h ( t) = t,tan t, ln t, types of output let s use a slightly nastier function. Define ( ) ( ) > h:=t-></t,tan(t),ln(t)>; h := t VectorCalculus:-`<,>` VectorCalculus:-`*`,, tan( t ), ln( t ) t > h(t); t e + tan( t ) e + y ln( t) e If we evaluate the limit of this function at a point where the function is continuous we see that Maple simply evaluates the function at the limit point, as it should. > limit(h(t),t=5); > limit(h(t),t=pi); 5 e + tan( 5 ) e + y π e + ln( 5 ) ln( π) e e If we evaluate the limit at 0 we see that the component is undefined, the y component is 0 and the component is negative infinity, as we would epect. Although the limit vector

9 itself is of little use it does inform us to the limit of each of the components of the function, which is really all we can epect. > limit(h(t),t=0); undefined e e If we further find the left and right hand limits of the function we get the corresponding outputs. > limit(h(t),t=0,left); > limit(h(t),t=0,right); e e e e The limits at π give similar output. > limit(h(t),t=pi/); + + π e undefined e ( ln( ) + ln( π) ) e y > limit(h(t),t=pi/,left); + + π e e ( ln( ) + ln( π) ) e y > limit(h(t),t=pi/,right); + π e e ( ln( ) + ln( π) ) e y Finally, as with real valued functions we can also find limits at infinity. > limit(h(t),t=infinity); undefined e + e y Derivatives of Vector Valued Functions and Tangent Vectors: Differentiation of vector-valued functions can be done with the diff command just like real valued functions. In fact, the same options apply; the only difference is that the output is in vector form. > diff(r(t),t);

10 cos( t ) e sin( t) e + y t e > diff(s(t),t); > diff(p(t),t); sin( t ) e + cos( t ) e + y e + t e + 3 t e y 4 sin( t ) cos( t) e > diff(p(t),t,t); e + 6 t e y > diff(p(t),t,t,t); 6 e > diff(p(t),t,t,t,t); 0 e > diff(p(t),t$3); 6 e > diff(p(t),t$4); 0 e Frequently we wish to define a function that is the derivative of another function. To do this we invoke the unapply command on the derivative. For eample, to define the function dr as the first derivative of r we use the command, > dr:=unapply(diff(r(t),t),t); dr := t rtable.. 3, { ( ) = cos( t ), ( ) = sin( t ), ( 3 ) = }, datatype = anything, t subtype = Vector column, storage = rectangular, order = Fortran_order, attributes = [ coords = cartesian ] > dr(t); cos( t ) e sin( t) e + y t e Now to find the derivative, which is the tangent vector, at any value of the parameter t we simply evaluate the derivative function at the desired value. For eample, the following are the tangent vectors to r at, π and 0 respectively.

11 > dr(); > dr(pi); cos( ) e sin( ) e + y e + π e e > dr(0); Error, (in dr) numeric eception: division by ero With the plotvvf command we created above and the arrow command we can create some very informative graphs of space curves. Recall that the command arrow(u,v); draws the vector v beginning at the point u. So if we let u be the point on the curve corresponding to a particular parameter value and v be the tangent vector at the same value for the parameter we can graph the tangent vector directly on the curve. For eample, to plot the curve r along with its tangent vector at t = π we would use the command, > display(plotvvf(r,0,0*pi,500,black,boed), arrow(r(pi),dr(pi),shape=arrow,color=black)); 3π Similarly, to plot the curve r along with its tangent vector at t = we would use the command, > display(plotvvf(r,0,0*pi,500,black,boed), arrow(r(3*pi/),dr(3*pi/),shape=arrow,color=black));

12 If we go a little further and use the seq function on the arrow command we can plot a series of tangent vectors along the curve. For eample, consider the following command, > display(plotvvf(r,0,0*pi,500,black,boed), seq(arrow(r(0*pi/40*n),dr(0*pi/40*n),shape=arrow, color=black),n=..40)); This is a long command but it is actually easy to understand if we break it down into manageable chunks. First the plotvvf portion plotvvf(r,0,0*pi,500,black,boed) simply plots the curve r from 0 to 0 π using 500 points in the plot, drawing the curve in black and putting boed aes around it. The seq portion is a little more involved, seq(arrow(r(0*pi/40*n),dr(0*pi/40*n), shape=arrow,color=black),n=..40)

13 In the arrow portion, the part r(0*pi/40*n),dr(0*pi/40*n) plots the tangent vector at 0*Pi/40*n for each n from to 40 starting at the point r(0*pi/40*n). That is, we are drawing 40 tangent vectors on the curve at t values of π π 4, 3 π π 5 π, 3 π 7 π 4,, π 9 π 5 π π 4, 3 π 3 π 7 π, 5 π 4 π 7 π 4,, 9 π 4,, 4,, 4,, 4,, 4,, 9 π 5 π π π 3 π 6 π 5 π 3 π 7 π 7 π 9 π 5 π 3 π,,,,,,,,,,,,, 8 π, π 7 π 35 π 9 π 37 π 9 π 39 π,,,,,,, 0 π As you can see from the above image this plot shows us the tangent vector structure along the curve. Hence it tells us the direction and magnitude of the rate of change, with respect to t, along the curve. So if we think about t as time these vectors are showing us the direction and speed of a particle that is traveling along the curve. Since this type of graph can be really useful we would like to streamline the command a little. It would be nice if we could create a command like plotvvf where we simply put in the vector-valued function and some other options and Maple would create the image. We will do this in two steps. First note that we can replace dr by unapply(diff(r(t),t),t), this will relieve us from having to find the derivative before we do the plot. Note that the following command produces the same image but dr is never used. > display(plotvvf(r,0,0*pi,500,black,boed), seq(arrow(r(0*pi/40*n),unapply(diff(r(t),t),t)(0*pi/40*n),shape=arrow,color=black),n=..40));

14 Now we are ready to create the command, we will call it plotvvftv for plot vector-valued function and tangent vectors. In the definition below note that r, lb, ub, pts, col and as are the same as in the plotvvf command we created earlier. We have added the parameters evf, nlb, nub, shp and vcol. Parameter evf nlb nub shp vcol Purpose This is the point evaluation function. That is, it is the formula, in terms of n, that is used to find the values of t where the tangent vector will be drawn. The formula must be in terms of n. The values that n takes on are all the integer values between nlb and nub. The lower bound for the integer range of n. The value must be an integer. The upper bound for the integer range of n. The value must be an integer. The shape of the arrows, harpoon, arrow, double_arrow, or cylindrical_arrow The color of the vectors. In the actual command you will notice that we have replaced plotvvf with the definition of plotvvf. This is so that we do not need to define the command plotvvf before defining this one. We also replaced all of the option values with their respective parameter names. Take some time to read through this command definition and really understand it. Maple is a system where if you take the time to create specialied functions like plotvvf and plotvvftv you can save yourself much more time in editing previous commands. > plotvvftv:=(r,lb,ub,pts,col,as,evf,nlb,nub,shp,vcol)-> display(spacecurve([seq(r(t)[n],n=..3)],t=lb..ub, numpoints=pts,color=col,aes=as),seq(arrow(r(evf), unapply(diff(r(t),t),t)(evf),shape=shp,color=vcol), n=nlb..nub)); plotvvftv := ( r, lb, ub, pts, col, as, evf, nlb, nub, shp, vcol ) display( spacecurve( [ seq ( r( t ), n =.. 3 ) ], t = lb.. ub, numpoints = pts, color = col, aes = as ), seq( n arrow ( r( evf ), unapply ( VectorCalculus:-diff ( r( t ), t ), t )( evf ), shape = shp, color = vcol ), n = nlb.. nub )) Now to create the sane plot as before we only need to input the following command. Take the time to match up the inputs below with the command definition and the last display command that produced this graph. > plotvvftv(r,0,0*pi,00,black,boed,0*pi/40*n,,40, arrow,black);

15 A few quick alterations to the calling command allow us to generate similar graphs of the other functions. > plotvvftv(s,0,*pi,00,black,boed,*pi/40*n,,40,arrow, black); > plotvvftv(q,-0*pi,0*pi,500,red,boed,-0*pi+0*pi/00*n, 0,00,arrow,black);

16 Zooming in on the previous image gives. > plotvvftv(q,0,3,50,blue,normal,3/0*n,0,0,arrow,black); > plotvvftv(p,0,3,50,blue,normal,3/0*n,0,0,arrow,black);

17 Integrating Vector Valued Functions: As with differentiation, integration of vector-valued functions is the same as with real valued functions. You can do both definite and indefinite integrals using the int command. For eample, the indefinite integral of p is, > int(p(t),t); t e t e y t 4 4 e and the definite integral from 0 to of p is, > int(p(t),t=0..); e e y 4 e Similarly, > int(q(t),t); ( cos( t) + t sin( t )) e + ( sin( t ) t cos( t) ) e + y t e > int(r(t),t); cos( t ) e + sin( t) e + y t ( 3/ ) 3 e

18 Also, as with differentiation, we can create functions of both definite and indefinite integrals by using the unapply command. For eample, the following command takes the indefinite integral of r and sets it to the function name ir. > ir:=unapply(int(r(t),t),t); ir := t rtable.. 3, { ( ) = cos( t ), ( ) = sin( t ), ( 3 ) = } 3 t( 3/ ), datatype = anything, subtype = Vector column, storage = rectangular, order = Fortran_order, attributes = [ coords = cartesian ] > ir(t); cos( t ) e + sin( t) e + y t ( 3/ ) 3 e Similarly, the following command takes the definite integral of r from a to b and sets it to the function name ir. > ir:=unapply(int(r(t),t=a..b),a,b); ir := ( a, b ) rtable.. 3, {( ) = cos( b ) + cos( a ), ( ) = sin( b ) sin( a ), ( 3 ) = 3/ } 3 a( / ) 3 b( ) datatype = anything, subtype = Vector column, storage = rectangular, order = Fortran_order, attributes = [ coords = cartesian ] > ir(a,b); ( cos( b ) + cos( a )) e + ( sin( b ) sin( a )) e + y b ( 3/ ) 3 a ( 3/ ) 3 3, e So finding the definite integral of r from 0 to can be done by, > ir(0,); ( cos( ) + ) e + sin( ) e + y 3 e Solving Equations & Finding Curve Intersections: Finding the point or points of intersection of two vector-valued functions can be difficult to do by hand and in fact even when using Maple the process can be less than straightforward. We will look at three eamples of increasing degrees of difficulty. An

19 easy eample is to find the points of intersection between the curves k ( t) = t, t t, t t and j ( t) = t, t t, t t. If we define and plot both of these curves together we see that there appears to be two points of intersection. > k:=t-><t^,t^3-t,t^4-t>; k := t VectorCalculus:-`<,>`( t, VectorCalculus:-`+` ( t 3, VectorCalculus:-`*` ( t, -)), VectorCalculus:-`+` ( t 4, VectorCalculus:-`*` ( t, -)) ) > j:=t-><t,t^-t,t^3-t>; j := t VectorCalculus:-`<,>`( t, VectorCalculus:-`+` ( t, VectorCalculus:-`*` ( t, -)), VectorCalculus:-`+` ( t 3, VectorCalculus:-`*` ( t, -) )) > display(plotvvf(k,-.5,.,00,red,boed),plotvvf(j,-.5,.,00,black,boed)); To find the values of t that correspond to the points of intersection we need to set the two curves, component by component, equal to each other and solve the resulting system of equations. To put the system in the correct format for the solve command we need to etract the first component of the first curve and set it equal to the first component of the second curve and then do the same for the other two components. Recall that to etract a component from a vector we simply need to place [n] after the vector, where n is the component to be etracted. So the solve command we need is > solve({k(t)[]=j(u)[],k(t)[]=j(u)[],k(t)[3]=j(u)[3]}, {t,u}); { u = 0, t = 0 }, { t =, u = } which can be streamlined by using the seq function to generate the list of equations. That is, > solve({seq(k(t)[n]=j(u)[n],n=..3)},{t,u}); { u = 0, t = 0 }, { t =, u = }

20 Checking the solutions shows that Maple has found the two points of intersection. > k(0); 0 e > j(0); 0 e > k(); e > j(); e Let s move to a moderately more difficult problem. Find the points of intersection s ( t) = cos t,sin t,sin t cos t (defined before) and between the curves ( ) ( ) ( ) ( ) h ( t ) = sin( t),cos( t), to be four points of intersection.. If we plot both of these curves together we see that there appears > h:=t-><sin(t),cos(t),/>; h := t VectorCalculus:-`<,>` sin( t ), cos( t ), VectorCalculus:-`*`, > display(plotvvf(h,0,*pi,500,black,boed), plotvvf(s,0,*pi,500,black,boed)); Setting up and evaluating the solve command for these curves gives a long and incomprehensible answer. Applying the evalf function to this result gives two solutions (note that the parameter values between the two curves is different.). Checking these

21 shows that that are legitimate solutions, with a minor bit of round-off error, which we will ignore. > solve({seq(h(t)[n]=s(u)[n],n=..3)},{t,u}); { t = arctan (, RootOf ( _Z 3, label = _L)), u = arctan ( RootOf ( _Z 3, label = _L)) }, { t = arctan ( RootOf (_Z 3, label = _L )) signum ( RootOf ( _Z π 3, label = _L)) π +, u = arctan ( RootOf ( _Z 3, label = _L) ) + signum ( RootOf ( _Z 3, label = _L)) π} > evalf(%); { t = , u = }, { u = , t = } > h( ); e e + y e > s( ); e e e y > h( ); ( ) e e + y e > s( ); ( ) e e e y So what has happened to the other two points of intersection? Note that the points of intersection that Maple found are those that have a positive y component. From our knowledge of symmetry we know that the other two points of intersection will be when the y component is negative. Fortunately, we can use Maple s solve command on one of the equations to find the other two solutions. If we solve the component equation we get, > solve(sin(t)^-cos(t)^=/); π,,, 3 π π 3 3 π 3 Note that the two numeric values for the parameter of the s curve that Maple gave us π π were the 3 and 3. If we evaluate s at the other two values we get > s(-/3*pi); e 3 + e y e

22 > s(-/3*pi); - + e 3 e y e which have the same and coordinates as the first two points but the y coordinate is negative. The only remaining question is what are the corresponding values for the parameter for the h curve? If we set the components equal to each other while π π substituting the 3 and 3 in for the s curve we obtain the two values that Maple found. > solve(sin(u)=cos(-pi/3)); π 6 > solve(sin(u)=cos(-*pi/3)); π 6 π π If we set the y components equal to each other while substituting the 3 and 3 in for 5π the s curve we get the value 3, which checks. > solve(cos(u)=sin(-pi/3)); > solve(cos(u)=sin(-*pi/3)); 5 π 6 5 π 6 > h(5/6*pi); e 3 + e y e From the graph and our knowledge of symmetry we know that the remaining parameter 7π value is 3 > h(7/6*pi); - + e 3 e y e Although Maple was a powerful tool in aiding the solution to this last problem it did not just hand us the solution. Recall the adage that says not only do you need the right tool for the right job but you also need to know how to use it.

23 We will do one final eample of solving for space curve intersections. Consider the two curves r ( t ) = sin( t),cos( t), t and s ( t) = cos( t),sin( t),sin ( t) cos ( t) defined. If we graph them we get the following image. > display(plotvvf(r,0,*pi,500,black,boed), plotvvf(s,0,*pi,500,black,boed));, both previously If we rotate the image around we see that there appears to be only one point of intersection between the two curves. You may want to increase the range of the parameter bounds in order to convince yourself that there will be only one point of intersection. Using the solve and evalf commands we get the following result. > solve({seq(r(t)[n]=s(u)[n],n=..3)},{t,u}); 4 { t = RootOf (_Z 4 sin ( arccos ( sin (_Z)) ) + 4 sin ( arccos ( sin (_Z) )) ) u = arccos ( RootOf ( _Z 4 sin ( arccos ( sin (_Z) )) + 4 sin ( arccos ( sin ( _Z))) ) )} > evalf(%); { u = I, t = I } 4, As you can see, the numeric approimations to the parameter values have a substantial I component, one that cannot be ignored because of round-off error. Now if we did not see the image we would be inclined to think that the two curves did not intersect. Since we know that they do we need to use another method to find the point of intersection. Recall that the fsolve command finds real solutions to equations and systems of equations using numeric methods but simply changing solve to fsolve gives us an error since the fsolve command needs to have the same number of variables as it does equations. When solving for the parameters of two curves in three dimensions we have three equations and only two variables.

24 > fsolve({seq(r(t)[n]=s(u)[n],n=..3)},{t,u}); Error, (in fsolve) number of equations, 3, does not match number of variables, One trick to get around this is to add a third, bogus, variable to the list of variables. Be careful that the variable you add has not already been defined as some value or epression. We will add a to the list. > fsolve({seq(r(t)[n]=s(u)[n],n=..3)},{t,u,}); { = , t = , u = } Checking the t and u values, ignoring, gives us the same point with a little round-off error. > r( ); e e e y > s( ); e e e y Eercises: t t. The following eercises concern the vector-valued function ( t) = e, tan ( t), t 3 b ( t) = e, tan t, t. t t a. Define the vector-valued function ( ) 3 b. b. Find the tangent vector function of b and define it as tb. c. Find the limit of the norm of the tangent vector function as t approaches infinity. d. If t represents time what does the answer in part c represent? e. If you are a bug that is walking along the curve, in step with t, describe your journey from t = 5 to t =. f. Plot the curve from t = 5 to t = 0 and on another graph plot the norm of the tangent vector function also from t = 5 to t = 0. Again, thinking of t as time relate the norm graph to the curve s graph in terms of position, direction and speed.. The following eercises deal with plots of the norm of the tangent vector function. a. Create a command called normplot that takes thee inputs, the curve, a lower bound for t and an upper bound for t and constructs the graph of the norm of the tangent vector function. That is, to plot the norm of the tangent vector function of a vector-valued function r from t = 5 to t = 0 the command would be, > normplot(r,5,0);

25 b. Create a command called normplotresy that takes 5 inputs, the curve, a lower bound for t and an upper bound for t, the lower bound for y and the upper bound for y and constructs the graph of the norm of the tangent vector function with the range restricted to the lower and upper bounds for y. That is, to plot the norm of the tangent vector function of a vector-valued function r from t = 5 to t = 0 and a restricted range from to 7 the command would be, > normplotresy(r,5,0,,7); c. Use these new commands to graph the tangent vector norm plots of the following functions. In each case give a detailed description the motion of a bug along each of the curves. i. r ( t ) = sin( t),cos( t), t ii. s ( t) = cos( t),sin( t),sin ( t) cos ( t) iii. k ( t) = 3 t, t t, t 4 t iv. h ( t) = t,tan( t), ln( t) v. q ( t ) = t cos( t), t sin( t), t vi. p 3 ( t ) = t, t, t 3. Find the angle of intersection (smallest angle) between the curves k ( t) = t, t t, t t and j ( t) = t, t t, t t at each of the intersections points. Plot an image of the curves and tangent vectors used to determine the intersection angle for each point of intersection. 4. Find the angle of intersection (smallest angle) between the curves r ( t ) = sin( t),cos( t), t and ( t) = cos( t),sin( t),sin ( t) cos ( t) s at each of the intersections points. Plot an image of the curves and tangent vectors used to determine the intersection angle for each point of intersection. 5. Finding the points of intersection. a. Define and plot the following two curves, plot them in different colors. c ( t) = sin( t),cos( 5t),cos( 3t) and d ( t) = cos( t),sin( 3t ),sin( 5t) b. From the graph, how many point of intersection are there? c. Define a command curveint that takes two vector-valued functions as input and does the solve command on them. So to find the points of intersection of the two curves above you simply need to type, > curveint(c,d); d. Use this new command to help find all the points of intersection between c and d. 6. Find all of the points of intersection between the two curves,

26 ( 4 + sin( 0t) ) cos( t),( 4 sin( 0t) ) sin( t),cos( t) f ( t) = + 0 and g ( t) = 4sin t,4cos t,sin 5t ( ) ( ) ( ) 7. Consider the two curves ( t) = sin( t),cos( t), sin( nt) r ( t) cos( t),sin( t), cos( nt) = r and, where n is a positive integer. Is there a relationship between the number of points of intersection of the curves and the value of n? If so, what is it? r ( t) sin t,cos t, sin nt and 8. Consider the two curves = ( ) ( ) ( ) r ( t) cos( t),sin( t), cos( mt) =, where n and m are positive integers. Is there a relationship between the number of points of intersection of the curves and the values of n and m? If so, what is it? 9. Create a new command called plotvvftv that does the same thing as plotvvftv ecept that instead of inputting evf, nlb, and nub you simply input a single number nv. The command should place nv evenly spaced tangent vectors along the plot of the curve starting at the beginning point of the curve and ending at the last point on the curve. So the beginning of the plotvvftv command should be, > plotvvftv:=(r,lb,ub,pts,col,as,nv,shp,vcol)-> 0. Use the command in the last eercise to plot the following curves along with a set of 30 evenly spaced tangent vectors along the portion of the curve that is being plotted. In each case, describe the motion of a particle moving along the curve considering t as time. Note that you may also need the normplot command to help describe the motion. 3 4 a. k ( t) = t, t t, t t b. h ( t) = t,tan( t), ln( t) c. c ( t) = sin( t),cos( 5t),cos( 3t) d. d ( t) = cos( t),sin( 3t ),sin( 5t) e. f ( t) = ( 4 + sin( 0t) ) cos( t),( 4 + sin( 0t) ) sin( t),cos( 0t) f. g ( t) = 4sin( t),4cos( t),sin( 5t)