Finally, we have all the tools laid out and ready, we can focus on doing something constructive, or at least interesting, with the math we have learned.

Often in architecture, engineering, physics, math, chemistry, design, artwork, or any number of creative fields, we need to build paths (lines and curves) to describe motion or outline objects.

We already have drawn lines and hooked them together to make straight “angular paths”, like we did in orienteering, going in straight lines from point to point.

We start with a bunch of points placed in the order that we want our path to go through, P₀, P₁, P₂, P₃, P₄.  We listed 5 here, but there could be any number of points on our path.  We will also start by looking at just the “x” coordinate of the first segment between the first two points, P₀ and P₁.  As we have seen, we can expand this into more dimensions (x, y, z) later.

Let’s look at the, now familiar, equation for the x coordinate of a line between the first two points P₀:(x₀, y₀) and P₁:(x₁, y₁):

    x(S) = x₀ + D_x0*S

The independent variable “S” marks the distance traveled along this segment from P₀, and the D_x0 is the x-direction of this segment between P₀ and P₁.  “S” starts at S=0 (at P₀) and goes to S=S₁ (at P₁), where “S₁“ is the distance of the line segment between P₀ and P₁.

Since this first segment of the path stops when the line reaches the end point (x₁, y₁) after traveling the distance between them, S₁:

     x₁ = x(S₁) = x₀ + D_x0*S₁

thus

     x₁ – x₀ = D_x0*S₁        and
     D_x0 = (x₁ – x₀)/S₁

Substituting this expression for D_x0 into the equation of the line we get:

     x(S) = x₀ + (x₁ – x₀)*S/S₁

Notice that the ratio S/S₁ can be thought of as a percentage and ranges from 0 to 1 as S goes from 0 to S₁.  To avoid calculating lengths of the segments, we can use the substitution:

        t = S/S₁

into the equation for a line and we get:

       x(t) = x₀ + (x₁ – x₀)*t     where t goes from 0 to 1.

This turns out to be a very convenient way to write the equation of a line between two points.  With this equation, x(0) = x₀ and x(1) = x₁, and “t” marks the “percentage” traveled along the line segment, from 0% (t=0) to 100% (t=1).

Expanding this into two dimensions we get the equations for all 4 segments between the 5 points:

Segment 1:      (x(t), y(t)) = (x₀, y₀) + (x₁-x₀, y₁-y₀)*t
Segment 2:      (x(t), y(t)) = (x₁, y₁) + (x₂-x₁, y₂-y₁)*t
Segment 3:      (x(t), y(t)) = (x₂, y₂) + (x₃-x₂, y₃-y₂)*t
Segment 4:      (x(t), y(t)) = (x₃, y₃) + (x₄-x₃, y₄-y₃)*t

In each case, “t” goes from 0 to 1.  This is written out in full detail to make sure that you can see the pattern.  In two dimensions you can find the length of Segment 1 by S₁² = (x₁-x₀)² + (y₁-y₀)².  The others are not written out, the pattern should be obvious.

Alright, we can draw lines between the points. Now let’s look at drawing curves between the points.  This becomes a bit harder to pin down because there literally are an infinite number of curves that could be drawn between these same 5 points.  We will have to specify the kind of curves we want.

We will learn how to draw curve segments between the points and put them together smoothly to form any continuous path we desire.  We could even use this math to find an equation that describes the path of a fly walking around on the ceiling.

One common way of doing this is to use a method called “cubic splines”.  With this method the coordinates of each curve are written as a “cubic” equation in t:

    x(t) = c₀ + c₁*t + c₂*t² + c₃*t³

where we can pick any numbers we want for c₀, c₁, c₂, and c₃.  The numbers we pick determine the shape and position of the curve. We can also write this expression using the “dot” product we have discussed before:

    x(t) = (c₀, c₁, c₂, c₃)^.(1, t, t², t³)

