/* D03PCF Example Program C Version. NAG Copyright 1997. */ #include #include /* Include the appropriate NAG C Header File */ #include /* Global data */ double alpha; /* Routine for PDE initial condition */ void uinit(double u[20][2], double *x, int npts) { int i; for (i=0; i