## Copyright (C) 1998-2003 Joao Cardoso.
## 
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by the
## Free Software Foundation; either version 2 of the License, or (at your
## option) any later version.
## 
## This program is distributed in the hope that it will be useful, but
## WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
## General Public License for more details.
##
## This file is part of plplot_octave.

function p3

  global pl_automatic_replot;
  t = pl_automatic_replot;
  pl_automatic_replot = 0;

  xscale = 6.;
  yscale = 1.;
  xoff = 0.;
  yoff = 0.;

  m=60;
  x = xoff + xscale * ((1:m) + 1) / m;
  y = yoff + yscale * (x.^2);

  n=5;
  xs = x((1:n) * 10 + 3);
  ys = y((1:n) * 10 + 3);
  xe = (1:n)./5;
  yel = ys - (1:n);
  yer = ys + 2*(1:n);

  ## Plot the data points */
  title("#frPLplot Example 3 - y=x#u2 #d(with error bars)");
  plot(x, y, [xs' xe'], [ys' yel' yer'], '~g; ;');

  pl_automatic_replot = t;

endfunction
