loom

1
float sigma;
2
float R;
3
float b;
4

5
float x;
6
float y;
7
float z;
8

9
float __dt_x;
10
float __dt_y;
11
float __dt_z;
12

13
func __ODE[*]()
14
{
15
    __dt_x = sigma*(y-x);
16
    __dt_y = R*x - y - x*z;
17
    __dt_z = -b*z + x*y;
18
}
19