avec = [1 0 0];
bvec = [0 1 0];

pt1 = 5;
pt2 = 5;

na=5;
nb=5;
ctr = 1;
for i=1:na
        M(ctr,:)=avec*i;
        ctr = ctr + 1;
end