/* Driver for routine df1dim */

#include <stdio.h>
#include "nr.h"
#include "nrutil.h"

#define NDIM 3

int ncom;   /* defining declarations */
float *pcom,*xicom;
void (*nrdfun)();

void dfunc(x,df)
float df[],x[];
{
	int i;

	for (i=1;i<=3;i++) df[i]=(x[i]-1.0)*(x[i]-1.0);
}

main()
{
	ncom=NDIM;
	pcom=vector(1,ncom);
	xicom=vector(1,ncom);
	nrdfun=dfunc;
	printf("\nEnter vector direction along which to\n");
	printf("plot the function. Minimum is in the\n");
	printf("direction 1.0 1.0 1.0 - enter x y z:\n\n");
	pcom[1]=pcom[2]=pcom[3]=0.0;
	scanf("%f %f %f",&xicom[1],&xicom[2],&xicom[3]);
	scrsho(df1dim);
	free_vector(xicom,1,ncom);
	free_vector(pcom,1,ncom);
	return 0;
}