1;

function res=reg_grid1(f,xl,xr,dx)
  n=(xr-xl)/dx;
  dx=(xr-xl)/n;
  res.dx=dx;
  res.xl=xl;
  res.xr=xr;
  res.n=n;
  res.x=(0:(n-1))*dx+xl;
  y=res.x;
  for i=0:n-1
    y=f(res.x(i));
  endfor
  res.y=y;
endfunction
