require 'numeric trig plot' NB. Brain Bradie [Numerical Analysis] 2006 NB. Example 11.7 NB. ----------------------- NB. * partial_diff_wav0 5 NB. verb type NB. * (f0 pdiff_wav0_conj g0) 5 NB. conjunction type NB. -----------3 : 0 type----------- partial_diff_wav0=: 3 : 0 NB. u 5 STR=: (% y)* i.>: y W0=: sin STR W01=: wav_sub0 STR W1=: (- 1{W0), (}. }: W01), _2{ W0 TMP=. (: ctr){W0 TMP=. TMP,: y W0=: u STR NB. f0 W01=: (u wav_sub0_conj v) STR NB. f1 W1=: (- 1{W0), (}. }: W01), _2{ W0 TMP=. (: ctr){W0 NB. f2 TMP=. TMP,{: y NB. W(n) WN0=: }.@}: >_2{ y NB. W(n-1) TMP2=. }.@}: +//. (}. WN) ,: }: WN NB.W(n)j-1 +W(n)j+1 TMP2-WN0 NB. W(n)j+- - W(n-1)j )