void square(double *x, int *nin, double *result)
{
    int n = nin[0];

    int i;
    for (i = 0; i < n; i++)
        result[i] = x[i] * x[i];
}