function [ uprime ] = LorenzModel( t, u, p )
%This function implements the Lorenz attractor
%p=[a,b,r]

uprime = [0;0;0];

uprime(1) = -p(1)*u(1) + p(1) *u(2);
uprime(2) = p(3)*u(1) - u(2) -u(1)*u(3);
uprime(3) = -p(2)*u(3) + u(1)*u(2);



end