Paddling workouts

I bought a kayak a couple of weeks ago. It’s one of those fold-up jobs by Oru, and I got it instead of a more traditional kayak because it’s very light and can easily fit inside my car instead of being strapped to the roof.

Folded kayak in SUV trunk

I’ve been taking it out to local lakes and paddling around to get myself used to it. The first couple of times, because I was focused on assembling it correctly and getting in without tipping over, I forgot to start the Paddling workout on my Apple Watch, so I didn’t have a complete record of either of those little trips. But I have recorded the more recent paddles and have learned a rather disappointing fact: the Workout app doesn’t keep track of where you went.

This was quite a shock. For years, I’ve been using Workout to track my walks, getting both a nice map of my path and a sense of how fast I was going during each stretch.

Hiking workout map

(The map can also show that you got in the car and drove away before turning the workout off.)

I assumed the Paddling workout would do the same, but here’s what I saw after paddling around at Busse Woods a few days ago:

Paddling workout map

It gives you a dot where you started the workout but nothing else. A search of the web to see if I’d done something wrong soon told me that I hadn’t and that lots of paddlers were unhappy about this deficiency.

Apple has always been concerned with battery life on the Watch, but I don’t understand why keeping track of my location while kayaking is going to be any more draining than doing the same thing while hiking. Or biking, for that matter, a workout type that also tracks location and presents a map.

There are third-party paddling apps that will track my location, so I guess I’ll have to start trying them out. My initial sense is that these apps are meant for people who are a lot more serious about kayaking than I ever expect to be and therefore have features I’ll never use. Oh, well. Every hobby has hidden costs.

Update 18 May 2024 5:42 PM
A few updates:

  1. On Mastodon, Steve told me that cross-country skiing also records just a single dot on the workout map. Just as stupid as how it treats kayaking.
  2. Joe Rosensteel wrote a nice followup post in which he pointed out another Apple Watch/Workout deficiency that you’d think could be easily fixed: while your watch will (a) notify you when it detects that you’ve started a walk without starting a workout for it, and (b) also notify you if it notices that you’re standing (mostly) still without pausing your current workout, it won’t do the obvious (c) notify you to resume your workout when you start walking again. I’ve been bitten by this bug several times.
  3. I’ve signed up for a free month of Strava rather than a paddling-specific app/service. I’m hoping this will make it easy to track both my walking and kayaking in a single place.

A factoring trick

In Matt Parker’s latest video, he gives two proofs to show that if you

you’ll get a number that is a multiple of your original number minus one. In words, this seems magical; in equation form,

b n1=(b1)(something)

not so much. Certainly you know that

b 21=(b1)(b+1)

because that comes up so often. You may also remember from algebra class that

b 31=(b1)(b 2+b+1)

So it isn’t tremendously surprising to hear that b1 is a factor no matter what n is. But what’s the proof?

Matt’s first proof is to notice that b=1 is a solution to

b n1=0

for any n. Therefore, (b1) has to be one of the factors of b n1. He considers this a rather prosaic proof, and quickly moves on to the fun one, which involves expressing the numbers in different bases. It is a cute proof, and you should watch the video to see it.

But I wanted to go back to

b n1=(b1)(something)

and see if there’s a solution for something that works for every n. Obviously there is, or I wouldn’t have written this post. And you’ve probably already guessed what that solution is, but let’s figure it out systematically.

We’ll set up the polynomial division of b n1 by b1 and run it out for a few steps:

Polynomial division

As you can see, after Step m in the process, the difference (i.e., the expression under the horizontal line) is

b nm1

So after m=n1 steps, the difference will be b1, which means that the last term in the quotient will be +1 and there will be no remainder. Therefore, something is

b n1+b n2++b+1=∑ m=0 n1b m

which is easier to remember than I would have guessed.


Gravity and Edgar Rice Burroughs

Yesterday, I was thinking about what textbook I would move on to after I finish Adventures in Celestial Mechanics. Because I’ve been thinking about gravity, I pulled Oliver Kellogg’s Foundations of Potential Theory off my shelf and started leafing through it. I soon came upon a gravity problem that I first heard about when I was a teenager but never saw the solution for. The problem has to do with Pellucidar, and it’s not as tricky as I expected.

Pellucidar is an imaginary land created by Edgar Rice Burroughs for a series of pulpy adventure books, starting with At the Earth’s Core. The conceit of the stories is that the Earth is actually a hollow shell and that there’s an inhabited land, Pellucidar, on the inside surface of the shell. Intrepid adventurers from our side of the shell dig through to Pellucidar, meet the strange people and creatures who live there, and have adventures intrepidly. Never one to let an idea go unexploited, Burroughs even had Tarzan visit Pellucidar.

