% Function to be optimized:
% Author: K. Passino, Version: 5/23/00
function J=nutrientsfunc1(theta,flag)

if flag==1  % Test to see if main program indicated nutrients
	J=0;
	return
end


% To test concentric pattern of groups of bacteria

J=-1000+1000*exp(-0.07*((theta(1,1)-15)^2+(theta(2,1)-15)^2));