OK, I’m simply reposting what I said in the other thread, although I was also thwarted by the fact that the difference was not really within the capabilities of my trusty HP-48GX calculator. One of these days I’ll make it over to one of my school’s computer labs to try Maple or Matlab (actually probably both, Maple for the symbolic math, and Matlab for the number crunching. I believe that Matlab uses a very lardge number of significant figures, which would be nice)
[/quote]
Or, we could always do this the hard way. With the tidal variation of 6.7 meters, the earth’s equator becomes an ellipse 6.7 meters larger in it’s major axis than its minor axis [I’m actually referring to 1/2 of each axis here, which is consistent with how the numbers are used later]. The equation in cartesian coordinates for an ellipse is:
x[sup]2[/sup]/a[sup]2[/sup] + y[sup]2[/sup]/b[sup]2[/sup] = 1
We will use 6378300 m for b and 6378306.7 m for a. This gives us an equation of:
x[sup]2[/sup]/(4.0682796359E13 m[sup]2[/sup]) + y[sup]2[/sup]/(4.068271089E13 m[sup]2[/sup]) = 1
Now we need to choose a point about mid-tide. We will take (for simplictity sake) the point where the ellipse crosses the line x = y. This point is at (4510146.63079 m, 4510146.63079 m).
Now we need to figure out where the person’s eyes are. We will assume that they are 2 meters off the ground. The slope of the tangent line at that point is given by:
m = - b[sup]2[/sup]x[sub]0[/sub]/a[sup]2[/sup]y[sub]0[/sub]
Since x[sub]0[/sub] = y[sub]0[/sub], they cancel, leaving us with, and we get a result of -.999997899137.
Since the atual location of the person’s eyes is two meters perpendicular to this, we need to find the slope for that line. It is the negative reciprocal of the previous slope, or 1.00000210087. A bit of calculation shows the person’s eyes to be at the point (4510148.045 m, 4510148.04501 m). At about this point, I start wishing my calculator kept track of more significant figures.
Now we are simply left with finding the two tangents that intersect this point. If this is a trivial problem, I’m unaware of it, and hence have the long way here:
we know that these lines will be tangent to the ellipse, so we can easily get the slopes:
m[sub]1[/sub] = b[sup]2[/sup]x[sub]1[/sub] / a[sup]2[/sup]y[sub]1[/sub]
m[sub]2[/sub] = b[sup]2[/sup]x[sub]2[/sub] / a[sup]2[/sup]y[sub]2[/sub]
We also know that they pass through the point (x[sub]0[/sub], y[sub]0[/sub]), and that the points (x[sub]1[/sub], y[sub]1[/sub]) and (x[sub]2[/sub], y[sub]2[/sub]) will be on the ellipse, and thus fulfill the basic ellipse equation of x[sup]2[/sup]/a[sup]2[/sup] + y[sup]2[/sup]/b[sup]2[/sup] = 1.
Thus we have two simultaneous equations for each point:
(y[sub]0[/sub]-y[sub]1[/sub]) / (x[sub]0[/sub]-x[sub]1[/sub]) = b[sup]2[/sup]x[sub]1[/sub] / a[sup]2[/sup]y[sub]1[/sub]
and
x[sub]1[/sub][sup]2[/sup]/a[sup]2[/sup] + y[sub]1[/sub][sup]2[/sup]/b[sup]2[/sup] = 1
and the same for (x[sub]2[/sub], y[sub]2[/sub])
After a bit of work, I reduced this to:
y[sub]1[/sub] = sqrt( (b[sup]2[/sup]/a[sup]2[/sup]) ( (a[sup]2[/sup]y[sub]0[/sub] / b[sup]2[/sup]x[sub]0[/sub]) + (a[sup]2[/sup]/x[sub]0[/sub]) )[sup]2[/sup] - b[sup]2[/sup])
Which is all well and good, except that there should be another possible answer for y[sub]2[/sub]. Perhaps someone who actually is good at this shit could find the other equation.
As for this one, we can now substitute:
a[sup]2[/sup] = 4.0682796359E13 m[sup]2[/sup]
b[sup]2[/sup] = 4.068271089E13 m[sup]2[/sup]
y[sub]0[/sub] = 4510148.04501 m
x[sub]0[/sub] = 4560148.045 m
At this point, I realize that my calculator seems unable to deal with the precision adequately. Maybe I’ll get around to figuring out how to solve it in Maple or Matlab or something. In the meantime, I invite others to give it a try.