function y=MyTaylorIntegral(df,x0,n,a,b)

y=0;
for k=0:n
    y = y + df(x0,k)/factorial(k+1)*((b-x0)^(k+1)-(a-x0)^(k+1));
end

end
