Numerical integration: Difference between revisions
Add J
m (→TI-89 BASIC: use small integral sign) |
(Add J) |
||
Line 834:
*Main> intSimpson (\x -> x*x) 0 1 10
0.3333333333333334
=={{header|J}}==
'''Solution:'''
<lang j>integrate=: adverb define
'a b steps'=. 3{.y,128
size=. (b - a)%steps
size * +/ 2 u\ a + size * i.>:steps
)
simpson=: adverb define
(6 %~ [: +/ 4&*@u@-:@(+/) , u) y
)
trapezium=: adverb define
-: +/ u y
)
rectangle=: adverb define
u -: +/ y
)</lang>
'''Example usage'''
<lang j>
*: simpson integrate 0 1p1 10
10.3354255601
</lang>
=={{header|Java}}==
| |||