sbiosimulate_wrapper(noRuns)

function data_end = sbiosimulate_wrapper(noRuns)

% Load simbiology project

sbioloadproject GeneRegulation m1 ;

data_end = zeros(noRuns, length(m1.species)) ;

% Performing Parameter Scan ...

for i=1:noRuns
    
    [t,x] = sbiosimulate(m1);             % Simulate the model
    data_end(i, :) = x(end,:) ;           % Find concentration@ t=1000
    
end