Reaction-Diffusion Hellmouth

on 29 Oct 2021

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

%% GRAY-SCOTT REACTION-DIFFUSION Feed-Kill

f=.055;k=.061;

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<1.5e3

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(w.*u);

colormap(flip(hot))

shading interp

axis off

view(37,75)

set(gcf,Color="k")

Load full remix tree