
function [c, ceq] = consconsletter1(x)
global a b0 b1 b2 maxbird s a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 bfmax s e  area in g in2 aa % nonlinear inequality constraints
c = [x-1; %% Set yeild greater less than or equal to 1 for each parcel
      -x;%%% Set yeild positive
      
%       x(1)-x(2);
%       x(2)-x(3);
%       x(3)-x(4);
%       x(4)-x(5);
%       x(5)-x(6);
%       x(6)-x(7);
%       x(7)-x(8);
%       x(8)-x(9);
%       x(9)-x(10);            
      
  -sum(1-((x.^a(1))+(-in(1)*x.^g(1)))-in2(1))+(e*100); %%% constriant that say population over landscape is greater than e for species 1.
  -sum(1-((x.^a(2))+(-in(2)*x.^g(2)))-in2(2))+(e*100);
  -sum(1-((x.^a(3))+(-in(3)*x.^g(3)))-in2(3))+(e*100);
  -sum(1-((x.^a(4))+(-in(4)*x.^g(4)))-in2(4))+(e*100);
   -sum(1-((x.^a(5))+(-in(5)*x.^g(5)))-in2(5))+(e*100);
  -sum(1-((x.^a(6))+(-in(6)*x.^g(6)))-in2(6))+(e*100);
    -sum(1-((x.^a(7))+(-in(7)*x.^g(7)))-in2(7))+(e*100);
  -sum(1-((x.^a(8))+(-in(8)*x.^g(8)))-in2(8))+(e*100);
  -sum(1-((x.^a(9))+(-in(8)*x.^g(9)))-in2(9))+(e*100);

  
  ];
     

ceq=[];%%% No equality constriants in this program. For teh target models we set the production equality constraint here. 
 
end



































































































































































































































































































































































































































































































































































