function [self_action, disease_left_neighbor, disease_right_neighbor] = disease_outcome( disease )
	%% for a given disease (defined below) returns possible outcomes/actions
	% disease probabilities matrix
	% notice that probabilities of  self actions 
	% i.e. stay ill, die, and heal should add up to 1
	% probabilities distributed in a vector where each positions corresponds to
	% 1 stay ill  \
	% 2 die       | sum must be 1
	% 3 heal      /
	% 4 disease left neighbor (this one should be <= 1)
	% 5 disease left neighbor (this one should be <= 1)

	% normal dead member  (death is non contagious) 
	d(1,:)=[1.0, 0.0, 0.0,    0.0, 0.0];

	% weakly contagious but high mortality, hard to heal
	d(2,:)=[0.2, 0.8, 0.0,    0.1, 0.1];

	% highly contagious but low mortality, easy to heal
	d(3,:)=[0.8, 0.1, 0.1,    0.4, 0.4];

	% healthy alive member  (life is a disease too) 
	d(4,:)=[1.0, 0.0, 0.0,    0.0, 0.0];




	% ready to roll dices
	% 1st do we disease anyone
	% roll for dices for left neighbor
	p=rand();
	if ( p  <= d(disease, 4) )
		disease_left_neighbor=true;
	else
		disease_left_neighbor=false;
	end

	p=rand(); % reroll dices for right neighbor
	if( p <= d(disease, 5) )
		disease_right_neighbor=true;
	else
		disease_right_neighbor=false;
	end

	% 2nd what about yourself
	p=rand();

	if ( p <= d(disease, 1) )
		self_action =  disease; % keep what you have
	elseif ( p <= (d(disease, 2) + d(disease,1)) )
		self_action = 1; %death
	else
		self_action = 4; %healing,  recovery to normal
	end
	
end
			




