VERB_code_2.3
Y.m
1 function res = Y(alpha)
2  T0 = 1.3802;
3  T1 = 0.7405;
4  res = 2.*(1.0 - sin(alpha)).*T0 + (T0 - T1).*(sin(alpha).*log(sin(alpha)) + 2.*sin(alpha) - 2.*sqrt(sin(alpha)));
5 return
double Y(double alpha)
Y(α) function.