布魯塞爾振子（英語：Brusselator）也稱布魯塞爾方程是一組模擬自催化反應的非線性微分方程：[1]
d ( x ( t ) ) d t = a − b ∗ x ( t ) + x ( t ) 2 ∗ y ( t ) − x ( t ) {\displaystyle {\frac {d(x(t))}{dt}}=a-b*x(t)+x(t)^{2}*y(t)-x(t)}
d ( y ( t ) ) d t = b ∗ x ( t ) − x ( t ) 2 ∗ y ( t ) {\displaystyle {\frac {d(y(t))}{dt}}=b*x(t)-x(t)^{2}*y(t)}
利用龍格－庫塔法可求得布魯塞爾振子的數值解，並利用Maple繪圖[1]
布魯塞爾振子模擬
布魯塞爾振子圖
布魯塞爾振子極限環場矢圖
布魯塞爾振子 3D 圖
