Solver SSA Triangle Solver2 - 2 sides and adjacent angle

Source code of 'SSA Triangle Solver2 - 2 sides and adjacent angle'

This Solver (SSA Triangle Solver2 - 2 sides and adjacent angle) was created by by tomlaidlaw(14)

: View Source, Show, Put on YOUR site
About tomlaidlaw: 74 year old man interested in mathematical recreations.

==section input
<img src=http://tomlaidlaw.com/goatandsilo/anytriangle.gif>
Input the following elements of a triangle:
 side a = *[input sidea=3], side c = *[input sidec=4], and Angle C = *[input anglec=50]
Find the remaining side and the other angles, Area, etc.
There may be 0, 1, or 2 triangles depending on your entries.

==section solution
perl
my $pi=3.1415926536;
my $a=$sidea;
my $ar=sprintf("%.3f", $a);
my $c=$sidec;
my $cr=sprintf("%.3f", $c);
my $C1=$anglec;
my $Cr=sprintf("%.3f", $C1);
my $C2=$C1*$pi/180;
my $C2r=sprintf("%.3f", $C2);
my $h=$a*sin($C2);
my $hr=sprintf("%.3f", $h);
my $sinC=sin($C2);
my $cminush = abs($c-$h);
print "$cminush";
if($cminush < .001)
{print "right triangle.";
 goto One}
if ($c<$h)
{
print "Given , side a = $a, side c = $c, Angle C = $C1.
h = $hr  NO TRIANGLE since side c is less than height, try again. ";
goto End}

if($c>=$a)
{print "Given  side a = $a, side c = $c, Angle C = $C1.
 ONE TRIANGLE since side c is greater than or equal to height.";
goto One}
if($c>$h && $c<$a)
{print "Given  side a = $a, side c = $c, Angle C = $C1.
 TWO TRIANGLES since side c is greater than height but less than side a. ";
}
One:
print
 "

Use the Law of Sines {{{sin(A)/a = sin(B)/b = sin(C)/c}}}";
my $sinA2=$a*$sinC/$c;
my $A2=atan2($sinA2,(sqrt(1-$sinA2**2)));
my $A1=$A2*180/$pi;

my $B1=180-$C1-$A1;
my $B2=$B1*$pi/180;
my $sinB=sin($B2);

my $b=$sinB*$c/$sinC;

my $Ar=sprintf("%.3f",$A1);
my $A2r=sprintf("%.3f",$A2);
my $Br=sprintf("%.3f",$B1);
my $B2r=sprintf("%.3f",$B2);
my $Cr=sprintf("%.3f",$C1);
my $C2r=sprintf("%.3f",$C2);

my $ar=sprintf("%.3f",$a);
my $br=sprintf("%.3f",$b);
my $cr=sprintf("%.3f",$c);

my $s = 1/2*($a+$b+$c);
  my $K=sqrt($s*($s-$a)*($s-$b)*($s-$c));
my $Kr=sprintf("%.3f", $K);

my $R=.5*$a/sin($A2);
my $Rr=sprintf("%.3f", $R);

my $r=2*$K/($s*2);
my $rr=sprintf("%.3f", $r);

my $ha=$c*sin($B2);
my $har=sprintf("%.3f", $ha);

my $hb=$a*sin($C2);
my $hbr=sprintf("%.3f", $hb);

my $hc=$b*sin($A2);
my $hcr=sprintf("%.3f", $hc);
print
"
side a = $ar
side b = $br
side c = $cr

AngleA = $Ar degrees or $A2r radians
AngleB = $Br degrees or $B2r radians
AngleC = $Cr degrees or $C2r radians

Area K = $Kr

Circumcircle Radius = $Rr __(R = .5a/sinA)
incircle radius = $rr __(r = 2K/(a+b+c))

height over side a = $har __(ha = csinB)
height over side b = $hbr __(hb = asinC)
height over side c = $hcr __(hc = bsinA)
h=$hr
";

if($c>=$a)
{goto End}
print
 "

second triangle";
my $A2s=$pi-$A2;
my $A1s=$A2s*180/$pi;
my $sinA2s=sin($A2s);
my $A2sr=sprintf("%.3f",$A2s);
my $A1sr=sprintf("%.3f",$A1s);

my $B2s=$pi-$C2-$A2s;
my $B1s=$B2s*180/$pi;
my $sinBs=sin($B2s);
my $B2sr=sprintf("%.3f",$B2s);
my $B1sr=sprintf("%.3f",$B1s);

my $bs=$sinBs*$c/$sinC;
my $bsr=sprintf("%.3f",$bs);
 my $ss = 1/2*($a+$bs+$c);
  my $Ks=sqrt($ss*($ss-$a)*($ss-$bs)*($ss-$c));
my $Ksr=sprintf("%.3f", $Ks);

my $hs=$a*$sinC;
my $hsr=sprintf("%.3f", $hs);
my $Rs=.5*$a/sin($A2s);
my $Rsr=sprintf("%.3f", $Rs);

my $rs=2*$Ks/($ss*2);
my $rsr=sprintf("%.3f", $rs);

my $has=$c*sin($B2s);
my $hasr=sprintf("%.3f", $has);

my $hbs=$a*sin($C2);
my $hbsr=sprintf("%.3f", $hbs);

my $hcs=$a*$sinBs;
my $sinBs5 = 5*$sinBs;
my $hcsr=sprintf("%.3f", $hcs);

print
"
side a = $ar
side b = $bsr
side c = $cr

AngleA = $A1sr degrees or $A2sr radians
AngleB = $B1sr degrees or $B2sr radians
AngleC = $Cr degrees or $C2r radians
Area K2 = $Ksr

Circumcircle Radius = $Rr __(R = .5a/sinA)
incircle radius = $rsr __(r = 2K/(a+b+c))

height over side a = $hasr __(ha = csinBs)
height over side b = $hbsr __(hb = asinC)
height over side c = $hcsr __(hc = bsinAs)

End:
print
"
Tom Laidlaw";
==section output
a,

==section check
  sidea=3, sidec=4 anglec=50