n=input("n=");
p=input("p=");
N=input("N=");
M=0;
for j=1:N
X=0;
for k=1:n
  x=floor(rand()+p);
  X=X+x;
end;
M=M+X;
end ;
m=M/N;
disp(m);
