%---------------------------------------------------------  
% ルンゲ−クッタ法による時間発展(上位ルーチン)
%---------------------------------------------------------
function [x,y]=rk4u(m,h,x,y,f)
  dx=h/m;
  for i=1:m
    [x,y]=rk4l(dx,x,y,f);
  end