function y=runge(x)
[m,n]=size(x)
uno = ones(m,n);
y = uno./(1+x.^2);
