Sponsored links: Algebra eBooks
 

Related

globalsolve-integrate-linsolve-true

globalsolve: true ;

equ1: -R1+(phi1-phi3)...

equ2: -R2+(phi3-phi2)...

Calculate

globalsolve-linsolve-sqrt-true

h: 45;

e3: 2;

l1: 2*sqrt(150*150+(4...

Calculate

globalsolve-linsolve

globalsolve: trueeq1:...

eq2:phi*(h2-h1)+phi1*...

linsolve([eq1,eq2],[p...

Calculate

globalsolve-linsolve-true

globalsolve: true;

eq1:h2(phi-phi2)+h1*(...

eq2:-phi2*h2-phi*h2+(...

Calculate

globalsolve-integrate-linsolve-true

globalsolve: true ;

equ1: -R1+(phi1-phi3)...

equ2: -R2+(phi3-phi2)...

Calculate

globalsolve-linsolve-phi-true

globalsolve: true;

l:325;

h1:150;

Calculate

globalsolve-linsolve-true

e1: 1.6;

e2: 2.5;

e3: 1.6;

Calculate

globalsolve-linsolve-programmode

e1: -I1-I2+I3=0;

e2: 69.1*I2+33.0*I3=...

e3: -115.4*I1-33.0*I...

Calculate

globalsolve-true

eq1:-R1+(Phi1-Phi3)*l=0;

eq2:R2+(Phi3-Phi2)*l=0;

eq3:(h1+h2)*Phi-Phi3*...

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
globalsolve:true;
(%o1)                                true
(%i2) l:325;
(%o2)                                 325
(%i3) b:650;
(%o3)                                 650
(%i4) h1:150;
(%o4)                                 150
(%i5) h2:225;
(%o5)                                 225
(%i6) e1:1.6;
(%o6)                                 1.6
(%i7) e2:2.5;
(%o7)                                 2.5
(%i8) e3:1.6;
(%o8)                                 1.6
(%i9) e4:1.6;
(%o9)                                 1.6
(%i10) nu:0.3;
(%o10)                                0.3
(%i11) E:70000;
(%o11)                               70000
(%i12) S1:400;
(%o12)                                400
(%i13) S2:1200;
(%o13)                               1200
(%i14) phi=24000000/(2*b*(h2+h1)/2);
                                        1280
(%o14)                            phi = ----
                                         13
(%i15) G:E/(2*(1+nu));
(%o15)                         26923.07692307692
(%i16) phi1:(1/h1)*(-phi*(h2-h1)+phi2*h2);
                               225 phi2 - 75 phi
(%o16)                         -----------------
                                      150
(%i17) phi3:(1/h1)*(phi*(h2+h1)-phi2*h2);
                              375 phi - 225 phi2
(%o17)                        ------------------
                                     150
(%i18) Ws1:(l*l*l)/(h1*h1*6*E*S1)*((((phi(h2+h1)-phi2*h2)^2)+((-phi*(h2-h1)+phi2*h2)^2)-2*((phi(h2+h1)-phi2*h2)*(-phi*(h2-h1)+phi2*h2))));
                                 2
(%o18) (2197 ((225 phi2 - 75 phi)  - 2 (phi(375) - 225 phi2)
                                                                  2
                       (225 phi2 - 75 phi) + (phi(375) - 225 phi2) ))/241920000
(%i19) Ws2:(l*l*l)/(h1*h1*6*E*S2)*(phi2-phi3)^2;
                                    375 phi - 225 phi2 2
                       2197 (phi2 - ------------------)
                                           150
(%o19)                 ---------------------------------
                                   725760000
(%i20) Wp1:phi1*phi1*l*h1/(G*e1*2);
                                                           2
(%o20)            2.5148809523809522E-5 (225 phi2 - 75 phi)
(%i21) Wp2:0;
(%o21)                                 0
(%i22) Wp3:0;
(%o22)                                 0
(%i23) W:2*Ws1+2*Ws2+2*Wp3+Wp1+Wp2;
                                 2
(%o23) (2197 ((225 phi2 - 75 phi)  - 2 (phi(375) - 225 phi2)
                                            2
 (225 phi2 - 75 phi) + (phi(375) - 225 phi2) ))/120960000
                                            2
 + 2.5148809523809522E-5 (225 phi2 - 75 phi)
                375 phi - 225 phi2 2
   2197 (phi2 - ------------------)
                       150
 + ---------------------------------
               362880000
(%i24) solve(diff(W,phi2)=0,phi2);

rat: replaced 0.011316964285714 by 507/44800 = 0.011316964285714
                              356413 phi + 2808 phi(375)
(%o24)                [phi2 = --------------------------]
                                       1701013
(%i25) 
Run Example
globalsolve:true;
(%o1)                                true
(%i2) k1(Temp):=10^(7.176-3.85/Temp);
                                      7.176 - 3.85/Temp
(%o2)                   k1(Temp) := 10
(%i3) k2(Temp):=10^(9.942-4.05/Temp);
                                      9.942 - 4.05/Temp
(%o3)                   k2(Temp) := 10
(%i4) OH(pH):=10^(pH-OH);
                                          pH - OH
(%o4)                         OH(pH) := 10
(%i5) ECH(t);
(%o5)                               ECH(t)
(%i6) diff(ECH,t):=k1*ECH+k2*ECH*OH;
(%o6)                 diff(ECH, t) := k1 ECH + k2 ECH OH
(%i7) solve(ECH);
(%o7)                              [ECH = 0]
(%i8) 

Related Help

Help for Globalsolve