switch * cot to / tan since cot not in PDL

Fourmilab · Sep 29, 2021 · ccb2f6c · ccb2f6c
1 parent 6400565
commit ccb2f6c
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/perl/fbench.pl b/src/perl/fbench.pl
@@ -184,8 +184,8 @@ sub transit_surface {
            my $sagitta = sin(($old_axis_slope_angle + $iang) / 2);
            $sagitta = 2 * $radius_of_curvature * $sagitta * $sagitta;
            $object_distance = (($radius_of_curvature * sin(
-              $old_axis_slope_angle + $iang)) *
-              cot($axis_slope_angle)) + $sagitta;
+              $old_axis_slope_angle + $iang)) /
+              tan($axis_slope_angle)) + $sagitta;
            return;
         }
         my $rang = -asin(($from_index / $to_index) *