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

#include "bnc.h"

#define DIM 10

int main()
{
	long int i;
	double ip;
	DVector x, b;

	x = init_dvector(DIM);
	b = init_dvector(DIM);

	for(i = 0; i < DIM; i++)
	{
		set_dvector_i(x, i, (double)(i + 1));
		set_dvector_i(b, i, (double)(DIM - i));
	}

	ip = ip_dvector(x, b);

	printf("Inner Produce: %25.17e\n", ip);

	free_dvector(x);
	free_dvector(b);

	return EXIT_SUCCESS;
}