% Douglas N. Arnold, dna@math.psu.edu, 11/96
% for use with ginterp
function runge(n)
global xmin xmax ymin ymax x y
xmin=-6;
xmax=6;
ymin=-.5;
ymax=1.5;
x=linspace(-5,5,n+1);
y=1./(1+x.^2);