辛普森的计划规则

时间:2013-09-19 20:54:29

标签: scheme rule

有人可以帮助我为辛普森的计划规则写一个方法的主体吗?以下是我到目前为止所做的,但我不知道是否应该递归地解决这个问题或使用不同的方法。

(define (sum-i term a b)
 (if (= a b)
    (term b)  
     (+ (term a) (sum-i term (+ a 1) b))))

(sum-i (lambda (x) x) 0 1)
(sum-i (lambda (x) (* x x)) 0 1)
(sum-i (lambda (x) (* x x x x)) 0 1)

(define (simpsons-rule f a b n)
   (let  ((h (/ (- b a) n)))
   (define (y_k k) (f (+ a (* k h))))

2 个答案:

答案 0 :(得分:1)

你的关闭我会告诉你你给出的函数和一个参数太多的情况,以及如何设置总和调用。

您可以递归设置,使用定义的递归函数而不是f,并在总和它的情况下调用它。

(define (simpsons-rule f a b n)
   (define h   (/ (- b a) n))
   (define Y_k (lambda (k) (f (+ a (* k h)))))
   (define (f x)
     (cond ((or (= x 0) (= x n)) ...)
           ((even? x) ...)
           (else ...)))
   (* (/ h 3)
      (sum f 0 plus1 n))) 

答案 1 :(得分:1)

#NOW我写了一个函数来计算simpons1 / 3rule的整数值:      #I想整合这个(x ^ 3-7x ^ 2 + 10)从0到100说

  ftn<-function(x){ 
  return((x^3-7*x^2+10))
 }
  simpons<-function(a,b,ftn,n){
  h<-(b-a)/n
   fa<-ftn(a)
   fb<-ftn(b)
  sum1<-0
  i<-1
  while(i<=n-1){
   sum1<-sum1+ftn(a+i*h)
 i<-i+2
    }
    j<-2
   sum2<-0
   while(j<=n-2){  
   sum2<-sum2+ftn(a+j*h)
   j<-j+2
   }
  s1<-h/3
  simp<-s1*(fa+fb+4*sum1+2*sum2)
   return(simp)
   }