首页 > 其他 > 详细

习题1.29 (积分方法的优化---simpson规则)

时间:2017-11-25 00:33:10      阅读:298      评论:0      收藏:0      [点我收藏+]
(define (sum term a next b)
  (if (> a b)
      0
      (+ (term a)
         (sum term (next a) next b))))

(define (integral f a b dx)  ;积分
  (define (add-dx x) (+ x dx))
  (* (sum f (+ a (/ dx 2.0)) add-dx b)
     dx))

  (define (cube x) (* x x x))
  
(integral cube 0 1 0.0001)



(define (simpson f a b n)
  
  (define h (/ (- b a) n))

(define (y k)
  (f (+ a (* k h))))

(define (factor k)
  (cond ((or (= k 0) (= k n)) 1)
        ((odd? k) 4)
        (else 2)))
  
(define (term k)
  (* (factor k) (y k)))

(define (next k)
  (+ k 1))

(if (not (even? n))
    (error "n can‘t be odd")
    (* (/ h 3)
       (sum term (exact->inexact 0) next n))))

(simpson cube 0 1 100)

  

0.24999999874993412
0.24999999999999992

  

习题1.29 (积分方法的优化---simpson规则)

原文:http://www.cnblogs.com/R4mble/p/7892731.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!