<jupyter_start><jupyter_text>Brusselator from [The Endeavour](https://www.johndcook.com/blog/2020/02/07/behold-the-brusselator/?utm_source=feedburner&utm_medium=email&utm_campaign=Feed%3A+TheEndeavour+%28The+Endeavour%29)<jupyter_code>from scipy import linspace
from scipy.integrate import solve_ivp
import matplotlib.pyplot as plt
A, B = 1,3
def brusselator(t, z):
  x, y = z
  return [A + x*x*y - (B+1)*x, B*x - x*x*y]
a, b = 0, 10
t = linspace(a, b, 1000)
for x0 in range (0,6):
  for y0 in [0, 3]:
             sol = solve_ivp(brusselator, [a,b], [x0, y0], t_eval=t)
             plt.plot(sol.y[0], sol.y[1], ":", color="tab:blue")
plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.17159763313609466, 'fraction_numerical': 0.04881656804733728, 'mean_word_length': 4.689075630252101, 'pattern_counts': {'":': 1, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '490372', 'license': 'MIT', 'n_tokens_mistral': 328, 'n_tokens_neox': 308, 'n_words': 58, 'path': 'Brusselator.ipynb', 'repo_name': 'supai-red/examples', 'size': '676'}