SICP問題1.29
simpsonの公式
(define (my-simpson f a b n) (define (h) (/ (- b a) n)) (define (y k) (f (+ a (* k (h))))) (define (term k) (cond ((or (= k 0) (= k n)) (y k)) ((even? k) (* 2 (y k))) (else (* 4 (y k))))) (define (next i) (+ i 1)) (* (/ (h) 3.0) (sum term 0 next n)) )
実行結果。なんか胡散臭い
(my-simpson cube 0 1 100) ; => 0.25 (my-simpson cube 0 1 1000) ; => 0.25
これで良いのか?