While I read one or two Tarzan books and several of the John Carter of Mars books, I never picked up any of the Pellucidar novels. But I did read somewhere—possibly in an Isaac Asimov science essay—that the concept of Pellucidar was flawed in a way apart from the many obvious ways: the inhabitants of Pellucidar would not “stick” to the inner surface of a hollow Earth because there’s no net gravitational force inside a spherical shell.

This seemed very odd to the teen-aged me, but Asimov (or whomever) didn’t explain it, and I didn’t have the analytical tools to check it out myself. But now I do, following more or less the derivation given by Kellogg.

Suppose we have a thin spherical shell of radius a. We want to know the gravitational attraction that the sphere exerts on a particle, P, which may be either inside or outside the sphere.

Sphere with coordinate systems

We’ll set up our coordinate system so that P is on the z axis. Every point on the sphere, Q, can be identified by two angular coordinates, θ and ϕ, which are the point’s latitude and longitude.

(If you’re wondering how we should handle the analysis if P isn’t on the z axis, don’t. One of the most important things to learn in mechanics is that the coordinate systems we use are up to us; there is no x-y-z stamped on the universe. Wherever P happens to be, we can draw a z axis that runs through it from the center of the sphere.)

The gravitational force that some small differential element at Q on the sphere exerts on particle P is

dF=GρdAr 2

where G is the universal gravitational constant, ρ is the areal density of the sphere (i.e., it has units of mass per area, not mass per volume), and r is the distance from Q to P. dA is the differential area of the element on the sphere, and dF is the differential force (per unit mass of P) that element exerts on particle P.

Here’s a sketch of the differential area at Q:

Differential area element

It’s bounded by lines of latitude that are dθ apart from one another and lines of longitude that are dϕ apart from one another. The area of the element is

dA=(adθ)(acosθdϕ)=a 2cosθdϕdθ

If you’re wondering why we can use the formula for the area of a rectangle when the edges and the surface itself are curved, I will remind you that this is how calculus works. These differential lengths are teeny tiny and therefore the difference between the area of the surface and the area of a rectangle is extra extra teeny tiny. In calculus, we keep track of the teeny tiny and throw away the extra extra teeny tiny.

To express the differential force in terms of θ and ϕ, let’s look at a side view of the sphere with P and Q. First, we’ll consider the case where P is “above” Q:

Geometry with P above Q

The distance between P and Q is

r=(acosθ) 2+(zasinθ) 2=a 2+z 22azsinθ

Although I’ve drawn the figure with P outside the sphere, the formula above also applies when P is inside the sphere. What matters is that P is above Q.

Now let’s consider the case where P is “below” Q:

Geometry with P below Q

Here, r is the hypotenuse of a different right triangle, but it turns out that

r=(acosθ) 2+(asinθz) 2=a 2+z 22azsinθ

is the same formula as before. That’s convenient. So the magnitude of the differential force is

dF=Gρa 2cosθa 2+z 22azsinθdϕdθ

To determine the force on P due to the entire sphere, we’re going to integrate on θ from π/2 to π/2 and on ϕ from 0 to 2π. But first, we need to think about the direction of the differential force.

Note that in each of the two side views, the gravitational force on P has two components: a component parallel to the z axis (vertical in our drawings) and a component perpendicular to the z axis (horizontal in our drawings). When we integrate over ϕ, we’re dealing with forces associated with a circle of constant latitude, and the components perpendicular to the the z axis (horizontal) will cancel each other out. For example, the horizontal component when Q is at ϕ=180° cancels the horizontal component when Q is at ϕ=0°, the horizontal component for ϕ=270° cancels the horizontal component for ϕ=90°, and so on.

The upshot is that we only have to worry about the “vertical” component of the force, i.e., the component in the z direction. To get that component, we multiply dF by the z direction cosine.

When P is above Q, the z direction cosine is

zasinθa 2+z 22azsinθ

where the minus sign out in front indicates that the component is acting in the negative z direction. Remember, we’re computing the gravitational force on P due to Q, so the force acts from P toward Q.

When P is below Q, the z direction cosine is

asinθza 2+z 22azsinθ

Compare this formula with the previous one and you’ll see that they are algebraically identical. As with the expressions for r, the signs work out to give the same equation for both cases.

So the z component of the differential force acting on P due to Q is

dF z=Gρa 2cosθ(asinθz)(a 2+z 22azsinθ) 3/2dϕdθ

and the total force that the sphere exerts on P is

F=F z=∫ π/2 π/2∫ 0 2πGρa 2cosθ(asinθz)(a 2+z 22azsinθ) 3/2dϕdθ

This may seem awful, but it isn’t. First, notice that ϕ doesn’t appear in the integrand, so the integral over ϕ is just

