::::::::::::::
max3.m
::::::::::::::
function r=max3(a,b,c)
r=max([a b c])*3;
::::::::::::::
close.m
::::::::::::::
function r=close(v,s)
new=abs(v-s);
t=min(new);
i=min(find(t==new));
r=v(i);

::::::::::::::
cumulative.m
::::::::::::::
function y=cumulative(x)
[r c]=size(x);
for i=1:c
  y(i)=sum(x(1:i));
end
::::::::::::::
initM.m
::::::::::::::
function r=initM()

r(1:10, 1:10)=0;
r([1 10],:)=6;
r(:,[1 10])=6;
r(1:5,6:10)=2;
r(6:10,1:5)=4;
r(1:2,1:2)=5;