function hsl=rgb2hsl(rgb_in)

%Usage
%       HSL = rgb2hsl(RGB)

zeroidx= (mx==mn);
S(zeroidx)=0;

lowlidx=L <= 0.5;
calc=(mx-mn)./(mx+mn);
idx=lowlidx & (~ zeroidx);
S(idx)=calc(idx);

hilidx=L > 0.5;
calc=(mx-mn)./(2-(mx+mn));
idx=hilidx & (~ zeroidx);
S(idx)=calc(idx);

hsv=rgb2hsv(rgb);
H=hsv(:,1);

hsl=[H, S, L];

hsl=round(hsl.*100000)./100000;
hsl=reshape(hsl, size(rgb_in));