The “c”s are called the coefficients and the powers of “t” are called the basis functions.  We can now take the derivative of the x-coordinate as we learned in the last post.

      dx/dt = (c₀, c₁, c₂, c₃)^.(0, 1, 2*t, 3*t²)

as you recall, we can use these derivatives to find the direction of the tangent line at any point “t” along the curve.  As it turns out, we now have all the tools necessary to draw a “cubic spline” curve between the points.

We can pick “t” to be a “percentage of the curve” where “t” starts at t=0, the beginning of the curve, and goes to t=1 at the end of the curve (similar to how we drew lines).  And so:

    x(0) = (c₀, c₁, c₂, c₃)^.(1, 0, 0, 0)  =  x₀ (Substitute t=0 into the equations above)
    x(1) = (c₀, c₁, c₂, c₃)^.(1, 1, 1, 1)  =  x₁ (Substitute t=1 into the equations above)

written out this means:

    c₀ + 0  + 0 +  0  = x₀           and
    c₀ + c₁ + c₂ + c₃ = x₁

and we have found that the coefficient c₀ = x₀ (the x-coordinate of the starting point) and the sum of the coefficients must add to x₁ (the x-coordinate of the ending point).  This meets the conditions that the curve will start at the beginning point and end at the ending point.  We still need some more conditions to determine what the other c’s will be.

It is natural to also place conditions on the direction that the curve leaves the starting point and the direction it arrives at the ending point.  This is determined by the tangent direction (derivatives) at these points:

dx/dt(0)   = (c₀, c₁, c₂, c₃)^.(0, 1, 0, 0) = S₀*D_x0
dx/dt(1)   = (c₀, c₁, c₂, c₃)^.(0, 1, 2, 3) = S₁*D_x1

We can set S₀, S₁ and the directions D_x0 and D_x1 to be whatever we want and we can immediately see:

                c₁                = S₀*D_x0
    c₁ + 2*c₂ + 3*c₃ = S₁*D_x1

We have determined c₀ and c₁ from the starting point and the direction of the tangent at that point.  We can use substitutions to find the rest.  We have 4 equations and 4 unknowns and so we do this:

      3*c₀ + 3*c₁ + 3*c₂ + 3*c₃ = 3*x₁   (multiply the sum by 3)
                       c₁ + 2*c₂ + 3*c₃ = S₁*D_x1  (from above)

Subtracting the second equation from the first and making substitutions for c₀ and c₁:

      3*c₀ +   2*c₁ +   1*c₂ + 0*c₃ = 3*x₁ – S₁*D_x1  (making subs)
      3*x₀ + 2*S₁*D_x0 + c₂              = 3*x₁ – S₁*D_x1   (solve for c₂)

          c₂ = 3*x₁ – 3*x₀ – S₁*D_x1 – 2*S₀*D_x0      or
         c₂ = 3*(x₁ – x₀) – (S₁*D_x1+2* S₀*D_x0)

          c₃ = x₁ – (c₀ + c₁ + c₂)     (rearranging the sum)
          c₃ = x₁ – [x₀ + S₀*D_x0 + 3*(x₁ – x₀) – (S₁*D_x1 + 2*S₀*D_x0)]
          c₃ = -2*(x₁ – x₀) + S₁*D_x1 + S₀*D_x0

and now we have found all the “c” coefficients.

We can now draw the curves between the points if we know the starting point (x₀, y₀), the ending point (x₁, y₁), the starting direction (D_x0,D_y0) and the ending direction (D_x1, D_y1). The S₀ and S₁ variables are arbitrary and can be varied to taste.  We now plug all of these coefficients into our original curve:

      x(t) = (c₀, c₁, c₂, c₃)^.(1, t, t², t³)
      dx/dt = (c₀, c₁, c₂, c₃)^.(0, 1, 2*t, 3*t²)   = S₀*D_x0 (tangent direction)

