function wdot=LorenzOneJet(t,w,sigma,r,b)

wdot=zeros(6,1);
wdot(1) = sigma*(w(2)-w(1));
wdot(2) = r*w(1)-w(2)-w(1)*w(3);
wdot(3) = w(1)*w(2)-b*w(3);
wdot(4) = sigma*(w(5)-w(4));
wdot(5) = r*w(4)-w(5)-w(1)*w(6)-w(4)*w(3);
wdot(6) = w(1)*w(5)+w(4)*w(2)-b*w(6);

