n=input('how many pairs of die do you want to roll?\n')
numtr=input('how many trials do you want')
A=ceil(rand(2*n,numtr)*6);
Occ=zeros(12,1);%number of times a 2,3,4,etc occurs in our rolls
for ii = 1:numtr
for jj =1:2:2*n
    sum=A(jj,ii)+A(jj+1,ii);
    Occ(sum)=Occ(sum)+1;
end
end
for kk=1:12
    Av(kk)=Occ(kk)/numtr;
    fprintf('tHE AVERAGE NUMBER OF TIMES %d occurs in %d rolls\n',kk,n)
    fprintf('is %f\n',Av(kk))
end