1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
|
(kill(all),
fpprec : 25,
load("hypergeometric"),
expand_hypergeometric : false,
max_fpprec : 1000,
domain : complex,
close_p(a,b, dig) := is(abs(a-b) < 10^-dig * abs(a)),
start : absolute_real_time (),
0);
0$
close_p(nfloat(sqrt(3), [], 10), 1.7320508075688772, 10);
true$
string(nfloat(sqrt(3), [], 25));
"1.732050807568877293527446b0"$
string(nfloat(sqrt(3), [], 75));
"1.73205080756887729352744634150587236694280525381038062805580697945193301691b0"$
close_p(nfloat(4/5,[],10), 0.8d0,10);
true$
close_p(nfloat(4/5,[], 42), 8.0b-1, 41);
true$
string(nfloat(%pi,[],50));
"3.1415926535897932384626433832795028841971693993751b0"$
string(nfloat(%e,[],50));
"2.7182818284590452353602874713526624977572470937b0"$
string(nfloat(%phi,[],50));
"1.6180339887498948482045868343656381177203091798058b0"$
close_p(nfloat(sin(6),[], 20), (block([fpprec : 20], bfloat(sin(6)))), 20);
true$
close_p(nfloat(x*cos(x),[x= 6.7], 10), (block([fpprec : 10], bfloat(6.7b0 * cos(6.7b0)))), 10);
true$
every(lambda([a,b], close_p(a,b,5)), nfloat([7.03!, exp(3.08), cos(9.13)], [], 5),
[5354.48837930151,21.75840239619708,-0.95686667438726]);
true$
string(nfloat(a < b, [a = 3, b = 5], 28));
"3.0b0 < 5.0b0"$
is(abs(nfloat(sqrt(2) + sqrt(3) - sqrt(5 + 2 * sqrt(6)),[], 47)) < 1/10^47);
true$
is(abs(nfloat(10^75*( 1-sin(x)^2-cos(x)^2),[x = 12.807], 50)) < 1/10^50);
true$
string(nfloat(exp(cos(1000)),[], 23));
"1.7548424420723242327897b0"$
string(makelist(nfloat(1 - sqrt(1 - exp(-x)),[x = 25 * k], 42),k,1,5));
"[6.94397193250611967043059826014130435829847b-12,9.64374923981958891509255092525427404224488b-23,1.33931848097919003989446756792808961123831b-33,0.0b0,0.0b0]"$
string(nfloat((sqrt(x+1) - 1)/x,[x = 1.0e-10], 42));
"4.99999999987500000000624999544558471101821b-1"$
every(lambda([a,b], close_p(a,b,18)), nfloat([log(x), cosh(x), sinh(x), csc(x), tan(x), asinh(x), atanh(x), acos(x), asin(x)],[x = 2.082 - %i], 19),
[8.371185747911492687b-1-4.477697799919425369b-1*%i,2.200428274786527453b0-3.322048884197778467b0*%i,
2.133063058316276022b0-3.426963992701597692b0*%i,6.283668414956201133b-1-2.68443764645675033b-1*%i,-
1.119096519392880926b0*%i-2.633126796132956158b-1,1.560007670336818004b0-4.142419418455600467b-1*%i,
3.940482433774016247b-1-1.354653476519420094b0*%i,1.502038751933996562b0*%i+4.876895240070315755b-1,
1.083106802787865044b0-1.502038751933996562b0*%i]);
true$
block([fpprec : 58], close_p(nfloat(acosh(x), [x = 28.07b0 + %i * 22/7], 42),
bfloat(acosh(28.07b0 + %i * 22/7)), 40));
true$
nfloat(abs(x),[x = 1.0],10);
1.0$
(declare(x,complex), nfloat(conjugate(x),[x = 1.3 + %i],10));
1.3 - 1.0 * %i$
nfloat(matrix([sqrt(2), sqrt(3)],[sqrt(8), sqrt(10)]),[],10);
matrix([1.4142135623730951, 1.7320508075688772], [2.8284271247461903, 3.1622776601683795])$
nfloat([a = b, a < b, a <= b, a > b, a >- b],[a = 2, b = 4],10);
[2.0 = 4.0, 2.0 < 4.0, 2.0 <= 4.0, 2.0 > 4.0, 2.0 > - 4.0]$
nfloat(map('rhs, solve(x^3 + x + 1,x)),[], 10);
[0.341163901914009 - 1.161541399997251 * %i, 1.161541399997251 * %i + 0.341163901914009, - 0.682327803828018]$
nfloat([acos(x), asin(x), atan(x), acosh(x), asinh(x), atanh(x)], [x=0.12],10);
[1.4505064444001086, 0.12028988239478806, 0.11942892601833845, 1.4505064444001083 *%i,
0.11971385039877737, 0.12058102840844402]$
nfloat('(sum(1/(k + 1) - 1/k,k,1,42)),[], 10);
-0.976744186046511$
string(nfloat(sin(10^50),[],53));
"-7.8967249342931008271028953991740775396008340462140272b-1"$
string(nfloat(sin(10^50),[],100));
"-7.89672493429310082710289539917407753960083404621402719145780873622189996980060989863343675758968847b-1"$
string(nfloat(%pi,[],50));
"3.1415926535897932384626433832795028841971693993751b0"$
string(nfloat(1/7,[],50));
"1.4285714285714285714285714285714285714285714285714b-1"$
string(nfloat(gamma(33/10),[],50));
"2.6834373819557687935963273147667112586281870043548b0"$
nfloat('(abs(f)),[f=28.07], 10);
28.07$
nfloat('(realpart(f)),[f=28.07 - %i], 10);
28.07$
nfloat('(imagpart(f)),[f=28.07 - %i], 10);
-1.0$
nfloat(a * b, [a = 6.7, b= 12],23);
8.0400000000000002131628b1$
close_p(nfloat(hypergeometric([2,2/3],[5/8], x), [x = 1.2], 20),
2.8881276944825158378b1 - 3.7728678507291344336b0 * %i,20);
true$
nfloat(a*b,[a = 1.3, b = 42], 16);
5.46b1$
nfloat(f!, [f= 4], 20);
2.4b1$
close_p(nfloat(atan2(x,y),[x=1.2b0, y = -8.1b0], 25),2.994514298201390738133141b0, fpprec-2);
true$
close_p(nfloat('(realpart(x)),[x= 1.2 + %i]),1.2,20);
true$
close_p(nfloat('(imagpart(x)),[x= 1.2 + %i]),1.0,20);
true$
nfloat('(a,b),[a=3,b=7]);
7.0b0$
nfloat([x,y],[x=1,y=2],23);
[1.0b0,2.0b0]$
block([fpprec : 25], close_p(nfloat('(product(sin(k*x),k,1,28)),[x = 0.1b0], 25),
1.909750489550340725818861b-6, 24));
true$
(time : (absolute_real_time ()- start), print(time), is(time < 4));
true$
nfloat(%i,[], 57);
1.0b0*%i$
nfloat(%i,[], 5);
1.0*%i$
/* nfloat bug - ID: 3358420 */
nfloat(1.0*%i,[],20);
1.0b0 * %i$
/* nfloat - ID: 3455326 */
close_p(nfloat((sqrt(1+x)-1)/x,[x=1/10^15],15),1/2, 15);
true$
(remvalue(start,time), remove(x,complex), reset(domain, fpprec, expand_hypergeometric, max_fpprec),
remfunction(close_p),0);
0$
|
