code:

$M = 1000000;
$N_MAX = 10000;
$X_MAX = 1000;
$M_LEN = length($M)-1;
@a = (1*$M);
@ans = (1*$M);

for($n = 1; $n <= $N_MAX; $n++){
　　$tmp_mod = 0;
　　for($x = 0; $x <= $X_MAX; $x++){
　　　　$tmp_a = int(($tmp_mod*$M + $a[$x]) / $n);
　　　　$tmp_mod = ($tmp_mod*$M + $a[$x]) % $n;
　　　　$a[$x] = $tmp_a;
　　　　$ans[$x] = $ans[$x] + $a[$x];
　　}
}

$tmp_adv = 0;
for($x = $X_MAX; $x >= 0; $x--){
　　$ans[$x]+= $tmp_adv;
　　$tmp_adv = int($ans[$x] / $M);
　　$ans[$x] = $ans[$x] - $tmp_adv * $M;
}

print "e = ".$tmp_adv.".";
for($x = 0; $x <= $X_MAX; $x++){
　　printf("%0".$M_LEN."d", $ans[$x]);
}

exit;

ans.