function area = auc(x,y)


% sortiere daten
[x,i] = sort(x);
y = y(i);

fast = 1;
if fast == 1
    xd   =  x(2:end) - x(1:end-1);
    ym   = (y(2:end) + y(1:end-1))/2;  % erster Teil ist Qudratisch + 
                                        % zweiter Teil hälfte vom Quadrat
    area = xd'*ym;
else
    % alternativer ansatz mit for schleife:
    area = 0;
    for i = 2:length(x)
        xstep = x(i) - x(i-1);
        ystep = y(i) - y(i-1);
        area = area + (xstep * y(i-1)) + (xstep * ystep / 2);
    end
end