#include <stdio.h>

double * gradiente(double (*pf) (double * r, int n), double *r, int n, double eps);
/* esta función es un hiperparaboloide (una parábola en R^n, bah...) o la 
distancia al origen al cuadrado, como más les guste */

double parab (double *r, int n);

int main()
{
	double *r, *g;
	int i,n;
	
	printf("Ingrese la dimension del dominio (>0)\n");
	scanf("%d",&n);
	if (n<1)
	{
		printf("La dimensión debe ser mayor que 0\n");
		exit(1);
	}
	
	if ((r = (double *) malloc(n*sizeof(double)))==NULL)
	{
		perror("Error allocando espacio para el punto r\n");
		exit(1);
	}
	
	/* inicializo el punto donde quiero calcular el gradiente */
	for (i=0; i<n; i++)
		r[i] = i;
		
	/* calculo el gradiente de la función parab */	
	g = gradiente(parab, r, n, 0.00001);
	if (g == NULL)
	{
		perror("Error allocando espacio para el gradiente g\n");
		exit(1);		
	}
	
/* imprimo el valor de la función, sólo para verlo... */
	printf("f(r,n) = %lf\n", parab(r,n));

/* ahora imprimo el gradiente y el punto donde lo calculé */
	for (i=0; i<n; i++)
		printf("r[%d] = %lf, g[%d] = %lf\n",i, r[i],i,g[i]);
		
	return 0;	
}

double parab (double *r, int n)
{
	double res = 0.;
	int i;
	

	for (i=0; i<n; i++)
		res += r[i]*r[i];
	
	return res;
}