function [z] = Ableitung (f,c,y0,a,b,n)


h = (b-a)/n;
x = a:h:b;

z = (f(x+h,c)-f(x-h,c))/(2*h);