the y(t) equation follows exactly the same pattern (just substitute “y” for “x” in all the equations).  Voila, we have found a “cubic” curve that passes through these points and has these tangent directions.

We can now do the same for all of the curves between each pair of points, drawing curves through all of the points.  There is a little bit of wiggle room in these curve equations.  The S₀, S₁, … and direction parameters for the tangents at the points can be set at will.  Since the direction vectors are the same at the points, the curve will smoothly pass through the points in the same direction.  These direction vectors show the “nose” direction (D_x,D_y) and “speed” (S) that a fly would be traveling if it were walking along the curve. It is fun to program these equations into a computer that draws them and watch the curves change as you change the location of the points and the directions and magnitudes of the tangents at the points.  This is how most modern, curve-fitting software works.

Of course, you can use another set of “basis functions” rather than just the powers of “t”, if you wish, to possibly get a better fit on your curves.  You might also want to impose other conditions, such as minimum curvature.  Just remember, if you use more basis functions, you will have more equations to solve for, you will need more conditions, and you might need some experience in “linear algebra” to solve for the coefficients.  Happy curving!

We covered the concept of differential calculus in just the last post.  It is time to slow down just a bit, catch our breath, and explore more in depth what this all means.

First, we started out with the concept of a function.  We used the example of a curve in space where the x-coordinate (and all of the others) depend on an index “t” that indicates where each point is along this curve in space.

     x(t) = x-coordinate of a point at the index t

If you specify a “t” value, the function x(t) will give the x-coordinate of a point along the curve that corresponds to “t”.

It might not have been perfectly clear, but we assumed that the curve is made out of points that are infinitely close to each other (same as a line).  There are no gaps or kinks in the curve, such that when you increase “t” by an infinitesimally small amount “dt” it will return the coordinates of a point that is infinitesimally close along the curve.  This is the concept of a “smooth” or “continuous” function.

An example of this is a car traveling down a straight road.  Suppose the function x(t) shows how far the car has gone down the road at a time “t”.  At each time, “t”, the function x(t) will tell you the exact position of the car. Now since a car can’t instantaneously disappear from one point and appear in a new point, it always takes time to travel continuously between any two points, which means that if you look at the position of the car a billionth of a second later, the position will be very close to where it was before.

    As we discussed, “dt” is much smaller than a billionth of a second, it is the smallest amount possible without being zero.  If “t” has increased by a small amount, “t+dt”, then the x-coordinate also will have changed by a very small amount “dx”. This is what we mean by the expression:

           dx = x(t+dt) – x(t)

The derivative dx/dt is defined to be:

     dx/dt  = [x(t+dt) – x(t)]/dt

which indicates how fast “x” is changing in relation to a change in “t” at the point “t”.  We can now hopefully see why differential calculus is the mathematics of motion.  In the example of the car above, dx/dt is velocity of the car at any instant “t” along its path – it moves a distance dx in the time dt.

    In the example where:

    x(t)       =      t²
 x(t+dt)   =  (t+dt)²

   dx = x(t+dt) – x(t)
         =  (t+dt)² –  t²
         = t²+2*t*dt + dt² – t²   (the t² cancels out)
         = 2*t*dt + dt²

and

    dx/dt = [x(t+dt) – x(t)]/dt
                 = 2*t       (when dt –> zero)

We say that the derivative of t² is 2*t, which also indicates how fast x is changing in ratio to the change in t at every point along the curve.

This function and differential notation is very powerful, we can get the answers to some questions without even knowing the exact function.

Suppose that x(t) is the sum of two functions of t, f(t), g(t):

    x(t) = f(t) + g(t)    where f(t) and g(t) are smooth

     df = f(t+dt) – f(t)
     dg = g(t+dt) – g(t)

dx =         x(t+dt)                –       x(t)
         = f(t+dt) + g(t+dt)     –    (f(t) + g(t))
         = [f(t+dt) – f(t)] + [g(t+dt) – g(t)]
         =             df           +          dg