∫ 0 2πdϕ=2π

We can pull that and the other constants out to get

F=2πGρa 2∫ π/2 π/2cosθ(asinθz)(a 2+z 22azsinθ) 3/2dθ

This still looks bad, but with tools like Mathematica and SymPy, integrations like this aren’t the pain in the ass they used to be. After a few seconds of chugging in Mathematica,

the answer pops out:

F=2πGρa 2z 2(sgn(az)1)

where the signum function (sgn) is +1 if the argument is positive and 1 if the argument is negative. Therefore, if P is outside the sphere (a<z),

F=2πGρa 2z 2(11)=4πGρa 2z 2

At this point, we recall that 4πa 2 is the area of the sphere, so that multiplied by the areal density, ρ , is the mass of the sphere, which we’ll call M. So for this case

F=GMz 2

which is just what we’d expect (remember that this is the force per unit mass of P, so there’s no mass term for the particle).

On the other hand, if P is inside the sphere (a>z),

F=2πGρa 2z 2(11)=0

which is the surprising (to me, anyway) result that the sphere exerts no net gravitational force on a particle inside it.

How can this be? Well, certainly it makes sense if P is at the center of the sphere. If we move P up from the center, there’s more mass below P but it is, on the average, farther away. It works out that the increasing mass and increasing distance balance each other to keep the net force at zero.

Oliver Kellogg didn’t have Mathematica. His suggestion was to do a change of variable from θ to r and then integrate. He doesn’t give any details of the process, just the two answers. But I went ahead and did it his way, too, using

r=a 2+z 22azsinθ dr=azcosθa 2+z 22azsinθdθ=azcosθrdθ

and

asinθz=a 2z 2r 22z

With those substitutions, the integrand becomes

12az 2a 2z 2r 2r 2

which is much easier to deal with than the original integrand. But we have to be careful with the limits. If z>a, the upper limit (θ=π/2) becomes

r=za

and the lower limit (θ=π/2) becomes

r=z+a

Therefore, the integral is

F=πGρaz 2&Integral; za z+aa 2z 2r 2r 2dr

where I got rid of the negative sign at the front of the integrand by exchanging the limits. With no trig functions to make life difficult, this works out very quickly to

F=4πGρa 2z 2

when P is outside the sphere. This is the same answer as before.

If z<a, the upper limit (θ=π/2) becomes

r=az

and the lower limit (θ=π/2) becomes

r=a+z

Therefore, the integral is

F=πGρaz 2&Integral; az a+za 2z 2r 2r 2dr

which works out to

F=0

when P is inside the sphere. Again, a repeat of the answer we got when integrating over θ.

While it was nice to see explicitly how Kellogg did the integration, it was faster to just let Mathematica do its thing. And I didn’t get any deeper physical insight doing the change of variable.

Anyway, now I understand the surprising thing I read about 45 or so years ago. And I’m downloading At the Earth’s Core.


More than 100% pregnant

Kieran Healy tweeted out a gem earlier this evening:

Yes, Professor Huberman, a probability of 120% is indeed “a different thing altogether.”1

While I find it fascinating that someone can run through an argument that clearly proves they are wrong and continue along as if they hadn’t, I don’t have any expertise in psychology, so I won’t try to explain it.

I do, however, have some expertise in probability, and Huberman is dancing around the edge of one of my favorite problems: If there’s a 1 in n chance of something happening in one try, what’s the probability of it happening in n tries?

Remarkably, the answer is not

n(1n)=1

It’s

1(11n) n

The term in the parentheses is the probability of it not happening in a single trial. That term raised to the nth power is the probability of it not happening in n trials (assuming independence between trials). And subtracting that from 1 gives us the probability of it happening (at least once) in n trials.

What’s great about this problem is that you can make a pretty decent estimate of the answer without doing any calculations. The reason is

lim n&rightarrow;(11n) n=1e&approx;0.37

and so

lim n&rightarrow;1(11n) n=11e&approx;0.63

Now, in the Huberman problem n is 5, which is not especially close to infinity, but luckily the convergence to the limit is pretty rapid.

Convergence of probability

I would have guessed the answer for n=5 was about 0.65; it’s actually 0.67. Not super accurate, but better than 1.

Update 10 May 2024 12:21 PM
John D. Cook also saw the Huberman clip and wrote a nice post this morning about how adding probabilities is a reasonably close approximation if both

  • the base probability for a single trial is small; and
  • the number of trials is small.

Absolutely true, but not as much fun as a solution that ends up with Euler’s constant popping up out of nowhere.


  1. The clip is from this 4 hour, 22 minute (!) video episode of the Huberman Lab Podcast. Imagine if Huberman were John Gruber’s guest on The Talk Show