Add reaction-diffusion in BQN, with video generation script

grayscott.bqn

@@ -0,0 +1,31 @@
+⟨PBM,PPM⟩←•Import"../bqn-graphics/pnm.bqn"
+⟨Inferno,Magma,Viridis⟩←•Import"../bqn-graphics/colormaps.bqn"
+
+# Parameters
+dt‿da‿db‿f‿k←⟨1.0,0.16,0.08,0.060,0.062⟩
+n←500  # image size
+m←10  # seed size at the center of the image
+
+# Discrete Laplacian
+Lapl←(+´∘∾1‿¯1(⌽⋈⌽˘)¨<)-4⊸×
+# Gray-Scott update
+UpdateA←{𝕨+dt×(da×Lapl 𝕨)+(-𝕨×𝕩⋆2)+f×1-𝕨}
+UpdateB←{𝕩+dt×(db×Lapl 𝕩)+(𝕨×𝕩⋆2)-𝕩×k+f}
+Update←(UpdateA⋈UpdateB)´
+
+# Initialization
+a←n‿n⥊1
+Mask←{𝕩‿𝕩↑(-𝕨‿𝕨+𝕩÷2)↑(2×𝕨‿𝕨)⥊1}
+b←(m Mask n)⌈0.1×n‿n•rand.Range 0
+
+# Run
+Step←{stepsize 𝕊 i‿a‿b:
+  •term.OutRaw"."
+  •term.Flush @
+  ap‿bp←Update⍟stepsize a‿b
+  ("out/step"∾(¯5↑∾•Fmt¨5-⊸↑i)∾".pnm")•file.Chars (PPM ⌊255×>Inferno¨1-a)∾@+10
+  (i+1)‿ap‿bp
+}
+
+50 Step•_while_(600⊸>∘⊑⊢) 0‿a‿b
+•Out""