15 October 2022 6:42:25.779 AM ellipse_test(): FORTRAN90 version Test ellipse(). ellipse_area1_test(): ellipse_area1() computes the area of an ellipse. R = 10.0000 Matrix A in ellipse definition x*A*x=r^2 5.00000 1.00000 1.00000 2.00000 Area = 104.720 ellipse_area2_test(): ellipse_area2() computes the area of an ellipse. Ellipse: 5.00000 * x^2 + 2.00000 * xy + 2.00000 * y^2 = 10.0000 Area = 104.720 ellipse_area3_test(): ellipse_area3() computes the area of an ellipse. Ellipse: x^2/ 10.0000 ^2 + y^2 / 3.33333 ^2 = 1 Area = 104.720 ellipse_aspect_ratio_test(): ellipse_aspect_ratio() computes the aspect ratio of an ellipse. A B Ratio 1.0 0.0 0.000000 1.0 0.1 0.100000 1.0 0.2 0.200000 1.0 0.3 0.300000 1.0 0.4 0.400000 1.0 0.5 0.500000 1.0 0.6 0.600000 1.0 0.7 0.700000 1.0 0.8 0.800000 1.0 0.9 0.900000 1.0 1.0 1.000000 ellipse_eccentricity_test(): ellipse_eccentricity() computes the eccentricity of an ellipse. A B Ecc 1.0 0.0 1.000000 1.0 0.1 0.994987 1.0 0.2 0.979796 1.0 0.3 0.953939 1.0 0.4 0.916515 1.0 0.5 0.866025 1.0 0.6 0.800000 1.0 0.7 0.714143 1.0 0.8 0.600000 1.0 0.9 0.435890 1.0 1.0 0.000000 ellipse_flattening_test(): ellipse_flattening() computes the flattening of an ellipse. A B Flat 1.0 0.0 1.000000 1.0 0.1 0.900000 1.0 0.2 0.800000 1.0 0.3 0.700000 1.0 0.4 0.600000 1.0 0.5 0.500000 1.0 0.6 0.400000 1.0 0.7 0.300000 1.0 0.8 0.200000 1.0 0.9 0.100000 1.0 1.0 0.000000 ellipse_point_dist_2d_test(): ellipse_point_dist_2d() is given a point P, and finds the distance to an ellipse in 2D. The ellipse is (X/R1)^2 + (Y/R2)^2 = 1 R1 = 3.000000 R2 = 2.000000 P DIST -1.2000 3.9000 2.0124 -0.8000 3.6000 1.6524 -0.4000 3.3000 1.3138 0.0000 3.0000 1.0000 0.4000 2.7000 0.7154 0.8000 2.4000 0.4654 1.2000 2.1000 0.2570 1.6000 1.8000 0.1000 2.0000 1.5000 0.0080 2.4000 1.2000 0.0000 2.8000 0.9000 0.0978 3.2000 0.6000 0.3115 3.6000 0.3000 0.6231 4.0000 0.0000 1.0000 4.4000 -0.3000 1.4164 4.8000 -0.6000 1.8568 5.2000 -0.9000 2.3125 ellipse_point_near_2d_test(): ellipse_point_near_2d() is given a point P, and finds the nearest point PN on an ellipse in 2D. The ellipse is (X/R1)^2 + (Y/R2)^2 = 1 R1 = 3.000000 R2 = 2.000000 P PN -1.2000 3.9000 -0.8237 1.9231 -0.8000 3.6000 -0.5835 1.9618 -0.4000 3.3000 -0.3094 1.9893 0.0000 3.0000 0.0000 2.0000 0.4000 2.7000 0.3450 1.9867 0.8000 2.4000 0.7239 1.9409 1.2000 2.1000 1.1326 1.8520 1.6000 1.8000 1.5623 1.7074 2.0000 1.5000 1.9959 1.4931 2.4000 1.2000 2.4000 1.2000 2.8000 0.9000 2.7198 0.8440 3.2000 0.6000 2.9087 0.4897 3.6000 0.3000 2.9842 0.2049 4.0000 0.0000 3.0000 0.0000 4.4000 -0.3000 2.9920 -0.1457 4.8000 -0.6000 2.9761 -0.2522 5.2000 -0.9000 2.9582 -0.3327 ellipse_points_2d_test(): ellipse_points_2d() returns points on an ellipse; Ellipse center at 5.00000 , -2.00000 radii R1 = 3.00000 R2 = 1.00000 and angle PSI = 0.523599 and area = 9.42478 7.59808 -0.500000 7.20897 -0.282767 6.48356 -0.326967 5.53230 -0.625872 4.50000 -1.13397 3.54382 -1.77392 2.80933 -2.44829 2.40835 -3.05441 2.40192 -3.50000 2.79103 -3.71723 3.51644 -3.67303 4.46770 -3.37413 5.50000 -2.86603 6.45618 -2.22608 7.19067 -1.55171 7.59165 -0.945594 ellipse_points_arc_2d_test(): ellipse_points_arc_2d() returns points on an elliptical arc. The ellipse has center 5.00000 -2.00000 radii R1 = 3.00000 R2 = 1.00000 and angle PSI = 0.523599 The arc extends from THETA1 = 1.57080 to THETA2 = 6.28319 Sample points: 4.50000 -1.13397 3.54382 -1.77392 2.80933 -2.44829 2.40835 -3.05441 2.40192 -3.50000 2.79103 -3.71723 3.51644 -3.67303 4.46770 -3.37413 5.50000 -2.86603 6.45618 -2.22608 7.19067 -1.55171 7.59165 -0.945594 7.59808 -0.500000 ellipse_test(): Normal end of execution. 15 October 2022 6:42:25.779 AM