T = 200000;
own_meinhardt_1;
X = euler(f, x0, T/20000, T);
plot3(X(1,:), X(2,:), X(3,:))