SimBiology first script -- modified to index names of species

20201002_simbio_first.m

% 2020-HS Intro Bio Computers
% SimBiology first script -- modified to index names of species

% The trick
index_of = containers.Map();
for i = 1 : length(BioComp_1.Species)
    index_of(BioComp_1.Species(i).Name) = i;
end

clear AB;
A_index = 1;

for A = 0:0.5:10
    B_index = 1;
    for B = 0:0.5:5
        BioComp_1.Species(index_of('A')).InitialAmount = A;
        BioComp_1.Species(index_of('B')).InitialAmount = B;
        [t, x] = sbiosimulate(BioComp_1);
        AB(A_index, B_index) = x(end, 3);
        B_index = B_index + 1;
    end
    A_index = A_index + 1;
end
figure;
surf(0:0.5:10, 0:0.5:5, AB');
xlabel('A initial');
ylabel('B initial');
zlabel('AB equilibrium');