Related

globalsolve-linsolve-true

globalsolve: true;

eq1:(phi-phi3)*h2+(ph...

eq2:(phi-phi3)*h2*b+(...

Calculate

globalsolve-linsolve

e1 : (0^3)*k + (0^2)*...

e2 : 3*(0^2)*k + 2*0*...

e3 : (10^3)*k + (10^2...

Calculate

globalsolve-linsolve-true

globalsolve: true ;

eq1:phi1 + phi2 + phi...

eq2:J*((phi3*h/e3) - ...

Calculate

globalsolve-linsolve-plot2d-true

globalsolve : true;

w10(x):=(q/6*x^3+c12/...

w11(x):=(q/2*x^2+c12*...

Calculate

globalsolve-linsolve-phi-true

globalsolve: true;

l:325;

h1:150;

Calculate

l1: 303.4;

l2: 601.684;

l3: 45;

Calculate

globalsolve-linsolve-true

globalsolve: true;

eq1:R1+(phi1-phi3)*32...

eq2:R2+(phi2-phi3)*32...

Calculate

globalsolve-linsolve-programmode

eq1:10*a+b=3*(a+b)+3;

eq2:a*a-b*b=2*(a-b)^2;

[globalsolve:false, p...

Calculate

globalsolve-linsolve-true

globalsolve: true;

eq1: phi1*h1-phi2*h2+...

eq2: (phi1-phi2)*h2+(...

Calculate

globalsolve-linsolve-true

globalsolve: true;

linsolve ([x + 3*y =...

x;

Calculate

globalsolve

Run Example
```(%i1)e1: x1*(sqrt(-d^2-c^2-b^2)) + x2*b + x3*c + x4*d = 0;
2    2    2
(%o1)          d x4 + c x3 + b x2 + sqrt(- d  - c  - b ) x1 = 0
(%i2) e2: x1*(-b) + x2*sqrt(-d^2-c^2-b^2) + x3*(-d) + x4*c = 1;
2    2    2
(%o2)          c x4 - d x3 + sqrt(- d  - c  - b ) x2 - b x1 = 1
(%i3) e3: x1*(-c) + x2*d + x3*sqrt(-d^2-c^2-b^2) + x4*(-b) = (-c*sqrt(-d^2-c^2-b^2) - bd)/(d^2 + c^2);
2    2    2
(%o3) - b x4 + sqrt(- d  - c  - b ) x3 + d x2 - c x1 =
2    2    2
- c sqrt(- d  - c  - b ) - bd
-----------------------------
2    2
d  + c
(%i4) e4: x1*(-d) + x2*(-c) + x3*b + x4*sqrt(-d^2-c^2-b^2) = (-d*sqrt(-d^2-c^2-b^2) + bc)/(d^2 + c^2);
2    2    2
(%o4) sqrt(- d  - c  - b ) x4 + b x3 - c x2 - d x1 =
2    2    2
bc - d sqrt(- d  - c  - b )
---------------------------
2    2
d  + c
(%i5)  [globalsolve: true, programmode: true, solvradcan: true, llinsolce_param: true ];
(%o5)                      [true, true, true, true]
(%i6) linsolv([e1, e2, e3, e4], [x1, x2, x3, x4]);
2    2    2
(%o6) linsolv([d x4 + c x3 + b x2 + sqrt(- d  - c  - b ) x1 = 0,
2    2    2
c x4 - d x3 + sqrt(- d  - c  - b ) x2 - b x1 = 1,
2    2    2
- b x4 + sqrt(- d  - c  - b ) x3 + d x2 - c x1 =
2    2    2
- c sqrt(- d  - c  - b ) - bd          2    2    2
-----------------------------, sqrt(- d  - c  - b ) x4 + b x3 - c x2 - d x1 =
2    2
d  + c
2    2    2
bc - d sqrt(- d  - c  - b )
---------------------------], [x1, x2, x3, x4])
2    2
d  + c
(%i7) ```
Run Example
```eq1:10*a+b=3*(a+b)+3;
(%o1)                      b + 10 a = 3 (b + a) + 3
(%i2) eq2:a*a-b*b=2*(a-b)^2;
2    2            2
(%o2)                        a  - b  = 2 (a - b)
(%i3) [globalsolve:false, programmode:true];
(%o3)                            [false, true]
(%i4) linsolve([eq1,eq2],[a,b]);
3 a              3 b
(%o4)                  [a = ---------, b = - ---------]
2 b + 7 a        2 b + 7 a
(%i5) ```
Run Example
```globalsolve: truel: 325;
(%o1)                                 325
(%i2) b: 650;
(%o2)                                 650
(%i3) h1: 150;
(%o3)                                 150
(%i4) h2: 225;
(%o4)                                 225
(%i5) e1: 1.6;
(%o5)                                 1.6
(%i6) e2: 2.5;
(%o6)                                 2.5
(%i7) e3: 1.6;
(%o7)                                 1.6
(%i8) e4: 1.6;
(%o8)                                 1.6
(%i9) E: 70000;
(%o9)                                70000
(%i10) Mx: 24000000;
(%o10)                             24000000
(%i11) S1: 400;
(%o11)                                400
(%i12) S2: 1200;
(%o12)                               1200
(%i13) l1: (b^2 + ((h2 - h1)/2)^2)^(1/2);
25 sqrt(2713)
(%o13)                           -------------
2
(%i14) A: b*h1+b*(h2-h1/2);
(%o14)                              195000
(%i15) G: E/2.6;
(%o15)                         26923.07692307692
(%i16) phi: Mx/(2*A);
800
(%o16)                                ---
13
(%i17) eq1: (phi1-phi)*h1+(phi-phi2)*h2=0;
800                       800
(%o17)              225 (--- - phi2) + 150 (phi1 - ---) = 0
13                        13
(%i18) eq2: (phi-phi2)*h2+(phi-phi3)*h1=0;
800                800
(%o18)              150 (--- - phi3) + 225 (--- - phi2) = 0
13                 13
(%i19) eq3: (R1+l*(phi3-phi1))=0;
(%o19)                     R1 + l (phi3 - phi1) = 0
(%i20) eq4: (-R2+l*(phi3-phi2))=0;
(%o20)                     l (phi3 - phi2) - R2 = 0
(%i21) [phi2,phi3]: linsolve([eq1,eq2],[phi2,phi3]);
26 phi1 + 800           13 phi1 - 1600
(%o21)          [phi2 : -------------, phi3 : - --------------]
39                       13
(%i22) phi2;
26 phi1 + 800
(%o22)                       phi2 : -------------
39
(%i23) phi3;
13 phi1 - 1600
(%o23)                      phi3 : - --------------
13
(%i24) Wl1: l^3*(phi3-phi1)^2/ ( 6* E * S1);
3            13 phi1 - 1600         2
l  ((phi3 : - --------------) - phi1)
13
(%o24)              --------------------------------------
168000000
(%i25) Wl2: l^3*(phi2-phi3)^2/ ( 6* E * S2);
3          26 phi1 + 800             13 phi1 - 1600 2
l  ((phi2 : -------------) - phi3 : - --------------)
39                         13
(%o25)      ------------------------------------------------------
504000000
(%i26) Wp1: phi1^2*h1*l/( 2*G*e1);
2
(%o26)                    0.0017410714285714 l phi1
(%i27) Wp2: phi2^2*h2*l/( 2*G*e2);
26 phi1 + 800 2
(%o27)           0.0016714285714286 l (phi2 : -------------)
39
(%i28) Wp3: phi3^2*l1*l/( 2*G*e3);
13 phi1 - 1600 2
(%o28)   1.4508928571428572E-4 sqrt(2713) l (phi3 : - --------------)
13
(%i29) Wtot: 2*Wl1+2*Wl2+2*Wp3+Wp2+Wp1;
3            13 phi1 - 1600         2
l  ((phi3 : - --------------) - phi1)
13
(%o29) --------------------------------------
84000000
13 phi1 - 1600 2
+ 2.9017857142857144E-4 sqrt(2713) l (phi3 : - --------------)
13
3          26 phi1 + 800             13 phi1 - 1600 2
l  ((phi2 : -------------) - phi3 : - --------------)
39                         13
+ ------------------------------------------------------
252000000
26 phi1 + 800 2                            2
+ 0.0016714285714286 l (phi2 : -------------)  + 0.0017410714285714 l phi1
39
(%i30) phi1 : solve(diff(Wtot,phi1)=0,phi1);

rat: replaced 0.0034821428571429 by 39/11200 = 0.0034821428571429

rat: replaced 0.0033428571428571 by 117/35000 = 0.0033428571428571

rat: replaced 5.8035714285714288E-4 by 13/22400 = 5.8035714285714288E-4
13 phi1 - 1600
(%o30) [phi1 = ((phi3 : - --------------)
13
2       d             13 phi1 - 1600       d           26 phi1 + 800
(l  (4 (----- (phi3 : - --------------)) - ----- (phi2 : -------------) - 3)
dphi1                 13           dphi1              39
d             13 phi1 - 1600
+ 73125 sqrt(2713) (----- (phi3 : - --------------)))
dphi1                 13
26 phi1 + 800    2    d           26 phi1 + 800
+ (phi2 : -------------) (l  (----- (phi2 : -------------)
39             dphi1              39
d             13 phi1 - 1600               d           26 phi1 + 800
- ----- (phi3 : - --------------)) + 421200 (----- (phi2 : -------------))))
dphi1                 13                   dphi1              39
2       d             13 phi1 - 1600
/(l  (3 (----- (phi3 : - --------------)) - 3) - 438750)]
dphi1                 13
(%i31) float(phi1);
(%o31) [phi1 = ((phi3 : - 0.076923076923077 (13.0 phi1 - 1600.0))
1.0
2         d
(l  (4.0 (-------- (phi3 : - 0.076923076923077 (13.0 phi1 - 1600.0)))
1.0
dphi1
1.0
d
- 1.0 (-------- (phi2 : 0.025641025641026 (26.0 phi1 + 800.0))) - 3.0)
1.0
dphi1
1.0
d
+ 3808822.868108335 (-------- (phi3 : - 0.076923076923077
1.0
dphi1
(13.0 phi1 - 1600.0)))) + (phi2 : 0.025641025641026 (26.0 phi1 + 800.0))
1.0
2    d
(l  (-------- (phi2 : 0.025641025641026 (26.0 phi1 + 800.0))
1.0
dphi1
1.0
d
- 1.0 (-------- (phi3 : - 0.076923076923077 (13.0 phi1 - 1600.0))))
1.0
dphi1
1.0
d
+ 421200.0 (-------- (phi2 : 0.025641025641026 (26.0 phi1 + 800.0)))))
1.0
dphi1
1.0
2         d
/(l  (3.0 (-------- (phi3 : - 0.076923076923077 (13.0 phi1 - 1600.0))) - 3.0)
1.0
dphi1
- 438750.0)]
(%i32) ```

Related Help

Help for Globalsolve