convert the following matlab code to a python :

n = 1000; p = 1/216; q = 215/216; trials = zeros(16,1);

for x=0:15 result = nchoosek(n,x) * (p^x) * (q^(n-x)); trials(x+1) = result; end