00001 #include "tbl.h"
00002
00003 double* jacobi(int , double) ;
00004 double* pointsgausslobatto(int) ;
00005
00006 Tbl jacobipointsgl(int n) {
00007
00008 double* pointsgl = pointsgausslobatto(n) ;
00009
00010 Tbl jj(n+1,n+1) ;
00011 jj.set_etat_qcq() ;
00012
00013 int i,k ;
00014
00015 for (i = 0 ; i < n+1 ; i++ ) {
00016 double* yy = jacobi(n,pointsgl[i]) ;
00017 for (k = 0 ; k < n+1 ; k++ ) {
00018 jj.set(k,i) = yy[k] ;
00019 }
00020 delete [] yy ;
00021 }
00022 delete [] pointsgl ;
00023 return jj ;
00024 }
