function du = lorenz(t,u)
% Indatafil för att lösa Lorenz kända system av differentialekvationer.
%
% x' = a*(y - x)
% y' = r*x - y -x*z
% z' = x*y - b*z
% där a, b och r är lämpligt valda parmetrar.
%
% Vi löser systemet genom att anropa Matlabfunktionen ode45 (i kommandofönstret)
% > [t,u] = ode45(@lorenz,[0 40], [20 30 60]);plot3(u(:,1),u(:,2),u(:,3))

a = 10; b = 8/3; r = 28;
x = u(1); y = u(2); z = u(3);

du = zeros(3,1); %Vi initierar utdata som en kolonnvektor av längd 3

du(1) = a*(y - x);
du(2) = r*x - y -x*z;
du(3) = x*y - b*z;