::::::::::::::
tmp.m
::::::::::::::
clear;
A=1:10;
B=A.^2;
B;
C=fliplr(B);
C=C;
Y=[C 0 B];
X=-10:10
plot(X,Y)

::::::::::::::
tmp2.m
::::::::::::::
function result=tmp2(N)
total=d6(N)+d6(N);
result=hist(total,11);
result=result/N;


function rolls=d6(N)
R1=rand(1,N)*6;
rolls=ceil(R1);
::::::::::::::
tmp3.m
::::::::::::::
function r=tmp3(N)

count=1;
num=3;
array=[2]

while(count<N)
  if(min(mod(num, array)~=0))
     array=[array num];
     count=count+1;
  end
  num=num+1;
end
r=array;