function P =  myBinocdf (x, n, p)

if isscalar(p)
	p = ones(size(n)) * p;
end
% Initialisiere P mit 1
P = ones (size(x));

% Bestimme die Wahrscheinlichkeiten für den Bereich, wo x < n
k = (x < n);
if (isscalar (n) && isscalar (p))
		P(k) = betainc (1 - p, n - x(k), x(k) + 1);
	else
		P(k) = betainc (1 - p(k), n(k) - x(k), x(k) + 1);
	end

end
