sub log10{ log( $_[0] ) / log( 10 ) }

printf("%.16e\n", log10($_))
   for 1e2, 1e3, 1e4;

##</code><code>##

2.0000000000000000e+000
2.9999999999999996e+000
4.0000000000000000e+000