function deriv=diffin(x,delta)

fmais=fl(1/fl(fl(x)+fl(delta)));
fmenos=fl(1/fl(fl(x)-fl(delta)));
deriv=fl(fl(fmais-fmenos)/fl(2*fl(delta)));

end