/* fichero formcuad.c */
/* Cálculo de una forma cuadrática */

#include <stdio.h>
#include <stdlib.h>

void main(void)
{
	double *x, **a;
	int    n, i, j;
	double fc(int n, double **a, double *x);
	
	printf("\nTeclee el número de filas y de columnas: ");
	scanf("%d", &n);
	
	a = calloc(n, sizeof(double*));
	for(i=0; i<n; i++)
		a[i] = calloc(n, sizeof(double));
	x = calloc(n, sizeof(double));
	
	for(i=0; i<n; i++)
		for(j=0; j<n; j++) {
			printf("\nTeclee a(%d, %d): ", i+1, j+1);
			scanf(" %lf", &a[i][j]);
		}
	for(i=0; i<n; i++) {
		printf("\nTeclee x(%d): ", i+1);
		scanf(" %lf", &x[i]);
	}	
	printf("\n\nLa forma cuadrática vale: %lf", fc(n, a, x));	
}

/* función para calcular el valor de la forma cuadrática */
double fc(int n, double **a, double *x)
{
	int    i, j;
	double fac=0.0;
	 
	for (i=0; i<n; i++)
		for (j=0; j<n; j++)
			fac += x[i]*a[i][j]*x[j];
	 
	return (fac);
}