/* ------------------------------------ */
/*  Save as :  c0bn.c               */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A         = rsymmetric_mR(i_mR(r,r),999.);
double **EignVector=               i_mR(r,r);
double **T_EignVector=             i_mR(r,r);
double **EigsValue =               i_mR(r,C1);

double **T1=             i_mR(r,r);
double **T2=             i_mR(r,r);

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

  clrscrn();   
  eignVectoruv_mR(A,EignVector); 
  printf(" EignVector:");
  p_mR(EignVector,S10,P4,C6);
   
  transpose_mR(EignVector,T_EignVector);   
  mul_mR(T_EignVector,A,T1);
  mul_mR(T1,EignVector,T2);
  p_mR(T2,S10,P4,C6);  
   
  eigsuv_mR(A,EigsValue);
  printf(" EigsValue :");
  p_mR(EigsValue,S10,P4,C1);   
   
  f_mR(A);
  
  f_mR(EignVector);
  f_mR(T_EignVector);
  
  f_mR(EigsValue);  
     
  f_mR(T1);
  f_mR(T2);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(2)+R2);

} while(stop_w());

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