I have the following formulae that I found:
h[S_, H_, D1_] := (S*H)/D1;
α[L_, D1_] := ArcTan[L/D1];
b[S_, B_, D1_, L_] := (S*B)/((D1^2 + L^2)^0.5);
c[S_, B_, D1_, L_] :=
S*Tan[ArcCos[((D1 - B/2)^2 + (D1 + B/2)^2 + 2 L^2 - B^2)/(
2 Sqrt[(D1 - B/2)^2 + L^2] Sqrt[(D1 + B/2)^2 +...