Jで区分求積
area =: 1 : '+/ @ ((-~/@}: % {:) * u@}:@({. + (-~/@}: % {:) * i.@>:@{:))'
副詞として書いた。v area a b nで、vのグラフで[a,b]間をn等分して求積する。単純な左端型なので収束は遅い。台形近似とかも書きたい。
以下サンプル。
*: area 1 2 1000000 2.33333 ^. area 0.5 1.5 1000000 _0.0452293
area =: 1 : '+/ @ ((-~/@}: % {:) * u@}:@({. + (-~/@}: % {:) * i.@>:@{:))'
副詞として書いた。v area a b nで、vのグラフで[a,b]間をn等分して求積する。単純な左端型なので収束は遅い。台形近似とかも書きたい。
以下サンプル。
*: area 1 2 1000000 2.33333 ^. area 0.5 1.5 1000000 _0.0452293