Sponsored links: Algebra eBooks
 

Related

keepfloat

? keepfloat;

Calculate

keepfloat-linsolve

? keepfloat;

eq1:-cx-81+bxi=0;

eq2:cz+28+45-bzi=0;

Calculate

keepfloat-linsolve

? keepfloat;

eq1:-cx-81+bxi=0;

eq2:cz+28+45-bzi=0;

Calculate

keepfloat-linsolve

eq1: 5*x + 1920.7229*...

eq2: 1920.7229*x + 73...

eq3: 738907.9375*x + ...

Calculate

keepfloat-linsolve

eq1: 5*x + 1920.7229*...

eq2: 1920.7229*x + 73...

eq3: 738907.9375*x + ...

Calculate

keepfloat-linsolve

? keepfloat;

eq1:-cx-81+bxi=0;

eq2:cz+28+45-bzi=0;

Calculate

keepfloat-solve

solve(1.0-x,x), keepf...

Calculate

keepfloat-solve-true

keepfloat:true;

float(solve (x=0.1734...

Calculate

keepfloat-solve

solve(1.0-x,x), keepf...

Calculate

keepfloat-linsolve

?keepfloat;

eq1:-cx-81+bxi=0;

eq2:cz+28+45-bzi=0;

Calculate

keepfloat

Run Example
(%i1)rat(x/2.0);

rat: replaced 0.5 by 1/2 = 0.5
                                       x
(%o1)/R/                               -
                                       2
(%i2)  rat(x/2.0), keepfloat;
(%o2)/R/                             0.5 x
(%i3)  solve(1.0-x,x), keepfloat;

rat: replaced 1.0 by 1/1 = 1.0
(%o3)                               [x = 1]
(%i4) 
Run Example
kill(all);
(%o0)                                done
(%i1) load("log10");
(%o1)          /usr/share/maxima/5.21.1/share/contrib/log10.mac
(%i2) keepfloat:true;
(%o2)                                true
(%i3) ratkeepfloat:true;
(%o3)                                true
(%i4) ratprint:false;
(%o4)                                false
(%i5) dB(x):= 20*log10(x);
(%o5)                        dB(x) := 20 log10(x)
(%i6) R1:600;
(%o6)                                 600
(%i7) R2:600;
(%o7)                                 600
(%i8) C1:13.85*10^-9;
(%o8)                              1.385E-8
(%i9) C2:26.82*10^-9;
(%o9)                        2.6820000000000001E-8
(%i10) C3:33.06*10^-9;
(%o10)                       3.3060000000000006E-8
(%i11) C4:9.21*10^-9;
(%o11)                       9.2100000000000011E-9
(%i12) C5:31.47*10^-9;
(%o12)                       3.1470000000000003E-8
(%i13) L1:12.88*10^-3;
(%o13)                              0.01288
(%i14) L2:26.49*10^-3;
(%o14)                              0.02649
(%i15) eq_v1: (v1-x)/R1 + C1*s*(v1-0) - (v2-v1)/(s*L1) = 0;
           v1 - x   77.63975155279502 (v1 - v2)
(%o15)     ------ + --------------------------- + 1.385E-8 s v1 = 0
            600                  s
(%i16) eq_v2: (v2-v1)/(s*L1) + C2*s*(v2-0) - C3*s*(v3-v2) = 0;
                                             77.63975155279502 (v2 - v1)
(%o16) - 3.3060000000000006E-8 s (v3 - v2) + ---------------------------
                                                          s
                                               + 2.6820000000000001E-8 s v2 = 0
(%i17) eq_v3: C3*s*(v3-v2) + C4*s*(v3-0) - (v4-v3)/(s*L2) = 0;
       37.75009437523594 (v3 - v4)
(%o17) --------------------------- + 3.3060000000000006E-8 s (v3 - v2)
                    s
                                               + 9.2100000000000011E-9 s v3 = 0
(%i18) eq_v4: (v4-v3)/(s*L2) + C5*s*(v4-0) + v4/R2 = 0;
       37.75009437523594 (v4 - v3)                                v4
(%o18) --------------------------- + 3.1470000000000003E-8 s v4 + --- = 0
                    s                                             600
(%i19) H(s):= (1/0.12248029984818)*second(first(solve([eq_v1,eq_v2,eq_v3,eq_v4], [v4,v1,v2,v3])/x)[1]);
                      1
(%o19) H(s) := ----------------
               0.12248029984818
                       solve([eq_v1, eq_v2, eq_v3, eq_v4], [v4, v1, v2, v3])
          second(first(-----------------------------------------------------) )
                                                 x
                                                                             1
(%i20) float(H(s));
                                                           6
(%o20) (5.6790255971822221E+42 s)/(2.6998598124854999E+18 s
                           5                           4
 + 4.6787878036331813E+23 s  + 4.5984440412858664E+28 s
                           3                           2
 + 2.9943547699319576E+33 s  + 1.2114544840586922E+38 s
 + 3.1027079484570518E+42 s + 3.5065979068283219E+46)
(%i21) dividethru(H(s),35065979068283216162754605785751246244066600000);
(%o21) dividethru((5.6790255971822221E+42 s)
                       6                             5
/(2699859812485500000 s  + 467878780363318124692500 s
                                  4                                       3
 + 45984440412858665627860199375 s  + 2994354769931957449785561558212640 s
                                            2
 + 121145448405869224811779225339451316138 s
 + 3102707948457051766887788569335640443000000 s
 + 35065979068283216162754605785751246244066600000), 
35065979068283216162754605785751246244066600000)
(%i22) plot2d(float(dB(cabs(H(2*%pi*%i*f)))),[f,10,20000],[logx,true]);
plotplot2d(float(dB(cabs(H(2*%pi*%i*f)))),[f,10,20000],[logx,true]);
Run Example
kill(all);
(%o0)                                done
(%i1) load("log10");
(%o1)          /usr/share/maxima/5.21.1/share/contrib/log10.mac
(%i2) keepfloat:true;
(%o2)                                true
(%i3) ev(keepfloat,true);
(%o3)                                true
(%i4) numer:true;
(%o4)                                true
(%i5) dB(x):= 20*log10(x);
(%o5)                        dB(x) := 20 log10(x)
(%i6) R1:600;
(%o6)                                 600
(%i7) R2:600;
(%o7)                                 600
(%i8) C1:13.85*10^-9;
(%o8)                              1.385E-8
(%i9) C2:26.82*10^-9;
(%o9)                        2.6820000000000001E-8
(%i10) C3:33.06*10^-9;
(%o10)                       3.3060000000000006E-8
(%i11) C4:9.21*10^-9;
(%o11)                       9.2100000000000011E-9
(%i12) C5:31.47*10^-9;
(%o12)                       3.1470000000000003E-8
(%i13) L1:12.88*10^-3;
(%o13)                              0.01288
(%i14) L2:26.49*10^-3;
(%o14)                              0.02649
(%i15) eq_v1: (v1-x)/R1 + C1*s*(v1-0) - (v2-v1)/(s*L1) = 0;
                                     77.63975155279502 (v1 - v2)
(%o15) 0.0016666666666667 (v1 - x) + ---------------------------
                                                  s
                                                            + 1.385E-8 s v1 = 0
(%i16) eq_v2: (v2-v1)/(s*L1) + C2*s*(v2-0) - C3*s*(v3-v2) = 0;
                                             77.63975155279502 (v2 - v1)
(%o16) - 3.3060000000000006E-8 s (v3 - v2) + ---------------------------
                                                          s
                                               + 2.6820000000000001E-8 s v2 = 0
(%i17) eq_v3: C3*s*(v3-v2) + C4*s*(v3-0) - (v4-v3)/(s*L2) = 0;
       37.75009437523594 (v3 - v4)
(%o17) --------------------------- + 3.3060000000000006E-8 s (v3 - v2)
                    s
                                               + 9.2100000000000011E-9 s v3 = 0
(%i18) eq_v4: (v4-v3)/(s*L2) + C5*s*(v4-0) + v4/R2 = 0;
       37.75009437523594 (v4 - v3)
(%o18) --------------------------- + 3.1470000000000003E-8 s v4
                    s
                                                    + 0.0016666666666667 v4 = 0
(%i19) H(s):= (1/0.12248029984818)*second(first(solve([eq_v1,eq_v2,eq_v3,eq_v4], [v4,v1,v2,v3])/x)[1]);
                      1
(%o19) H(s) := ----------------
               0.12248029984818
                       solve([eq_v1, eq_v2, eq_v3, eq_v4], [v4, v1, v2, v3])
          second(first(-----------------------------------------------------) )
                                                 x
                                                                             1
(%i20) float(H(s));

rat: replaced 1.385E-8 by 1/72202166 = 1.3850000012465E-8

rat: replaced 77.63975155279502 by 12500/161 = 77.63975155279503

rat: replaced 0.0016666666666667 by 1/600 = 0.0016666666666667

rat: replaced 2.6820000000000001E-8 by 1/37285608 = 2.6819999824060803E-8

rat: replaced 77.63975155279502 by 12500/161 = 77.63975155279503

rat: replaced -3.306000000000001E-8 by -1/30248034 = -3.30599998664376E-8

rat: replaced 9.2100000000000011E-9 by 1/108577633 = 9.2100000006447007E-9

rat: replaced 3.3060000000000006E-8 by 1/30248034 = 3.3059999866437601E-8

rat: replaced 37.75009437523594 by 99849/2645 = 37.75009451795841

rat: replaced 0.0016666666666667 by 1/600 = 0.0016666666666667

rat: replaced 3.1470000000000003E-8 by 1/31776295 = 3.1469999885134501E-8

rat: replaced 37.75009437523594 by 99849/2645 = 37.75009451795841
                                                           6
(%o20) (5.6790255971822221E+42 s)/(2.6998598124854999E+18 s
                           5                           4
 + 4.6787878036331813E+23 s  + 4.5984440412858664E+28 s
                           3                           2
 + 2.9943547699319576E+33 s  + 1.2114544840586922E+38 s
 + 3.1027079484570518E+42 s + 3.5065979068283219E+46)
(%i21) H(s);

rat: replaced 1.385E-8 by 1/72202166 = 1.3850000012465E-8

rat: replaced 77.63975155279502 by 12500/161 = 77.63975155279503

rat: replaced 0.0016666666666667 by 1/600 = 0.0016666666666667

rat: replaced 2.6820000000000001E-8 by 1/37285608 = 2.6819999824060803E-8

rat: replaced 77.63975155279502 by 12500/161 = 77.63975155279503

rat: replaced -3.306000000000001E-8 by -1/30248034 = -3.30599998664376E-8

rat: replaced 9.2100000000000011E-9 by 1/108577633 = 9.2100000006447007E-9

rat: replaced 3.3060000000000006E-8 by 1/30248034 = 3.3059999866437601E-8

rat: replaced 37.75009437523594 by 99849/2645 = 37.75009451795841

rat: replaced 0.0016666666666667 by 1/600 = 0.0016666666666667

rat: replaced 3.1470000000000003E-8 by 1/31776295 = 3.1469999885134501E-8

rat: replaced 37.75009437523594 by 99849/2645 = 37.75009451795841
                                                        6
(%o21) (5.6790255971822221E+42 s)/(2699859812485500000 s
                             5                                  4
 + 467878780363318124692500 s  + 45984440412858665627860199375 s
                                       3
 + 2994354769931957449785561558212640 s
                                            2
 + 121145448405869224811779225339451316138 s
 + 3102707948457051766887788569335640443000000 s
 + 35065979068283216162754605785751246244066600000)
(%i22) plot2d(float(dB(cabs(H(2*%pi*%i*f)))),[f,10,20000],[logx,true]);
plotplot2d(float(dB(cabs(H(2*%pi*%i*f)))),[f,10,20000],[logx,true]);

Related Help

Help for Keepfloat