Mercurial > hg > jgplsrc

annotate test/g200.ijs @ 0:e0bbaa717f41 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
lol J
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Mon, 25 Nov 2013 11:56:30 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
1 NB. ^y ------------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
3 x =: 0.001 * _900+?4 5$2000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
4 y =: 0.001 * _900+?4 5$2000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
5 a =: x j. (*y)*(o.1)|y
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
6 a =: a,0,*/\4$0j1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
7
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
8 exp=: (%!i._50) & (+/ .*) @ ((i._50)&(^~/))
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
10 5e_13 > | (^a)%~(^a)-exp a
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
11 a =!.5e_13 ^.^a
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
12 a =!.5e_13 ^^.a
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
13
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
14 1 -: ^ 0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
15 2.71828182845904523536 -: ^ 1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
16 7.38905609893065022723 -: ^ 2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
17 148.4131591025766 -: ^ 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
18
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
19 rou =: [:^ 0j2p1&% * i. NB. roots of unity
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
20 ([: */ 1: = rou ^ ])"0 i.4 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
21
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
22 0 < ^ _744.440 NB. _1074*^.2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
23 0 = ^ _744.441
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
24 _ > ^ 709.782 NB. 1024*^.2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
25 _ = ^ 709.783
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
26
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
27 0 = ^ __
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
28 0 = ^ __j5 __j_5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
29
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
30 'domain error' -: ^ etx 'abc'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
31 'domain error' -: ^ etx <'abc'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
32
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
33 0 = ^_744.441j2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
34 0 = ^_744.441j_2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
35
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
36 'limit error' -: ^ etx _744.440j2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
37 'limit error' -: ^ etx _744.440j_2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
38
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
39 (j./ ^710+^. 2 1 o.0.01) = ^ 710j0.01
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
40 (j./ ^710+^. 2 1 o.1.57) = ^ 710j1.57
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
41 (j./_1 1 * ^710+^.|2 1 o.1.58) = ^ 710j1.58
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
42 (j./_1 _1 * ^710+^.|2 1 o.4.71) = ^ 710j4.71
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
43 (j./ 1 _1 * ^710+^. 2 1 o.1.57) = ^ 710j_1.57
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
44 (j./ 1 _1 * ^710+^. 2 1 o.0.01) = ^ 710j_0.01
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
45
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
46
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
47 NB. x^y -----------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
48
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
49 pow =: ^.@[ ^@* ]
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
50
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
51 (x=:x+0=x=:_50+?4$100) (^-:pow) y=: _50+?4 10$100
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
52 (x=:x+0=x=:_50+?4 10 2$100) (^-:pow) y=: 0.02*_50+?4 10$100
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
53 (x=:x+0=x=:o._50+?4 10$100) (^-:pow) y=:o. 0.01*_50+?4 10$100
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
54 (x=:x+0=x=:r./_50+?2 4$100) (^-:pow) y=:j./0.01*_50+?2 4 10$100
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
55
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
56 (a^ i.30) -: */\1,29$ a=: 2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
57 (a^-i.30) -: */\1,29$%a=: 2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
58 (a^ i.30) -: */\1,29$ a=:_2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
59 (a^-i.30) -: */\1,29$%a=:_2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
60 (a^ i.30) -: */\1,29$ a=: 2.45
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
61 (a^-i.30) -: */\1,29$%a=: 2.45
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
62 (a^ i.30) -: */\1,29$ a=:_2.45
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
63 (a^-i.30) -: */\1,29$%a=:_2.45
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
64 (a^ i.30) -: */\1,29$ a=: 2j1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
65 (a^-i.30) -: */\1,29$%a=: 2j1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
66 (a^ i.30) -: */\1,29$ a=: 2j_1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
67 (a^-i.30) -: */\1,29$%a=: 2j_1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
68 (a^ i.30) -: */\1,29$ a=:_2j1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
69 (a^-i.30) -: */\1,29$%a=:_2j1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
70 (a^ i.30) -: */\1,29$ a=:_2j_1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
71 (a^-i.30) -: */\1,29$%a=:_2j_1.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
72
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
73 a -: a^1 [ a=: ?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
74 a -: a^1 [ a=:-?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
75 a -: a^1 [ a=:o._1e9+?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
76 a -: a^1 [ a=:j./_1e9+?2 20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
77
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
78 (20$1) -: a^0 [ a=: ?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
79 (20$1) -: a^0 [ a=:-?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
80 (20$1) -: a^0 [ a=:o._1e9+?20$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
81 (20$1) -: a^0 [ a=:r./_1e7+?2 20$2e7
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
82
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
83 1 0 1 1 -: 0 0 1 1 ^ 0 1 0 1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
84
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
85 1 -: 0^0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
86 (9$1) -: (-:_5+i.9)^0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
87
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
88 _ _ _ _ 1 0 0 0 0 -: 0 ^ _4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
89 _ _ _ _ 1 0 0 0 0 -: ({.0 4 )^ _4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
90 _ _ _ _ 1 0 0 0 0 -: ({.0 4.5)^ _4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
91 _ _ _ _ 1 0 0 0 0 -: ({.0 4j5)^ _4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
92
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
93 _ _ _ _ 1 0 0 0 0 -: 0 ^ -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
94 _ _ _ _ 1 0 0 0 0 -: ({.0 4 )^ -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
95 _ _ _ _ 1 0 0 0 0 -: ({.0 4.5)^ -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
96 _ _ _ _ 1 0 0 0 0 -: ({.0 4j5)^ -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
97
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
98 _ _ _ _ 1 0 0 0 0 -: 0 ^ }. 0j1, -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
99 _ _ _ _ 1 0 0 0 0 -: ({.0 4 )^ }. 0j1, -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
100 _ _ _ _ 1 0 0 0 0 -: ({.0 4.5)^ }. 0j1, -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
101 _ _ _ _ 1 0 0 0 0 -: ({.0 4j5)^ }. 0j1, -:_4+i.9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
102
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
103 a=:1.41421356237309504880
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
104 a -: 2 ^ 0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
105 (%a) -: 2 ^_0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
106 (%a) -: 0.5 ^ 0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
107 a -: 0.5 ^_0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
108 (0, a) -: +. _2 ^ 0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
109 (0,%-a) -: +. _2 ^_0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
110 (0,% a) -: +. _0.5 ^ 0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
111 (0, -a) -: +. _0.5 ^_0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
112
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
113 1.25992104989487316476 -: 2^%3
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
114
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
115 x=: j./_50+2 30 ?@$ 100
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
116 n=: 30 ?@$ 50
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
117 x (^ -: pow) n
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
118 x (^ -: pow) -n
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
119
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
120 'domain error' -: 3 ^ etx 'abc'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
121 'domain error' -: 3 ^~etx 'abc'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
122 'domain error' -: 2 ^ etx 2;4 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
123 'domain error' -: 2 ^~etx 2;4 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
124
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
125 'length error' -: 2 3 ^ etx 4 5 6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
126 'length error' -: 2 3 ^~etx 4 5 6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
127 'length error' -: 2 3 ^ etx i.4 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
128 'length error' -: 2 3 ^~etx i.4 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
129
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
130
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
131 NB. x^0 -----------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
132
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
133 e.&'1 ' ": x ^ 0 [ x=: ?200$2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
134 e.&'1 ' ": x ^ 0 [ x=: _5e6+?200$1e7
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
135 e.&'1 ' ": x ^ 0 [ x=: 0.001*_5e6+?200$1e7
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
136 e.&'1 ' ": x ^ 0 [ x=: j./0.001*_5e6+?2 200$1e7
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
137
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
138
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
139 NB. 0^y -----------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
140
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
141 e.&'0 ' ": 0 ^ 10$1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
142 e.&'0 ' ": 0 ^ 1+?10$1e6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
143 e.&'0 ' ": 0 ^ 0.1*1+?10$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
144 e.&'0 ' ": 0 ^ j./ 1+?2 10$2e3
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
145
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
146 e.&'0 ' ": (2-2) ^ 10$1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
147 e.&'0 ' ": (2-2) ^ 1+?10$1e6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
148 e.&'0 ' ": (2-2) ^ 0.1*1+?10$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
149 e.&'0 ' ": (2-2) ^ j./ 1+?2 10$2e3
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
150
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
151 e.&'0 ' ": (o.0) ^ 10$1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
152 e.&'0 ' ": (o.0) ^ 1+?10$1e6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
153 e.&'0 ' ": (o.0) ^ 0.1*1+?10$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
154 e.&'0 ' ": (o.0) ^ j./ 1+?2 10$2e3
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
155
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
156 e.&'0 ' ": (-~2j1) ^ 10$1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
157 e.&'0 ' ": (-~2j1) ^ 1+?10$1e6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
158 e.&'0 ' ": (-~2j1) ^ 0.1*1+?10$2e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
159 e.&'0 ' ": (-~2j1) ^ j./ 1+?2 10$2e3
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
160
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
161 (,'_') -: ": 0^_5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
162 (,'_') -: ": 0^_5.4
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
163
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
164 (,'0') -: ": 0^3j4
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
165 '_ _ 0 0 0 0 0 0' -: ": 0^_3 _1 0j3 0j_3 3j4 3j_4 _3j4 _3j_4
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
166
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
167
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
168 NB. 0^0 -----------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
169
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
170 (,'1') -: ": ^~0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
171 (,'1') -: ": ^~-~2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
172 (,'1') -: ": ^~-~2.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
173 (,'1') -: ": ^~-~2j5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
174
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
175
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
176 NB. x^0.5 ---------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
177
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
178 x -: *~ x^0.5 [ x=:?2000$2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
179 (%:x) -: x^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
180 (j.x^0.5) -: (-x)^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
181
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
182 x -: *~ x^0.5 [ x=:?2000$10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
183 (%:x) -: x^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
184 (j.x^0.5) -: (-x)^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
185
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
186 x -: *~ x^0.5 [ x=:0.01*?2000$10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
187 (%:x) -: x^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
188 (j.x^0.5) -: (-x)^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
189
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
190 x -: *~ x^0.5 [ x=:j./0.01*_5000+?2 2000$10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
191 (%:x) -: x^0.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
192
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
193
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
194 NB. x ^!.s y ------------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
195
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
196 (3 ^!.4 [ 6) -: */3+4*i.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
197 (3j4^!._1 [ 6) -: */3j4-i.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
198
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
199 f =: 4 : '*/x+(i.y)*/s'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
200
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
201 (12 ^!.s 6) -: 12 f 6 [ s=:1.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
202 (12 ^!.s 6) -: 12 f 6 [ s=:1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
203 (12 ^!.s 6) -: 12 f 6 [ s=:_1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
204
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
205 'domain error' -: ^!.4 etx 5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
206 'domain error' -: 3^!.4 etx 5.6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
207 'domain error' -: ex '^!.''a'''
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
208 'domain error' -: ex '^!.(<4) '
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
209
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
210
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
211 NB. x m&|@^ y -----------------------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
212
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
213 6 = 2 (10&|@^) 1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
214 2 = 2 (10&|@^) 1+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
215 4 = 2 (10&|@^) 2+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
216 8 = 2 (10&|@^) 3+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
217
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
218 6 = 2x (10&|@^) 1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
219 2 = 2x (10&|@^) 1+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
220 4 = 2x (10&|@^) 2+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
221 8 = 2x (10&|@^) 3+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
222
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
223 6 = 2 (10&|@^) 10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
224 2 = 2 (10&|@^) 1+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
225 4 = 2 (10&|@^) 2+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
226 8 = 2 (10&|@^) 3+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
227
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
228 6 = 2x (10&|@^) 10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
229 2 = 2x (10&|@^) 1+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
230 4 = 2x (10&|@^) 2+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
231 8 = 2x (10&|@^) 3+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
232
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
233 6 = 10&|@(2 &^) 1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
234 2 = 10&|@(2 &^) 1+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
235 4 = 10&|@(2 &^) 2+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
236 8 = 10&|@(2 &^) 3+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
237
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
238 6 = 10&|@(2x&^) 1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
239 2 = 10&|@(2x&^) 1+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
240 4 = 10&|@(2x&^) 2+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
241 8 = 10&|@(2x&^) 3+1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
242
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
243 6 = 10&|@(2 &^) 10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
244 2 = 10&|@(2 &^) 1+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
245 4 = 10&|@(2 &^) 2+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
246 8 = 10&|@(2 &^) 3+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
247
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
248 6 = 10&|@(2x&^) 10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
249 2 = 10&|@(2x&^) 1+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
250 4 = 10&|@(2x&^) 2+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
251 8 = 10&|@(2x&^) 3+10^40x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
252
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
253 f=: 3 : 0
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
254 assert. 2 (y&|@^ -: y&|@^&x:) 10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
255 assert. 3 (y&|@^ -: y&|@^&x:) 10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
256 assert. 5 (y&|@^ -: y&|@^&x:) 10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
257 assert. 64 (y&|@^ -: y&|@^&x:) 10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
258 assert. 17393 (y&|@^ -: y&|@^&x:) 10000
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
259 1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
260 )
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
261
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
262 f 1e9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
263 f"0 (i:2)+<.%:<:2^IF64{31 63
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
264 f"0 (i:2)+<.%:<:2^53
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
265
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
266 x=: 7700892415753674751x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
267 1 = 2 x&|@^ -:x-1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
268
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
269 (3 (17 &|)@^ 9.5) -: 17|3 ^9.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
270 (3 (17x&|)@^ 9.5) -: 17|3 ^9.5
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
271 (3.7 (17 &|)@^ 9 ) -: 17|3.7^9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
272 (3.7 (17x&|)@^ 9 ) -: 17|3.7^9
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
273
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
274 h=: 7927
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
275 y=: 1e5 4 2
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
276 ((895 h&|@^ y),h|895 ^_1) -: 895 h&|@^ y,_1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
277 ((895x h&|@^ y),h|895x^_1) -: 895x h&|@^ y,_1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
278
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
279 (10000&|@(2&^) x) -: 2 (10000&|@^) x=: 5 10 ?@$ 1e6
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
280
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
281
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
282 NB. x^y on infinite arguments -------------------------------------------
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
283
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
284 0 -: 0 ^ _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
285 _ -: 0 ^ __
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
286 _ -: 2.5 ^ _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
287 1 -: 1 ^ _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
288 0 -: _0.5 ^ _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
289 0 -: 0.5 ^ _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
290
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
291 'domain error' -: _1 ^ etx _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
292 'domain error' -: _1 ^ etx __
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
293 'domain error' -: _5.15 ^ etx _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
294 'domain error' -: _0.15 ^ etx __
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
295
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
296 0 0 0 1 _ _ _ _ -: ({._ 1x)^ __ _5 _4 0 1 2 3x _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
297 0 0 0 1 __ _ __ -: ({.__ 1x)^ __ _5 _4 0 1 2 3x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
298 0 1 _ _ _ -: 0 1 2 3x _ ^ {. _ 1x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
299 0 0 0 _ 1 0 0 0 -: __ _5 _4 0 1 2 3x _ ^ {. __ 1x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
300
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
301 'domain error' -: _5x ^ etx {._ 1x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
302 'domain error' -: __ 1x ^ etx {._ 1x
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
303
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
304 0 0 0 1 _ _ _ _ -: ({._ 1r1)^ __ _5 _4 0 1 2 3r1 _
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
305 0 0 0 1 __ _ __ -: ({.__ 1r1)^ __ _5 _4 0 1 2 3r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
306 0 1 _ _ _ -: 0 1 2 3r1 _ ^ {. _ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
307 0 0 0 _ 1 0 0 0 -: __ _5 _4 0 1 2 3r1 _ ^ {. __ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
308
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
309 0 -: _1r2 ^ {. _ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
310 0 -: 1r2 ^ {. _ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
311 _ -: 1r2 ^ {. __ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
312
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
313 'domain error' -: _5r1 ^ etx {._ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
314 'domain error' -: __ 1r1 ^ etx {._ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
315 'domain error' -: _1r2 ^ etx {.__ 1r1
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
316
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
317
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
318 4!:55 ;:'a exp f h pow rou s x y'
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
319
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
diff changeset
320