テスト勉強中に教科書を見ていて、f(x+dx)-f(dx)/dxだとx+dx/2近傍になっちゃうことに気がついた。 f(x+dx)-f(x-dx)/2dxとしないとだめなのね。

sub dx(){ 1.0e-4 } =head3 integral integral( func, from, to ) funcを積分範囲fromからtoで積分 =cut sub integral{ my( $func, $f, $t ) = @_; my $x = $f; my $r = 0; while( $x <= $t ){ # r = Sigma(x=f->t)f(x)dx $r += $func->($x); $x += dx; } re…