function us=shiftp(typ,u)
//{\it decale le tableau u de +1 ou -1}
//{\it avec condition aux limites de Neumann}
select typ
   case '+1' then          //{\it shift +1}
           [m,n] = size(u) ;
           us(1:n-1) = u(2:n) ;
           us(n) = u(n-1) ;
   case '-1' then          //{\it shift -1}
           [m,n] = size(u) ;
           us(2:n) = u(1:n-1) ;
           us(1) = u(2) ;
   else 
           error('shift de +1 ou -1 ?')
end;

endfunction

//function us=shiftp(typ,u)
////{\it decale le tableau u de +1 ou -1}
////{\it avec condition aux limites periodiques}
//select typ
//   case '+1' then          //{\it shift +1}
//           [m,n] = size(u) ;
//           us(1:n-1) = u(2:n) ;
//           us(n) = u(1) ;
//   case '-1' then          //{\it shift -1}
//           [m,n] = size(u) ;
//           us(2:n) = u(1:n-1) ;
//           us(1) = u(n) ;
//   else 
//           error('shift de +1 ou -1 ?')
//end;