Related

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-solve-true

keepfloat:true;

float(solve (x=0.1734...

Calculate

keepfloat-linsolve

eq1: 5*x + 1920.7229*...

eq2: 1920.7229*x + 73...

eq3: 738907.9375*x + ...

Calculate

? keepfloat;

Calculate

? keepfloat;

Calculate

? keepfloat;

Calculate

keepfloat-solve-true

keepfloat:true;

float(solve (x=0.1734...

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

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
(%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]);
plot```
Run Example
```kill(all);
(%o0)                                done
(%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]);
plot```

Related Help

Help for Keepfloat