NB. my misc NB. linear algebra require 'numeric plot' NB. Kimura P7 Example D0=: _2 1 1 1,_4 1 3 3,: _2 2 1r2 1 cr=: %. }:"1 NB. Cramer method rheval=:+/ .*&>/@|.&(128!:0) NB. Roger Hui rh=: [: (<0 1)&|: rheval ^:(1000) NB. ^:_ mseval=:(>&(1&{) +/ .* >&(0&{))&(128!:0) NB. M. Shimura ms=:[: (<0 1)&|: mseval ^:(100) NB. SVD NB. from Clausen NB. column robbery scam destroy NB. Oslow Northan Southtern area D1=:395 2456 1758 , 147 153 916 ,: 694 327 1347 NB. Beautiful JHSchool math D2=: 2 2 2 2,1 _1 1 _1,: _1 1 _1 1 D3=: 3 6 1 8,5 2 1 1,: 7 8 2 2 D4=:4 1 5,:3 2 1 NB. NB. Kimura p167 D5=: 2 _1 _1 1, _3 5 1 2,: 1 0 2 1 D6=: 2 1 _2 _1, 3 2 1 2,: _1 _1 1 1 NB. ============================ NB. spline data S0=: _10 0,_5 5,0 0,5 _5,: 10 0 S1=: |: /:~ |: ?. 2 10 $ 100 plot_spline=: 4 : 0 NB. draw 2 line pd 'reset' pd x pd y pd 'show' ) spline=: 3 : 0 N=: # y tmp0=.> N # < 1 4 1 ,N#0 (-i.10) |."0 1 > (N=. 10 )# < 1 ) NB. ============================== NB. keyPoint of Linear math /Satsuma&Yotsuya NB. Ex.1.2 K0=: 2 1 1 15,4 2 5 39,: 8 8 9 83 K1=: 2 1 1 15,4 6 3 41,: 2 9 3 37 NB. example of K0bad matrix NB. Satoh Nakamura p. 71 SN0=: 2 14 2 12 8,3 23 7 24 10,4 31 17 47 20,: 4 29 8 32 19 NB. Example for choleski C0=: 6 4 1 7,4 8 2 6,:1 2 1 8 C1=:1 0 1,0 2 0,: 1 0 3 NB. C1 is from MATLAB page NB. iteration method NB. Satoh Nakamura SN1=: 4 1 1 4,2 6 1 _8,: 1 2 5 7 T0=: 3 1 _1 0,1 _4 2 24,:2 _1 5 14 NB. NB. ========================== NB. SVD NB. wolflam Example W0=: 1 0 _1,: _2 1 4