function d = derivata2(f,x,h)

d = (feval(f,(x+h))-feval(f,(x)))/h;