function [x,y] = myAuswertung(f,a,b,n)

h=(b-a)/n;
x=zeros(n+1,1);
y=zeros(n+1,1);
for i=1:n+1
    x(i)=a+(i-1)*h;
    y(i)=f(x(i));
end