```/* nag_nd_shep_interp (e01zmc) Example Program.
*
* Copyright 2014 Numerical Algorithms Group.
*
* Mark 24, 2013.
*/
#include <math.h>
#include <nag.h>
#include <nag_stdlib.h>
#include <nage01.h>

int main(void)
{
/* Scalars */
Integer  exit_status = 0;
Integer  d, i, j, liq, lrq, m, n, nq, nw;
/* Arrays */
double   *f = 0, *q = 0, *qx = 0, *rq = 0, *x = 0, *xe = 0;
Integer  *iq = 0;
/* NAG types */
NagError fail;

INIT_FAIL(fail);

printf("nag_nd_shep_interp (e01zmc) Example Program Results\n");

/* Skip heading in data file */
scanf("%*[^\n] ");

/* Input the number of nodes. */
scanf("%ld %ld%*[^\n] ", &d, &m);

liq = 2 * m + 1;
lrq = (d + 1) * (d + 2)/2 * m + 2 * d + 1;
if (!(x = NAG_ALLOC(d*m, double)) ||
!(f = NAG_ALLOC(m, double)) ||
!(iq = NAG_ALLOC(liq, Integer)) ||
!(rq = NAG_ALLOC(lrq, double)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}

/* Input the data points X and F. */
for (i = 0; i < m; i++)
{
for (j = 0; j < d; j++) scanf("%lf", &x[i*d+j]);
scanf("%lf%*[^\n] ", &f[i]);
}

/* Generate the interpolant using nag_nd_shep_interp (e01zmc):
Interpolating functions, modified Shepard's method, d variables. */
nq = 0;
nw = 0;
nag_nd_shep_interp(d, m, x, f, nw, nq, iq, rq, &fail);
if (fail.code != NE_NOERROR)
{
printf("Error from nag_nd_shep_interp (e01zmc).\n%s\n", fail.message);
exit_status = 1;
goto END;
}

/* Input the number of evaluation points and allocate arrays with lengths
based on this. */
scanf("%ld%*[^\n]", &n);
if (!(xe = NAG_ALLOC((d)*(n), double)) ||
!(q = NAG_ALLOC((n), double)) ||
!(qx = NAG_ALLOC((d)*(n), double)))
{
printf("Allocation failure\n");
exit_status = -2;
goto END;
}

/* Input the evaluation points.*/
for (i = 0; i < n; i++)
{
for (j = 0; j < d; j++) scanf("%lf", &xe[i*d+j]);
scanf("%*[^\n] ");
}
/* Evaluate the interpolant using nag_nd_shep_eval (e01znc), at given
* interpolated values, where interpolant previously computed by
* nag_nd_shep_interp (e01zmc). */
nag_nd_shep_eval(d, m, x, f, iq, rq, n, xe, q, qx, &fail);
if (fail.code != NE_NOERROR)
{
printf("Error from nag_nd_shep_eval (e01znc).\n%s\n", fail.message);
exit_status = 2;
goto END;
}

/* Print evaluation points and interpolated function values at those points */

printf("\n%4s|%39s%23s|%8s", "", "Interpolated Evaluation Points", "",
"Values\n");
printf(" ---|--------------------------------------------------------------");
printf("+-------\n");
printf(" i  |  ");
for (i = 0; i < d; i++) printf("XE(%1ld,i)   ", i);
printf("|  q[i]\n");
printf(" ---|--------------------------------------------------------------");
printf("--------\n");

/* Results */
for (i = 0; i < n; i++)
{
printf(" %1ld ", i);
for (j = 0; j < d; j++) printf("%10.4f", xe[i*d+j]);
printf("  %10.4f\n", q[i]);
}

END:
NAG_FREE(f);
NAG_FREE(q);
NAG_FREE(qx);
NAG_FREE(rq);
NAG_FREE(x);
NAG_FREE(xe);
NAG_FREE(iq);

return exit_status;
}
```