Numerical integration: Difference between revisions
Content added Content deleted
m (ββTI-89 BASIC: use small integral sign) |
(Add J) |
||
Line 834: | Line 834: | ||
*Main> intSimpson (\x -> x*x) 0 1 10 |
*Main> intSimpson (\x -> x*x) 0 1 10 |
||
0.3333333333333334 |
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}}== |
=={{header|Java}}== |