function y=fun(x);

y = 1./(1+x.^5);