- /
-
Storm Tide
on 10 Oct 2021
- 6
- 17
- 0
- 0
- 275
n=200;
[x,y]=meshgrid(1:n);
f=interp1([1 n],[.03 .06],x);
k=interp1([1 n],[.075 .04],y);
a=ones(n);
b=rand(n)>0.98;
m=1./[20 5 20;5 -1 5;20 5 20];
for i=1:500
A=a+(conv2(a,m,'same')-a.*b.*b+f.*(1-a));
B=b+(conv2(b,m,'same')/2+a.*b.*b-(k+f).*b);
a=A;b=B;
end
imagesc(a)
axis off
colormap(fliplr(hot))