function TABLEAU=RKHeun()
  TABLEAU.ns = 2 ;
  TABLEAU.c  = [0,1] ;
  TABLEAU.b  = [1/2,1/2] ;
  TABLEAU.a  = [ 0 0 ; ...
                 1 0 ] ;
end