function TABLEAU=RK3()
  TABLEAU.ns = 3 ;
  TABLEAU.c  = [0 1/3 2/3] ;
  TABLEAU.b  = [1/4 0 3/4] ;
  TABLEAU.a  = [ 0   0   0 ; ...
                 1/3 0   0 ; ...
                 0   2/3 0 ] ;
end