Facebook
From Swiftrage, 7 Years ago, written in Scheme.
Embed
Download Paste or View Raw
Hits: 268
  1. (define trapez
  2.   (lambda (f a b)
  3.     (* (/ (+ (f a) (f b)) 2) (- b a))
  4.   )
  5. )
  6.  
  7. (define x1
  8.   (lambda (a b n)
  9.     (+ a (/ (- b a) n))
  10.   )
  11. )  
  12.  
  13. (define calka-tmp
  14.   (lambda (f a b n w)
  15.      (if (> n 1)
  16.          (calka-tmp f
  17.                     (x1 a b n)
  18.                     b
  19.                     (- n 1)
  20.                     (+ w (trapez f a (x1 a b n)))
  21.           )
  22.          w
  23.       )
  24.   )
  25. )
  26.  
  27. (define calka
  28.   (lambda (f a b n)
  29.     (calka-tmp f a b n 0)
  30.   )
  31. )
  32.  
  33. (display (calka (lambda (x) (* x x)) 0. 1. 100000))