NB. polymult.ijs 2010/5/25 NB. power of pm function NB. Original by Y. Nakano NB. Modified by T. Nishikawa pm =: +/ /. @(*/) npm =: 4 : 0 if. x = 0 do. 1 return. end. if. x = 1 do. y return. end. pm~ y if. x > 2 do. y (pm~^:(x-1)) y end. ) NB. 0 npm 1 2 NB. 1 NB. 1 npm 1 2 NB. 1 2 NB. 2 npm 1 2 NB. 1 4 4 NB. 3 npm 1 2 NB. 1 6 12 8 NB. 4 npm 1 2 NB. 1 8 24 32 16 NB. 5 npm 1 2 NB. 1 10 40 80 80 32 NB. ------Nakano--------------------------------------- NB. power of pm function (polynomials multiplication) wr=: 1!:2&2 bm=: 3 : 0 : if. x = 0 do. 1 return. end. im=.x i=. 1 ypm =. y *1 if. i = im do. goto_e. end. label_1. i=.i+1 ypm =. ypm pm y if. i = im do. goto_e. end. goto_1. label_e. NB. wr i ypm ) NB. ------------------------------------- NB. Lang変換チェック by Y.NAKANO NB. 2010/Jun/20 langchk0=: 3 : 0 : s=. x wr dx=.0 { x wr X =. }. x wr y wr ym=. # y i=.0 wr ' i= ', ":i wr p=. i{X while. i < ym do. s=.p wr a=.i{y if. (p = a) *(i>1) do. goto_e. end. wr t=. (i+1){X p=.t + (s -a)*dx NB. if. p = a do. goto_e. end. i=.i+1 wr ' i= ', ":i NB. if. p = a do. goto_e. end. end. label_e. wr ' Happy End !' p-a )