function res = RD_spots(ncells)
%

interval = 19999;
value_switch = 1;
stripe_flag = 0;

STRIPES = 1;
SPOTS = 2;
MAX = 200;
tri = zeros(MAX,3);
ngh = zeros(MAX*4,1);

p = zeros(MAX,MAX);
a = zeros(MAX,1);
b = zeros(MAX,1);
c = zeros(MAX,1);
da = zeros(MAX,1);
db = zeros(MAX,1);
speed = 1.0;
sim = 1;

xsize = 30;
ysize = 30;

beta_init = 12;
beta_rand = 0.1;
a_steady = 4.0;
b_steady = 4.0;
diff1 = 0.25;
diff2 = 0.0625;
p1 = 0.2;
p2 = 0.0;
p3 = 0.0;
sim = SPOTS;
value_switch = 2;

iterations = 20000;

for i = 1:xsize
    a(i) = a_steady;
    b(i) = b_steady;
    c(i) = beta_init + frand(-beta_rand,beta_rand);
end


Diff1 = 