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