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

#include "bnc.h"

int main()
{
	DVector x, y;
	DMatrix a;

	a = init_dmatrix(2, 2);
	x = init_dvector(2);
	y = init_dvector(2);

	set_dmatrix_ij(a, 0, 0, 1.0);
	set_dmatrix_ij(a, 0, 1, 2.0);
	set_dmatrix_ij(a, 1, 0, 3.0);
	set_dmatrix_ij(a, 1, 1, 4.0);

	set_dvector_i(x, 0, 1.0);
	set_dvector_i(x, 1, 2.0);

	print_dmatrix(a);
	print_dvector(x);

	mul_dmatrix_dvec(y, a, x);

	print_dvector(y);

	free_dmatrix(a);
	free_dvector(x);
	free_dvector(y);

	return EXIT_SUCCESS;
}