NB. ------------------------------- NB. Bradie Example 10.1 A0=: 2 _1 0 ,_1 2 _1 ,: 0 _1 2 jx1=: (i.4)*1r4 fx1=: 3 : '2 * 1 o. 2p1 * y' NB. 2sin(2 pi x) MATF=: (=/~i.3)- 1r4 * A0 NB. FTCS YX=: }. fx1 jx1 MATB=: (=/~i.3)+A0*1r4 NB. BTCS calc_ftcs=: 1 : 0 NB. y is lambda xi 'Lambda IX'=: y ) NB. ------------------------ loop_btcs=: 4 : 0 'MATB YX'=: y ANS=.: # IX) YX NB. calc system ) NB. Bradie Example 9.1 NB. ------------------------ MATB=: (=/~i.3)+A0*1r4 NB. BTCS loop_btcs=: 4 : 0 NB. Usage:( 1 x: ) 3 loop_btcs MATB;YX 'MATB YX'=. y ANS=.<'' NB. YX=. YX %. L:0 MATB for_ctr. i. x do. YX=. YX %. L:0 MATB ANS=. ANS,: +: Lamda TMP1=.;("1),.(SIZE - 2) #