FFroehlich

/* -----------------------------------------------------------------
 * Adjoint sensitivity example problem.
 * The following is a simple example problem, with the coding
 * needed for its solution by CVODES. 
 *    dy/dt = -p*y
 * on the interval from t = 0.0 to t = 4.0, with initial
 * conditions: y = 1.0. The reaction rates are:
 * p=0.05 
 * 
 * -----------------------------------------------------------------*/

FFroehlich

% goal: plot barplots sorted according to their y value

xx =     {'C','C','C','B','B','B','A','A','A'};
yy =     [ 2 , 1 , 3 , 3 , 2 , 1 , 1 , 3 , 2 ];
type =   [ 1 , 2 , 3 , 1 , 2 , 3 , 1 , 2 , 3 ];
subset = type==1;

subx = xx(subset);
suby = yy(subset);
disp(['expected outcome' subx(order) ': ' num2str(suby(order))])

FFroehlich

% this .csv was exported from pokeminer
T = readtable('./pokeloc.csv');

spawnpoints = T{:,[6,7]};
spawnpoints = unique(spawnpoints,'rows');

spawnrates = zeros([size(spawnpoints,1),150]);
spawnnumber = zeros([size(spawnpoints,1),1]);