Sponsored links: Algebra eBooks
 

Related

ratprint

ratprint:off;

Calculate

ratprint

ratprint:off;

Calculate

ratprint

Run Example
(%i1)rat(0.1234);

rat: replaced 0.1234 by 617/5000 = 0.1234
                                     617
(%o1)/R/                             ----
                                     5000
(%i2) rat(0.12345);

rat: replaced 0.12345 by 2469/20000 = 0.12345
                                     2469
(%o2)/R/                             -----
                                     20000
(%i3) rat(0.12345);

rat: replaced 0.12345 by 2469/20000 = 0.12345
                                     2469
(%o3)/R/                             -----
                                     20000
(%i4) mod(12345, 1000)/100000, numer;
(%o4)                               0.00345
(%i5) lastdigits(fnum, lastn):=block([m:10*fnum, ans, n:lastn], /* it doesn't help when lastn is bigger than the total digits in fnum it returns the whole digits in fnum anyway without erroring! and fnum=%pi stuff works only upto 8 digits cause num(rat(fnum)) has somuch only! */ratprint:false, k:rat(fnum), ans:mod(m, (10^-n))/(mod(m, 10)*denom(k)), return(float(ans)));
(%o5) lastdigits(fnum, lastn) := block([m : 10 fnum, ans, n : lastn], 
                                                   - n
                                          mod(m, 10   )
ratprint : false, k : rat(fnum), ans : -------------------, return(float(ans)))
                                       mod(m, 10) denom(k)
