       Trying values other than random integers

• To: mathgroup at smc.vnet.net
• Subject: [mg85629] Trying values other than random integers
• From: shama shahbaz <shammashahbaz at yahoo.com>
• Date: Sat, 16 Feb 2008 03:28:28 -0500 (EST)

```Dear All
MY problem is actually

i have a matrix
a1*x1 +a2*x2 +a3*x3  +a4*x4   +a5*x5  +a6*x6=0
b1*x1 +b2*x2 +b3*x3  +b4*x4   +b5*x5  +b6*x6=0
c1*x1 +c2*x2 +c3*x3  +c4*x4   +c5*x5  +c6*x6=0
where a1..a6 and b1....b6, c1...c6 are all known and x1 x2...x6 are
all unknowns .....my problems is more like an undetermined linear
equations but i have some constraints to the solutions of the
x1......x6
1)all the x1 ....x6 must be positive and non zero
2)x1^2 +x2^2........+x6^2=6

is actually
i have a matrix
a1*x1 +a2*x2 +a3*x3  +a4*x4   +a5*x5  +a6*x6=0
b1*x1 +b2*x2 +b3*x3  +b4*x4   +b5*x5  +b6*x6=0
c1*x1 +c2*x2 +c3*x3  +c4*x4   +c5*x5  +c6*x6=0
where a1..a6 and b1....b6, c1...c6 are all known and x1 x2...x6 are
all unknowns .....my problems is more like an undetermined linear
equations but i have some constraints to the solutions of the
x1......x6
1)all the x1 ....x6 must be positive and non zero
2)x1^2 +x2^2........+x6^2=6

i was responsed by Sir Daniel to try to solve this problem like
In:= coeffs = RandomInteger[{-10,10}, {3,6}]
Out= {{4, 5, -5, -4, 5, -7}, {-1, 3, -1, -4, 4, 2},
{-9, 3, 10, -1, -2, 8}}

vars = Array[x,6];
eqns = Map[#==0&, Join[coeffs.vars, {vars.vars-6}]];
ineqs = Map[#>0&, vars];

In:= InputForm[system = Join[eqns, ineqs]]
Out//InputForm=
{4*x + 5*x - 5*x - 4*x + 5*x - 7*x == 0,
-x + 3*x - x - 4*x + 4*x + 2*x == 0,
-9*x + 3*x + 10*x - x - 2*x + 8*x == 0,
-6 + x^2 + x^2 + x^2 + x^2 + x^2 + x^2 == 0,
x > 0, x > 0, x > 0, x > 0, x > 0, x > 0}

In:= InputForm[soln = Simplify[FindInstance[system, vars]]]
Out//InputForm=
{{x -> (-15179 + 562*Sqrt)/533792,
x -> (-30679 + 71*Sqrt)/266896,
x -> (87391 + 292*Sqrt)/533792,
x -> 7/32, x -> 3/8,
x -> (-38667 + 214*Sqrt)/533792}}

MY problem still remains when i try to use the matrix
{{-0.2162,0.1283,-0.09057,0.07089},{0.1283,-0.09057,0.07089,-0.05793},{-0.09057,0.07089,-0.05793,0.04855}}

it gives me no solution however i know that solution exsists .......can anyone help me ,should i go for non linear programming ,but it stills doesnt gaurantee me a feasible solution ...what should i do

THnax

```

• Prev by Date: Sub-division of surface plotting by PlotPoints->{m,n}
• Next by Date: Helpfile-like action in ordinary notebook
• Previous by thread: Re: Sub-division of surface plotting by PlotPoints->{m,n}
• Next by thread: Helpfile-like action in ordinary notebook