function y=runge(x)
n=length(x);
y=ones(n,1)./(1+x.^2);