Numerical integration: Difference between revisions

Line 3,252:
(define (integrate f a b steps meth)
(define h (/ (- b a) steps))
(* h (for/sum ([i (+ steps 1)])
(meth f (+ a (* h i)) h))))
Line 3,276:
<lang racket>
CUBED
left-rect: 0.2550250000000000624502500000000005
mid-rect: 0.2601382512500000624998750000000006
right-rect: 0.2653280100000000625502500000000006
trapezium: 0.260176505250025
simpson: 0.260151002525
 
RECIPROCAL
left-rect: 4.6559810575146865499105751468
mid-rect: 4.605752058870831604762548678376
right-rect: 4.55797007838401955698105751468
trapezium: 4.606975567949346605986057514676
simpson: 4.606159895230326605170384957133
 
IDENTITY
left-rect: 1250000212499997.5
mid-rect: 1250000512500000.00000050
right-rect: 1250000712500002.5000015
trapezium: 1250000512500000.00000050
simpson: 1250000512500000.00000050
 
IDENTITY
left-rect: 1800000317999997.000000004
mid-rect: 1800000617999999.00000049999999993
right-rect: 1800000918000003.000001002000000004
trapezium: 1800000617999999.00000049999999993
simpson: 1800000617999999.00000049999999993
</lang>
 
Anonymous user