Morphogenesis (Mimicking)

on 29 Oct 2021

%-----------------------------------------------------

%% GRAY-SCOTT REACTION-DIFFUSION Feed-Kill

f=.055;k=.061 %

k = 0.0610

N=128; % gridsize (N=128)

M=64;

%% initial conditions

x=meshgrid(1:N);

u=ones(N);

v=0*u;

r=50:63;

v(r,r+5)=1; % center seeding

v(r-6,r-1)=1;

%% petri-dish filter

w=sqrt((x-M).^2)<2*M;

%j = 0;

while j<8e3

u=u+w.*((del2(u,2))-u.*v.^2+f.*(1-u));

v=v+w.*((del2(v,2))/2+u.*v.^2-(f+k).*v);

j=j+1;

end

% plotting

h=surf(x,x.',w.*u);

%h.CDataMapping='scaled';

colormap(turbo)

shading interp

axis off

view(45,90)

Load full remix tree