(%i6) lastdigits(0.12345, 3);
(%o6)                        2.0251113811259619E-8
(%i7) 
Run Example
load(basic);
define: warning: redefining the built-in function prog1
define: warning: redefining the built-in function symbolcheck
define: warning: redefining the built-in function push
define: warning: redefining the built-in function pop
define: warning: redefining the built-in function tr_ev
(%o1)           /usr/share/maxima/5.21.1/share/macro/basic.mac
(%i2) load(functs);
(%o2)      /usr/share/maxima/5.21.1/share/simplification/functs.mac
(%i3)  ratprint:false;
(%o3)                                false
(%i4) setify(makelist(1/k, k, 1, ceiling(rat(0.1))));
(%o4)                                 {1}
(%i5) 
Run Example
kill(lastdigits,pidigits, a, b);
(%o1)                                done
(%i2) ratprint:false;
(%o2)                                false
(%i3) pidigits(prec):=block([fpprec:prec+1], ratprint:false, bftorat:false, ratepsilon:2.0e-16, float2bf:false, ans:float((%pi-mod(%pi, 10^(-prec)))), return(ans));
(%o3) pidigits(prec) := block([fpprec : prec + 1], ratprint : false, 
bftorat : false, ratepsilon : 2.0E-16, float2bf : false, 
                             - prec
ans : float(%pi - mod(%pi, 10      )), return(ans))
(%i4) lastdigits(fnum, lastn):=block([m:rat(fnum), n:lastn, ans], /* it doesn't help when lastn is bigger than the total digits in fnum it returns the whole digits in fnum anyway without erroring! and fnum=%pi stuff works only upto 8 digits cause num(rat(fnum)) has somuch only! */ratprint:false, ratio(num):=block([k, a], ratprint:false, for k:1 thru fpprec do (a:10^k/ratdenom(num), if(integerp(a)) then return(a))), ratepsilon:2.0e-16, bftorat:false,  ans:mod(ratio(fnum)*num(m), (10^n))/(ratio(fnum)*denom(m)), return(float(ans)));
(%o4) lastdigits(fnum, lastn) := block([m : rat(fnum), n : lastn, ans], 
ratprint : false, ratio(num) := block([k, a], ratprint : false, 
                                 k
                               10
for k thru fpprec do (a : -------------, if integerp(a) then return(a))), 
                          ratdenom(num)
                                                                       n
                                             mod(ratio(fnum) num(m), 10 )
ratepsilon : 2.0E-16, bftorat : false, ans : ----------------------------, 
                                                 ratio(fnum) denom(m)
return(float(ans)))
(%i5) lastdigits(4.1459, 5);
(%o5)                          4.14590001103631
(%i6) lastdigits(3.1459, 3);
(%o6)                               0.0459
(%i7) b:pidigits(5);
(%o7)                               3.14159
(%i8) lastdigits(3.14159, 3);
(%o8)                               0.00159
(%i9) lastdigits(pidigits(5), 3);
(%o9)                               0.00159
(%i10) pidigits(6);
(%o10)                             3.141592
(%i11) lastdigits(pidigits(6), 3);
(%o11)                       5.9199999999999997E-4
(%i12) lastdigits(4.14592, 3);
(%o12)                              0.00592
(%i13) lastdigits(b, 3);
(%o13)                              0.00159
(%i14) a:[0, float(%pi)];
(%o14)                      [0, 3.141592653589793]
(%i15) a;
(%o15)                      [0, 3.141592653589793]
(%i16) delete(pidigits(16), a);
(%o16)                                [0]
(%i17) delete(pidigits(fpprec), a);
(%o17)                                [0]
(%i18) delete(pidigits(fpprec), [0, float(%pi)]);
(%o18)                                [0]
(%i19) intervals(items, [granul]):=block([ start:first(items), stop:last(items), ans:[], granul:(if not(emptyp(granul)) then first(granul) else 1) ], ratprint:false, if (emptyp(items)) then [] elseif(listp(items)) then    if(emptyp(rest(items))) then []     else     while(float(rat(start))<
float(rat(stop))) do    (ans:delete([],             delete([float(rat(stop)), float(rat(stop+granul))],             if(float(rat(start+granul))>
float(rat(stop)))            then return() else            append(ans, [cons(float(rat(start)),  [float(rat(start+granul))])]))),      start:float(rat(start+granul))),     return(ans));
(%o19) intervals(items, [granul]) := block([start : first(items), 
stop : last(items), ans : [], granul : 
if not emptyp(granul) then first(granul) else 1], ratprint : false, 
if emptyp(items) then [] elseif listp(items)
 then (if emptyp(rest(items)) then [] else (while float(rat(start)) < 
float(rat(stop)) do (ans : delete([], delete([float(rat(stop)), 
float(rat(stop + granul))], if float(rat(start + granul)) > float(rat(stop))
 then return() else append(ans, [cons(float(rat(start)), 
[float(rat(start + granul))])]))), start : float(rat(start + granul))))), 
return(ans))
(%i20) intervals([0, 2*%pi], %pi);
(%o20) [[0.0, 3.141592653589793], [3.141592653589793, 6.283185307179585]]
(%i21) delete(pidigits(fpprec), intervals([0, 2*%pi], %pi));
(%o21) [[0.0, 3.141592653589793], [3.141592653589793, 6.283185307179585]]
(%i22) 
[append,block,cons,delete,do,emptyp,first,float,last,listp,ratprint,rest,return] [assume,diff,ev,false,float,fpprintprec,ic1,inf,kill,lambda,limit,map,numer,ode2,ratprint,ratsimp,rhs,solve,subst] [at,block,denom,float,mod,numer,rat,ratprint] [bfloat,bftorat,block,delete,ev,false,float2bf,if,integerp,mod,numer,ratdenom,ratepsilon,ratprint] [bfloat,bftorat,block,delete,false,if,integerp,mod,ratdenom,ratepsilon,ratprint] [bfloat,bftorat,block,emptyp,equal,first,fpprec,if,is,mod,not,numer,ratepsilon,ratprint,return,round] [bfloat,bftorat,block,emptyp,first,fpprec,if,mod,not,numer,ratepsilon,ratprint,return] [bfloat,bftorat,block,ev,float2bf,makelist,mod,ratepsilon,ratprint,return] [bfloat,bftorat,block,false,if,integerp,mod,ratdenom,ratepsilon,ratprint] [bfloat,bftorat,block,float2bf,makelist,mod,ratepsilon,ratprint,return] [bfloat,block,denom,false,float,mod,ratexpand,ratprint] [bfloat,block,denom,false,float,mod,ratprint,ratsimp] [bftorat,block,cons,delete,do,emptyp,first,float2bf,fpprec,lambda,last,listp,map,maplist,mod,ratepsilon,ratprint,rest,return,second] [bftorat,block,cons,delete,do,emptyp,first,float2bf,fpprec,lambda,last,listp,map,maplist,mod,ratepsilon,ratprint,rest,return] [bftorat,block,cons,delete,do,emptyp,first,float2bf,fpprec,lambda,last,listp,map,mod,ratepsilon,ratprint,rest,return] [bftorat,block,cons,delete,emptyp,false,first,float2bf,fpprec,if,integerp,kill,lambda,listp,map,mod,ratdenom,ratepsilon,ratprint,rest,return] [bftorat,block,delete,ev,false,float,float2bf,if,integerp,mod,ratdenom,ratepsilon,ratprint] [bftorat,block,delete,ev,false,float,if,integerp,mod,ratdenom,ratepsilon,ratprint] [bftorat,block,delete,float,if,integerp,mod,rat,ratdenom,ratepsilon,ratprint] [bftorat,block,emptyp,equal,first,float,fpprec,if,is,mod,not,numer,ratepsilon,ratprint,return,round] [bftorat,block,emptyp,first,float,fpprec,if,integerp,mod,not,numer,ratepsilon,ratprint,return] [bftorat,block,ev,false,float2bf,fpprec,if,integerp,kill,mod,ratdenom,ratepsilon,ratprint,return] [bftorat,block,false,float2bf,fpprec,if,integerp,kill,mod,ratdenom,ratepsilon,ratprint,return] [bftorat,block,float,fpprec,if,integerp,kill,makelist,mod,rat,ratdenom,ratepsilon,ratprint] [bftorat,block,float,if,integerp,mod,num,ratdenom,ratepsilon,ratprint] [bftorat,block,float2bf,fpprec,if,integerp,mod,ratdenom,ratepsilon,ratprint,return] [block,cot,equal,errormsg,false,if,is,limit,plot2d,ratprint,solvetrigwarn] [block,denom,false,float,mod,ratprint] [block,false,float,mod,ratprint] [coeff,false,float,notequal,part,ratprint,round,solve] [concat,denom,false,floor,length,makelist,map,mod,num,ratprint,transpose] [cos,false,numer,ratprint,solve] [define,diff,expand,false,integrate,kill,logabs,ratprint,reset,true] [define,diff,false,integrate,kill,logabs,ratprint,reset,true] [define,false,integrate,kill,logabs,ratprint,reset,true] [diff,exp,false,ic1,ode2,plot2d,ratprint,ratsimp,rhs,subst] [diff,exp,false,ic1,ode2,plot2d,ratprint,ratsimp,rhs] [diff,false,ic1,ode2,plot2d,ratprint,ratsimp,rhs] [diff,false,ic1,ode2,plot2d,ratprint,remvalue,rhs] [diff,false,kill,load,logabs,ratprint,reset,true] [false,float,kill,lhs,plot2d,ratprint,rhs,solve] [false,float,kill,max,min,plot2d,ratprint,rhs,solve] [false,float,kill,plot2d,ratprint,solve] [false,integrate,kill,load,quad_qags,ratprint,sin] [false,integrate,numer,plot2d,ratprint,rhs,solve] [false,kill,numer,ratprint,realroots,solve] [false,kill,numer,ratprint,solve] [false,kill,quad_qags,ratprint,sin] [false,numer,ratprint,solve] [false,quad_qags,ratprint,sin]

Related Help

Help for Ratprint