%makmx(vrs)
% this is a routine to make the desert tortoise matrix from 
% a set of values for the vital rates. order of rates:
%WARNING: this function uses whatever the current values 
	%are of vrs to build the mx. 

%s2 s3 s4 s5 s6 s7 g2 g3 g4 g5 g6 f5 f6 f7
%1  2  3  4  5  6  7  8  9  10 11 12 13 14
dims = 8;
mx = zeros(dims,dims);
for ii=3:7
	mx(ii,ii) = vrs(ii-2)*(1-vrs(ii+4));
	mx(ii+1,ii)=vrs(ii-2)*vrs(ii+4);
end
mx(8,8) = vrs(6);
mx(2,2) = mx(3,3);
mx(3,2)= mx(4,3);
mx(2,1) = vrs(1);
mx(1,6) = vrs(12);
mx(1,7) = vrs(13);
mx(1,8) = vrs(14);