function RBF1plot( x,y,yout,mu,sigma )
% RBF1plot

[d,M] = size( mu );
plot( x, y, '.' )
hold on
plot( x, yout, 'r-' )
meany = mean(y);
rangey = max(y)-min(y);
for j=1:M
  y0 = meany + mod(j,2)*rangey/20;
  y1 = rangey/40;
  plot( mu(j)+[-sigma -sigma sigma sigma], [y0+y1 y0 y0 y0+y1] ,'g' )
end
xlabel( 'x' )
ylabel( 'y' )
hold off