function [b,g]=makebspareshare9  %%% function to make differnt yeild density curves. "Hump" indicates intermediate species, Share indicates sharing species, spare indicates sparing species. Values fill this equation 1-((x.^a(1))+(-in(1)*x.^g(1)))-in2(1))
global b g in in2 aa
b=0;
g=0;
in=0;
bhump=[5.88834211835907;6.43475162245372;1.76192089776104;6.48025513683412;4.79415547735246;1.58524242999646;2.67098931320229;4.28128911522990;6.74504101260579;6.78933121119566;1.94567849006529;6.82355669056369;6.74300168945767;3.91225389233705;5.80168281333280;1.85131803176329;3.53056769575765;6.49441315113440;5.75324397735733;6.75695455835742;4.93444419493952;1.21427007144514;6.09477583521266;6.60395948654530;5.07241092914664;5.54644078347000;5.45879480874950;3.35336211720501;4.93286734106534;2.02712012686937];
bshare=[5.88834211835907;6.43475162245372;1.76192089776104;6.48025513683412;4.79415547735246;1.58524242999646;2.67098931320229;4.28128911522990;6.74504101260579;6.78933121119566;1.94567849006529;6.82355669056369;6.74300168945767;3.91225389233705;5.80168281333280;1.85131803176329;3.53056769575765;6.49441315113440;5.75324397735733;6.75695455835742;4.93444419493952;1.21427007144514;6.09477583521266;6.60395948654530;5.07241092914664;5.54644078347000;5.45879480874950;3.35336211720501;4.93286734106534;2.02712012686937];
bspare=[0.276935725388087;0.414985389386328;0.325975472178428;0.654440208531975;0.525959964012456;0.416493556356697;0.847745765106662;0.626737682037452;0.594751247462026;0.925474297446829;0.357255116938336;0.781480206199649;0.778356184850646;0.442401262277821;0.611039476652699;0.168268860606757;0.148555106799946;0.577717797708075;0.801250507091810;0.940609615806265;0.216915587626357;0.611941294784973;0.522451576952385;0.110711862551117;0.403410379958993;0.245964077373918;0.814856086615516;0.380093537840324;0.575679821955591;0.249083856549803];
bfifit=vertcat(bshare(1:5,1),bspare(1:4,1));
bmostshare=vertcat(bshare(1:8,1),bspare(1,1));
bmostspare=vertcat(bshare(1,1),bspare(1:8,1));
bsharesparehump=vertcat(bshare(1:3,1),bspare(1:3,1),bhump(1:3,1));
bsharehump=vertcat(bshare(1:5,1),bhump(1:4,1));
bsparehump=vertcat(bspare(1:5,1),bhump(1:4,1));
b=vertcat(bhump(1:9,1),bshare(1:9,1),bspare(1:9,1),bmostshare,bmostspare,bsharesparehump,bsharehump,bsparehump);

ghump=[0.0368268860606757;0.048555106799946;0.05717797708075;0.801250507091810;0.940609615806265;0.216915587626357;0.611941294784973;0.522451576952385;0.110711862551117;0.403410379958993;0.245964077373918;0.814856086615516;0.380093537840324;0.575679821955591;0.249083856549803;0.641783747261473;0.336674156086130;0.688671188629104;0.720293052826007;0.773336433541339;0.505487438652248;0.175439240197239;0.306079271845137;0.922003625351503;0.237140217072301;0.843235279740593;0.584508191734051;0.996521244964197;0.170357975877865;0.498410442797902];
gshare=zeros(9,1);
gspare=zeros(9,1);
gfifit=zeros(9,1);
gmostshare=zeros(9,1);
gmostspare=zeros(9,1);
gsharesparehump=vertcat(gshare(1:3,1),gspare(1:3,1),ghump(1:3,1));
gsharehump=vertcat(gshare(1:5,1),ghump(1:4,1));
gsparehump=vertcat(gspare(1:5,1),ghump(1:4,1));
g=vertcat(ghump(1:9,1),gshare,gspare,gmostshare,gmostspare,gsharesparehump,gsharehump,gsparehump);

ihump=ones(9,1);
ishare=zeros(9,1);
ispare=zeros(9,1);
ififit=zeros(9,1);
imostshare=zeros(9,1);
imostspare=zeros(9,1);
isharesparehump=vertcat(ishare(1:3,1),ispare(1:3,1),ihump(1:3,1));
isparehump=vertcat(ispare(1:5,1),ihump(1:4,1));
isharehump=vertcat(ishare(1:5,1),ihump(1:4,1));
in=vertcat(ihump,ishare,ispare,imostshare,imostspare,isharesparehump,isharehump,isparehump);

      v=ones(100,72);
   for j=1:1:72  
   for n=0.01:.01:1
    q=round(n*100);
    v(q,j)=1-((n.^b(j))+(-in(j)*(n.^g(j))))-in(j);
   end
   end
    m=max(v);
    in2=m'.*in;
    
aa=[0.717204491478464;0.717532166537288;0.718324567782782;0.722643319743731;0.725458682500445;0.735771320547784;0.742867278760385;0.750081888915348;0.751309478254026;0.758038015468520;0.759227367193345;0.764057116964345;0.766671534064359;0.770450590513484;0.779902604564297;0.781931882413222;0.789919265486780;0.800756287409354;0.809753481648948;0.810889792074482;0.818685896011558;0.833048040410344;0.838956632025114;0.842370147862913;0.843959206341231;0.854838817547240;0.857122710322430;0.857327048668507;0.857489340819000;0.878005523930996;0.886661371990248;0.891266981155529;0.900497831642498;0.901214103977678;0.920461989126686;0.937309129401156;0.945231907667531;0.954500055828284;0.958124434797750;0.959434902270077;0.963321983875662;0.972278625236152;0.975309296907959;0.982653024725164;0.987953284687636;0.993165384282148;0.993338646352100;0.993940833107334;0.996504361983562;1.00001346135412;1.00028297449291;1.00471497079667;1.01115696550632;1.01268149848240;1.01298990547857;1.02808343124338;1.04603290936881;1.04711503661406;1.05693764440517;1.06170585362925;1.06591998905354;1.07059983375307;1.07443605290421;1.08939489562741;1.09225440920114;1.09576315174498;1.10748132451945;1.10783677082640;1.10918314248944;1.11924749940088;1.12672946826021;1.12761668300735;1.12902224704042;1.13273649914904;1.13346375542011;1.13819851771327;1.14644455622048;1.16743103393877;1.17771032415113;1.18019834521144;1.18060877366184;1.18329365471781;1.19052825524757;1.19088913231577;1.19518827724123;1.21566538338773;1.23110080492149;1.23190715984566;1.23455350259847;1.24223233633379;1.24283334284042;1.24292858120794;1.24797209658354;1.25402778526195;1.25731248368683;1.27785312357215;1.28378473285832;1.28720838978470;1.28784902701365;1.29278920189698];



