function lngamma1 input Real x1; input Real x2; input Real x3; input Real T; output Real lngamma; end lngamma1;