function lipo = lipomodelode_abs(t,Y,par)
par(2)=abs(par(2));
lipo=Y;
% from week1 number10 slides but creating bimodality by using absolute
% value parameters
lipo(1) = -(par(1)+par(2))*Y(1)+par(2)*Y(2);
lipo(2) =  par(2)*(Y(1)-Y(2));
