dsqr.c

/* dsqr(x) - square of a real number
             Argument must be type double
*/
double dsqr(x)
double x;
{
   return(x*x);
}