As with most any problem, this too may be approached from a multitude of directions.
The tack I have chosen distills to the determination of the discrete among a locus of points at the coping intersection of the two cylinders at right angles.
For reference and context please see ALL prior posts associated with A Precursory Study of Pythagorean Theorem Generalization a Consequence of Three-Dimensions
Example:
Catheti
a = 3
b = 4
Hypotenuse
c = √(a2 + b2) = 5
To a requisite number of decimal places, exhaustive iteration of D1z against the defined a,b will produce values for D1x,D1y and allow for the application of Heron’s formula to lateral face ABD1.
Apex Singularity X,Y Coordinates
D1z = 5.58670988566614
D1y = -√ ((2a)2 – D1z2) = -2.18830360174273
D1x = -√ ((2b)2 – D1z2) = -5.72613941966140
Lateral Face ABD1 Edge Lengths
BD1 = √ ((a + |D1x|)2 + D1y2 + D1z2) = 10.5898776749955
AD1 = √ (D1x2 + (b + |D1y|)2 + D1z2) = 10.1141040862422
ABD1 Semi-Perimeter
s = (AD1 + BD1 + c) / 2 = 12.8519908806189
ABD1 Area
√ (s (s – AD1) (s – BD1) (s – c)) = a2 + b2 = 25.00000000000000
ABCD1 Volume
VP = (.5abD1z) / 3 = 11.1734197713323
Where this area converges with that of c2, the coordinates of the apex singularity for a Pythagorean compliant tetrahedron are known. (Additionally, pyramidal volume is derived.)
While this approach is conceptually simple, the large number of terms resulting from the expanded equation have prevented me from solving for D1z should that be possible.
Alternatively, the apex singularity can be viewed as residing on a hyperbolic paraboloid surface for a given D1x. D1zon the cylindrical surface is simply a projection of D1zhp from the hyperbolic paraboloid surface. An obvious deficiency in this approach remains the lack of relation between a and D1x.
Apex Singularity Coordinates on Hyperbolic Paraboloid Surface
D1x = -5.72613941966140
D1y = ±√ (-(4b2 – 4a2) + D1x2) = -2.18830360174273
D1zhp = D1x2 – D1y2 = 4b2 – 4a2 = 28.00000000000000
Apex Singularity Coordinate on Cylindrical surface
D1z = √ (4b2 – D1x2) = 5.58670988566614
BCD1 Lateral Face Height/Altitude
hBCD1 = 2a = 6.00000000000000
Lateral Face ABD1 Edge Lengths
CD1 = √ (D1 D1x2 + hBCD12) = 8.29389369677477
BD1 = ±√ (a2 + CD12 ±2√ (a2 CD12 – 4a2)) = 10.5898776749955
AD1 = ±√ (b2+CD12 ±2√ (b2 CD12 – 4b2)) = 10.1141040862422
Whether the position of this vertex or its usefulness in determining volume, the consequence of knowing the apex singularity thus allows for further examination of Harmonious Pythagorean Tetrahedra.