clear all
close all

s=10/7;
a=-5; b=5;
x=a:s:b;

y=1./(1+x.^2), x=x', y=y'

[n,m]=size(x)

z=a:0.01:b
egy=interp1(x,y,z,'pchip')

egymin=min(egy);
egymax=max(egy);

axis([a,b,egymin,egymax]);
hold on
plot(x,y,'o')
plot(z,egy)
hold off


