NB. J and complex number NB. Written by Giichiro Suzuki NB. JAPLA Nov./2008 NB. * complex 3 NB. * 2j1 min 3j4 NB. * 2j1 com_d 3j4 NB. --Suzuki------------------ complex=:3 :0 r=.j.^:(i.4)y+h=.0 while.h<<.-:y do.s=.j.^:(i.4)j./h,y-h=.h+1 r=.r,~.s,j.^:(i.4)j./(-h),y-h end. /:~r ) less=:3 :'/:~;complex&.><"0}.i.>:y' min=:[:<./[:+/@+.&>; int=:(2:=[:+/0:=]-<.)"1@+."0 com_d=:4 :0 Z=.less x min y (X=.(int x%Z)#Z);Y=.(int y%Z)#Z /:~~.((X e.Y)#X),(Y e.X)#Y ) NB. 左右の引数の公約数のセットを出力 NB. ----sombrero (JWiki)-------------------------- NB. ---function plot 3D------------------- NB. 3Dplot sombrero f=: 4 : '(cos r) %1 + r=. x +&: *: y' NB. plot _4 4 100;_4 4 100;'f' sombrero0=: [: (1&o. % ]) [: %: [: +/~ *: NB. 'surface' plot sombrero0 i:20j99 NB. 20j99 is 20 % 99 dyasombrero0=: (4 : '(1&o. % ]) %: +/ *: x , y')"0/ NB. plot _25 25 100; _15 15 100;'dyasombrero0' NB. ------------------------------------------ NB. Norman Thomson Vector Vol23 NR4 ln_nt=:(^.@{. , }.)@*. NB. log length and angle det=: -/ . * NB. determinant r21tom=: ,._1 1 &*@|. NB. matrix from real 2-list c21tom=: ,.+@|.@(*&1 _1) NB. matrix from complex 2-list times=: +/ . *&.> NB. -----my script-------------------------------- mk_cmat=:3 : '2 2 $ 1 _1 1 1 * 1 |. 2#( +.) y' NB. ----------------------------- fsr=: 4 : ' -: +/ y,x%y' fr=: -:@ +/@:( ],[ % ]) NB. -------------------- taltaria=: 3 : 0 'P0 Q0'=: y T0=: %: +/ (^&2 -: Q0) , ^&3 P0 % 3 T1=: -(-: Q0) + T0 T2=: -(-: Q0) - T0 +/(3 %: T1) , 3 %: T2 ) NB. ------------------------ sin_d=: 1&o.@(*&1r180p1) cos_d=: 2&o.@(*&1r180p1) degree_sub=: 3 : 0 TMP=. |. *. y (cos_d {. TMP),{: TMP ) NB. -----DATA thomson-------------------- 'I i j k'=: c21tom &.> R10=:1 0;0 1;0j1 0;0 0j1 NB. -----DATA 4 quadric-------------------- I2=: 4 4 $(_1 1 _1 1)(1 4 11 14)}16#0 J2=: 4 4 $ (_1 1 1 _1)(2 7 8 13)}16#0 K2=: 4 4 $ (_1 _1 1 1)( 3 6 9 12)}16#0