Tartalomjegyzék

Brussellator

TELJES MATLAB SCRIPT KIEGÉSZÍTŐ FÜGGVÉNYEKKEL

file:   d2017_09_20_brussellator.m
author: Peter Polcz <ppolcz@gmail.com>

Created on 2017. September 20.

Output:

┌d2017_09_20_brussellator
│   - Persistence for `d2017_09_20_brussellator` initialized [run ID: 9219, 2017.09.20. Wednesday, 21:38:13]
│   - Script `d2017_09_20_brussellator` backuped

syms t real
pcz_generateSymStateVector(2);

a = 1;
b = 3;

f_sym = [
    a + x1^2*x2 - (b + 1)*x1
    b*x1 - x1^2*x2
    ];

f_fh = matlabFunction(f_sym, 'vars', {t x});

xeq = [
    a
    b/a
    ]

[tt,xx] = ode45(f_fh, [0,100], [1;2.99]);

plot(xx(:,1),xx(:,2))

Output:

xeq =
     1
     3

End of the script.

Output:

└ 0.61835 [sec]

Published with MATLAB® R2015a, and sanitarized by BeautifulSoup (Python script using BSoup is written by Polcz).