/* ------------------------------------ */
/*  Save as :  c0a6.c               */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(void)
{
double ab[R6*C6]={
+754,+750,-983,-721,-737,+394,
+750,+886,+930,+794,-809,-245,
-983,+930,+990,-747,-37,-647,
-721,+794,-747,+886,-301,+682,
-737,-809,-37,-301,-491,+816,
+394,-245,-647,+682,+816,-33
};

int r = RC6;
	
double **A         = ca_A_mR(ab,i_mR(r,r));
double **EignVector=            i_mR(r,r);
double **EigsValue =            i_mR(r,C1);

  clrscrn();
  printf(" Copy/Past into the octave windows \n\n");
  p_octave_mR(A,"a",S1,P0);
  printf(" [V, E] = eigs (a,10) \n\n");
  stop();

  clrscrn();   
  eignVectoruv_mR(A,EignVector); 
  printf(" EignVector:");
  p_mR(EignVector,S5,P6,C6);
   
  eigsuv_mR(A,EigsValue);
  printf(" EigsValue :");
  p_mR(EigsValue,S13,P6,C1);   
  stop();
     
  f_mR(A);
  f_mR(EignVector);
  f_mR(EigsValue);
}
/* ------------------------------------ */
int main(void)
{

  fun();

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */