//  Programme qui resoud   -div(gradU) = graP dans un rectangle 
// de dimension lx et ly et des cl de Neuman sur le bord

clear  // detruit les anciennes variables
//
//  function scilab  que l'on va utiliser 
//
   getf('maillage_2D.sci');  getf('trace_maillage_2D.sci');  
   getf('gauss.sci');
   getf('cond_limit_2D.sci');
   getf('Q2.sci'); getf('Q1.sci');
   getf('assemblage.sci');   getf('cal_sol_lin.sci');
//
// pour tracer
//   
   getf('init_trace_2D.sci');  
   getf('trace_u_2D.sci'); 
//
// pour calculer le debit (Ecrit dans la seance3)
//
   getf('cal_debit_2D.sci');

//   calcul du systeme elementaire sur l'element 
//   (Ecrit en Seance 3)
//
    getf('cal_Melm_2D.sci'); 

//  on utilise des variables globales pour 
//    COOR, CONNEC, ADRESS, NUMER et  NI
//  et
//    XG et PG
//
   global CONNEC COOR NUMER ADRESS UC Nelt Nno Npt NI;
   global XG PG Ng
//
//  nbre de points pour l'integration numerique
   Ng = 4;
   gauss(Ng); 
//
   nx = 10;  ny = 10;  // nbre d'element dans les directions x et y
//
// maillage
   maillage_2D(1,1,nx,ny);
   init_trace_2D()
   trace_maillage_2D(nx,ny);
   cond_limit_2D();
//
//     Debut Cacul
//
// valeur initiale
   U = zeros(Npt,1); gradP = 1;  Q = 1.;
//
//  calcul de la solution pour gradP fixe
//
   U = cal_sol_lin(cal_Melm_2D,gradP,U);

   q = cal_debit_2D(U) 
   
   trace_u_2D(U,1);



      



