#include <stdlib.h>
#include <stdio.h>
#define RND ((double)rand()/RAND_MAX)
#include "real.h"
#include "jacobi.c"
#include "rotate.c"

int main(int argc, char** argv){
	int n=(argc>1? atoi(argv[1]):10);

	real** A=(real**)calloc(n,sizeof(real*));
	for(int i=0;i<n;i++) A[i]=(real*)calloc(n,sizeof(real));

	for(int i=0;i<n;i++) for(int j=i;j<n;j++){
		A[i][j]=RND; A[j][i]=A[i][j];}

	real** V=(real**)calloc(n,sizeof(real*));
	for(int i=0;i<n;i++) V[i]=(real*)calloc(n,sizeof(real));

	int rotations=jacobi(A,n,V);

	printf("Ann=%g\n",A[n-1][n-1]);
	return 0;
}