[交流] 【求助】帮忙求解四元方程组 已有4人参与

Originally posted by zxzj05 at 2010-10-09 15:39:31:
试试maple

(%i1) eq1: a*b/d=0.005754;
                                a b
(%o1)                           --- = 0.005754
                                 d
(%i2) eq2: a*c/(1-d)=4.169*10^(-6);
                          a c
(%o2)                    ----- = 4.1689999999999994E-6
                         1 - d
(%i3) eq3: (0.2-a)/(0.25-b)=4/d;
                                 0.2 - a    4
(%o3)                            -------- = -
                                 0.25 - b   d
(%i4) eq4: (0.2-a)/(0.03-c)=24/(7*d);
                                0.2 - a    24
(%o4)                           -------- = ---
                                0.03 - c   7 d
(%i5) solve([eq1,eq2,eq3,eq4],[a,b,c,d]);

rat: replaced -0.005754 by -433/75252 = -0.005753999893691

rat: replaced -4.168999999999999E-6 by -37/8875030 = -4.168999992112703E-6

rat: replaced 0.2 by 1/5 = 0.2

rat: replaced 0.25 by 1/4 = 0.25

rat: replaced 0.2 by 1/5 = 0.2

rat: replaced 0.03 by 3/100 = 0.03
(%o5) [[a = 0.18571751923947, b = 0.22416573626101,
c = - 1.3996884948583785E-4, d = 7.235228539576366],
[a = 0.01412678085481, b = 0.22439931350114, c = 1.3252708943556937E-4,
d = 0.55092802172929], [a = 1.3977106364338493E-4, b = 0.24969694088565,
c = 0.029646433041302, d = 0.0060654157927555]]

Originally posted by snoopyzhao at 2010-10-09 18:14:33:
Maxima 的结果。

[code]
(%i1) eq1: a*b/d=0.005754;
                                a b
(%o1)                           --- = 0.005754
                                 d
(%i2) eq2: a*c/(1-d)= ...