function y=MyTaylorPolynom(x,df,x0,n)

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

end