Thus the derivative is also the sum of the derivatives:

    dx/dt =  df/dt  + dg/dt.   

Suppose that x(t) is the product of two functions of t:

    x(t) = f(t)*g(t)    where f(t) and g(t) are smooth

   dx =           x(t+dt)          –       x(t)
         = f(t+dt)*g(t+dt)     –    f(t)*g(t)
         =[df + f(t)]*[dg + g(t)] – f(t)*g(t)
         = df*dg + df*g(t) + f(t)*dg + f(t)*g(t) – f(t)*g(t)
         = df*dg + df*g(t) + f(t)*dt

Thus the derivative is:

dx/dt = df*(dg/dt) + (df/dt)*g(t) + f(t)*(dg/dt)
            = (df/dt)*g(t) + f(t)*(dg/dt)     because df –> 0

dx/dt = (df/dt)*g(t) + f(t)*(dg/dt)

Notice that we must always divide differentials by differentials to get derivatives, and all differentials standing alone (dx, df, dg, dt) go to zero.

We can now use the product rule to get the derivative of:

    x(t) = t³ = t² * t   {f(t) = t²,  g(t) = t}

we already know that df/dt = 2*t from above and

    dg/dt = [t+dt – t]/dt = 1

therefore since

dx/dt = (df/dt)*g(t) + f(t)*(dg/dt)

dx/dt =    (2*t)*t      +    t²*1
                  = 3*t²

By extending this we can show that for any power of t:

If x(t) = tⁿ  =  t^(n-1) * t

then

  dx/dt = n * t^(n-1)

We can see where some of the derivative tables in your calculus book come from.

