Solving a system in a matrix

Posted on 2008-11-20 10:10 By KaiTaylor (1)

Categories:

How do I ...? with Maple

Hi,

I'm trying to solve a system of equations in a matrix. The matrix is of the form <(1/5)*sqrt(6)*c[3, 1]+(1/7)*sqrt(6)*c[3, 3],0,(2/7)*c[3, 3]*c[3, 1]+(2/9)*c[3, 3]^2>, which I want to equate to <0,0,1> to simultaneously find numerical values for c[3,1] and c[3,3]. I then need to take these values and imput them into a bigger matrix and repeat; I think I can do that but I'm caught up on this part. I've tried various commands from linalg, LinearAlgebra and Student[LinearAlgebra], with very little success.

Any help would be very much appreciated.

Many thanks in advance.

Kai

nonlinear

Comment on 2008-11-20 10:33 from acer ( Maple Rating 4

1410)

The equations are not all linear in those variables.

> eqnlist := Equate(
>   <(1/5)*sqrt(6)*c[3,1]+(1/7)*sqrt(6)*c[3, 3], 0,
>     (2/7)*c[3, 3]*c[3, 1]+(2/9)*c[3, 3]^2>,
>   <0,0,1>
>                   ):

> fsolve({eqnlist[1],eqnlist[3]},{c[3, 1],c[3, 3]});
                {c[3, 1] = 5.303300866, c[3, 3] = -7.424621212}
 
> _EnvExplicit:=true:

> solve(eqnlist,[c[3, 1],c[3, 3]],AllSolutions);
                  1/2                1/2
              15 2               21 2
[[c[3, 1] = - -------, c[3, 3] = -------],
                 4                  4
 
                   1/2                  1/2
               15 2                 21 2
    [c[3, 1] = -------, c[3, 3] = - -------]]
                  4                    4
 
> evalf(%);
[[c[3, 1] = -5.303300858, c[3, 3] = 7.424621200],
 
    [c[3, 1] = 5.303300858, c[3, 3] = -7.424621200]]
acer

Solving a system in a matrix

nonlinear

Comment viewing options