>function entropy (s) ...
$  v=strtochar(s);
$  m=getmultiplicities(unique(v),v);
$  m=m/sum(m);
$  return sum(-m*logbase(m,2))
$endfunction
>entropy("1223334444")
 1.84643934467