We will end off by finding the derivative for a ratio of functions:

    x(t) = f(t)/g(t)

   dx/dt = [f(t+dt)/g(t+dt)        –    f(t)/g(t)]/dt
               = [(df+f(t))/(dg+g(t)) –    f(t)/g(t)]/dt
               = [[g(t)*(df+f(t))-f(t)*(dg+g(t))]/(g(t)*(dg+g(t))]/dt
  = [g(t)*(df/dt) + f(t)*g(t)/dt – f(t)*(dg/dt) – f(t)*g(t)/dt]/ (g(t)*dg + g(t)*g(t))

and when df and dg –> 0,

    dx/dt  = [g(t)*(df/dt) – f(t)*(dg/dt)]/(g(t)*g(t))

When we learn the “algebra” of differentials, we can find derivatives for a vast array of smooth functions.

As we have seen, when we know the derivatives of a curve, we can find the direction of the tangent line at any point along the curve. Recall from the last post in 2-D:

    (D_x, D_y) = (dx/dt, dy/dt)/((dx/dt)²+(dy/dt)²)^1/2

The nice thing about this notation is that it is easy to go to 3-D, 4-D, and as many -D’s that you want just by continuing the pattern.

Calculus deals with the mathematics of movement. In order to move from a starting point (x(S), y(S)), in 2-D, to the point next to it, you need to know how to move a tiny distance “dS” in some direction (D_x, D_y).  If we are moving along a line, then we already have a way to express this.  If (x(S+dS), y(S+dS)) is the function notation for the adjacent point along the line

(x(S), y(S)) = (x₀, y₀) + (D_x, D_y)*S

(we are writing this in function notation) then the line segment that connects them can be written (it works out):

       (x(S+dS), y(S+dS)) = (x(S), y(S)) + (D_x, D_y)*dS

Looking at just the x-coordinate this means:

     x(S+dS) = x(S) + D_x*dS   or

     x(S+dS) – x(S) = D_x*dS

Now it is pretty easy to see that when dS goes to zero, the distance between the x-coordinate of the ending point and the starting point  “x(S+dS) – x(S)” also goes to zero.  We can use our differential notation to express this:

       dx = x(S+dS) – x(S)

where “dx” is defined as the difference between the x-coordinates.  We now can write the differential expression:

       dx = D_x*dS

Recall that the small “dx” means differential (the infinitesimal distance moved, and the large “D_x” means the direction it moves.  See the diagram above.

Now we are ready to make the leap.  Even though dx and dS get as close as possible to zero, their ratio in this example is fixed:

                    dx/dS = D_x

This is because the Direction of the line segment is fixed and does not depend on the distance between points, no matter how small.

The concept of calculus is used to travel a tiny distance “dS” to arrive at an adjoining point.  No matter how small dS is, you can always cut it in half (dS/2) to find a point that is closer.  Recall that this concept allows us to “zoom in” to a single point throughout all eternity and never get to the point.  In calculus “dS” is called a “differential” and it is “infinitesimally” small (as close as possible to zero, without being zero).  This is the concept that makes Calculus a bit strange but also makes it interesting.

We can use the same concepts to describe what happens along a smooth non-linear curve (x(t), y(t)).  If we then move along the curve by increasing “t” by “dt”, an infinitesimal amount to the point (x(t+dt), y(t+dt)) then these two points on the curve are infinitesimally close to one another.  We define the differentials:

     (dx, dy) = (x(t+dt)-x(t), y(t+dt)-y(t))  or

        dx = x(t+dt) – x(t)

        dy = y(t+dt) – y(t)

Note that the variable “t” no longer denotes a distance, it is just an “index” that determines the position of the point on the curve.

Let’s consider an example, with just the x-coordinate.  Suppose

     x(t) = 1+ t²   then

     x(t+dt) = 1+ (t+dt)² = 1+ t² +2*t*dt + dt²

      dx = x(t+dt) – x(t) = 2*t*dt + dt²

      dx/dt = 2*t + dt

as dt goes to zero, we can ignore the “dt” term.

        dx/dt = 2*t

The line between these two points can be written as:

       (x(t+dt), y(t+dt)) = (x(t), y(t)) + (D_x, D_y)*dS

      (x(t+dt)-x(t), y(t+dt)-y(t)) = (D_x, D_y)*dS

Or   (dx(t), dy(t)) = (D_x, D_y)*dS

now dS² = dx(t)² + dy(t)² or

          dS² = [(dx(t)/dt)² + (dy(t)/dt)²]*dt²

which shows that “dS” also scales with “dt”.

Once we know these relationships then we can solve for the direction (D_x(t), D_y(t)) of the line that connects two infinitesimally close points on the curve (x(t), y(t)) at the index “t”.  We use the function notation for the direction parameters to indicate that the direction of this “tangent” line depends on where it is along the curve.

Just as in the case of a line, even if dS, dt, dx, and dy go to zero, as the points get closer together, they could have ratios between them that do not go to zero.  Everything depends on the shape of the curve.  Certainly, we expect that the direction of the line will not necessarily change as the points converge.  Such a curve is called “differentiable”.

Let’s look at an example.  Suppose:

   (x(t), y(t)) = (t, t²)              then

   (x(t+dt), y(t+dt)) = (t+dt, (t+dt)²) = (t+dt, t²+2*t*dt + dt²)

   (dx(t), dy(t)) = (x(t+dt)-x(t), y(t+dt)-y(t)) = (dt, 2*t*dt+dt²)

    dS² = dx² + dy² = dt²+(2*t*dt)² = [1+(2*t)²]*dt²

since    (dx(t), dy(t)) = (D_x, D_y)*dS

              (dt, 2*t*dt+dt²) = (D_x, D_y)* [1+(2*t)²]^1/2*dt

So    (D_x, D_y) = (dt, 2*t*dt+ dt²)/ [1+(2*t)²]^1/2*dt

and  (D_x, D_y) = (1, 2*t + dt)/ [1+(2*t)²]^1/2

And so, the direction of the line is well defined when “dt” goes to zero, and we add function notation:

     (D_x(t), D_y(t)) = (1, 2*t)/ [1+(2*t)²]^1/2

This line direction is called the tangent to the curve at the point (x(t), y(t)).  At this point, this line will only touch the curve at this one singular point.

Also

   dS/dt = [1+(2*t)²]^1/2

and all is solved for every location “t” on the curve.  Perhaps there are many questions remaining.  One of them being, for example, how to calculate the distance along the curve “S” from one point to another.  These questions will be left for other posts.

So far, we have only described lines where the distance between any two points on the line is well defined by a “metric”.  Before we look at how to draw curves, let’s remind ourselves what a line looks like in math language:

     (x, y) = (x₀, y₀) + (D_x, D_y)*S

The starting point on the line is (x₀, y₀) and the direction the line takes (D_x, D_y) in a cartesian coordinate system, and “S” is the distance from the starting point. Note that “S” serves as an index or indicator for every point on the line. In math language, “S” is known as the “independent variable”. It determines the exact point on a line a distance “S” from the starting point. Of course, we assume we already know the starting point (x₀, y₀) and direction (D_x, D_y) of the line.

In “function notation” we put the independent variable (S) next to the “dependent variables” (x and y) in parenthesis to remind ourselves that if we pick a number for “S”, it will give us an exact point (x(S), y(S)) on the line:

     (x(S), y(S)) = (x₀, y₀) + (D_x, D_y)*S

This “function” notation is extremely helpful. We can use this notation to replace the variable “S” with any number or symbol we wish. For example, let us replace the variable “S” with “t” (like a “cut-paste” of “S” with “t”):

     (x(t), y(t)) = (x₀, y₀) + (D_x, D_y)*t

We can even replace “S” with a mathematical statement:

     (x(S+3), y(S+3)) = (x₀, y₀) + (D_x, D_y)*(S+3)

The rule is that whatever goes in the parenthesis in function notation it is substituted into the “function” at the appropriate place.

Now we can start drawing curves by considering functions that are “non-linear” (they don’t draw a line).  We have already looked at how to draw a circle.  You may recall that we use the same equation for a line, but instead of varying the distance “S” we varied the direction (D_x, D_y) = (Cos(t), Sin(t)), where “t” is the angle off the x-axis.  So, the equation for a circle of radius S=1, centered at (x₀, y₀) = (0, 0) would be written:

        (x(t), y(t)) = (Cos(t), Sin(t))

where “t” goes from 0 to 360 degrees.  We could also draw just any arc of the circle by determining the range that the angle “t” spans, for example, a quarter circle is drawn when t goes from 0 to 90 degrees.  You get the idea.

We should mention at this point that we chose “t” as the “index” that indicates where any point is on the curve.  It doesn’t tell us how much distance traveled along the curve, or even the time that it takes for a point to travel along the curve.  It is just taken as a general parameter.  In the case of the circle of radius 1, we know that the circumference is 2*Pi over 360 degrees and so the distance traveled along the circle from 0 to the angle “t” is S = 2*Pi * (t/360).

Now we can have fun drawing a bunch of curves, using this function notation.

An oval that has major axis 3 and minor axis 2 has this equation:

        (x(t), y(t)) = (3*Cos(t), 2*Sin(t))

A spiral that starts at (0, 0) and goes to (1, 0) is drawn by this equation:

        (x(t), y(t)) = ((t/360)*Cos(t), (t/360)*Sin(t))

We can also draw curves in 3-dimensions just as easily; here is a spiral around a cone:

      (x(t), y(t), z(t)) = ((t/360)*Cos(t), (t/360)*Sin(t), t/360)

We can use any non-linear equation to draw a curve, like:

(x(t), y(t)) = (-2*t+3*t², 2*t²-5*t³+4*t⁴)

We could keep going all day, curves are always visually interesting.

Later we will see how we can hook curves together to form any shape we want, but first we need another tool to analyze curves, differential calculus.  It really is not as difficult as you might think, so keep tuned.

As we enter dimensions beyond our physical experience, it might help our understanding to look back and see patterns in smaller dimensions that still apply in larger dimensions.  Is there such a pattern when finding directions (D_x, D_y, D_z, D_t) in 4 dimensions or larger.  One obvious pattern, as dimensions increase, is the way we measure distance:

   1-D:         D_x² = 1
   2-D:         D_x² + D_y² = 1
   3-D:          D_x² + D_y² + D_z² = 1
   4-D:          D_x² + D_y² + D_z² + D_t² = 1

The pattern for higher dimensions is obvious. We expect this pattern to continue to any number of dimensions; we just keep adding on perpendicular triangles.  See the diagram above.

There is one problem with setting directions in higher dimensions that is not so obvious. For example, in 1-D, D_x² = 1 means that D_x = 1 or D_x = -1, this defines the directions to be along the positive or negative x-axis.  So far so good. In 2-D, however, when D_y² = 1, it forces D_x² = 0, which does not indicate either a positive or negative direction along the x-axis.  This ambiguity is due to the rotational symmetry of perpendicular lines, as we have mentioned before, you can “flip” negative and positive directions around a perpendicular direction, and nothing changes.

This ambiguity of direction exists whenever we add a new dimension.  In 3-D, for example, when D_z² = 1, it forces D_x² + D_y² = 0 and there are ambiguities as to where to place the direction of both the x-axis and y-axis.  To a man standing exactly on the north pole, for example, every direction is south, there is no east and west.  I have also heard it said that no matter how you comb the fuzz flat on a tennis ball, there will always be a part (singularity).

To deal with these ambiguities when any of our direction parameters is zero, it is helpful to fix the angles of the x-axis, y-axis, z-axis, etc. so that directions will not be ambiguous when you pass through the “singularities” where one or more of the direction parameters are forced to be zero.  So a man standing on the north pole, where every direction is south, still knows which angle to face if he wants to point to Paris, France.  This makes defining directions consistent.

Let us look at a way to define directions as we are increasing dimensions. We defined the angle A_xy to be the angle from the x-axis in the x-y plane.  Then we defined the angle A_z to be the angle off of the x-y plane in the z direction.  Listing the direction parameters for the dimensions we know:

   2D: D₂ = (D_2x, D_2y) = (Cos(A_xy), Sin(A_xy))
   3D: D₃ = (D_3x, D_3y, D_3z) = (D_2x*Cos(A_z), D_2y*Cos(A_z), Sin(A_z))
                 = (D_2x, D_2y, 0)*Cos(A_z) + (0, 0, 1)*Sin(A_z)

We added the “2” and “3” subscripts to indicate the dimension.

It might take you a few minutes to see the pattern.  Expanding into 4-D, we define the angle A_t to be the angle off the x-y-z space toward the t-axis, then continuing the pattern, the direction in 4-D is:

    D₄ = (D_4x, D_4y, D_4z, D_4t) = (D_3x, D_3y, D_3z, 0)*Cos(A_t) + (0, 0, 0, 1)*Sin(A_t)

    So putting them together and stacking them we get:

         D_4x = Cos(A_xy)*Cos(A_z)*Cos(A_t)
         D_4y = Sin(A_xy)*Cos(A_z)*Cos(A_t)
         D_4z =                   Sin(A_z)*Cos(A_t)
         D_4t =                                   Sin(A_t)

And we have found the Direction parameters in 4-D, and hopefully we can see the pattern that will allow us to use the same scheme to find directions in any dimension. An interesting side note: The angle in the t-direction, A_t, determines how close the 4-D direction is to the t-axis. When it gets close to 90 degrees, the direction is close to the t-axis and Cos(A_t) is close to zero, and D_4x, D_4y and D_4z are close to zero, this restricts motion in the x-y-z space. The faster time is moving, the more restricted the motion is in 3-D space. How strange. We will see just how strange in the following post.