Solving algebraic equations with some conditions

• To: mathgroup at smc.vnet.net
• Subject: [mg92007] Solving algebraic equations with some conditions
• From: Kay-Michael Voit <kay at voits.net>
• Date: Tue, 16 Sep 2008 19:26:36 -0400 (EDT)

```Hello,
I am trying to solve a system of 3 PDEs using a Fourier ansatz:
nL[x_, t_] =
Sum[cn[p, l] Exp[I (p k x + l \[Omega] t)], {p, -1, 1}, {l, -1, 1}];
EscL[x_, t_] =
Sum[cE[p, l] Exp[I (p k x + l \[Omega] t)], {p, -1, 1}, {l, -1, 1}];
iL[t_] = Sum[ci[l] Exp[I l \[Omega] t], {l, -1, 1}];

So I have a system of up to 55 algebraic equations (Not every Exp[I (p k
x + l \[Omega] t)] occures after inserting the ansatz, thus some will be
trivial) for the 55 Fourier coefficient. The equations are not linear,
but linear in every coefficient. I was told that this would take some
time but should be possible.

Additionally, I know that the functions are all real, thus I get more
equations:
K[p,l]=Conjugate[K[p,l]]

Unfortunately, I have the impression that Mathematica cannot handle
these in Solve[], as simple examples produce things like Conjugate^-1[I].
A possibility I considered is constructing the coefficients with real
and imaginary part:

iL[t_]
= Sum[(ci[l, re] + ci[l, im] I) Exp[I l \[Omega] t], {l, -1, 1}];

Then I had to give Mathematica assumtions that these are all real so
that an expression like

-H Wges+(I cn[0,0,im]+cn[0,0,re])/\[Tau]==0

will be correktly solved to

{cn[0,0,im] -> 0, cn[0,0,re]->\[Tau] H Wges}
(also assuming that \[Tau], H, Wges are real)
I am not sure how to achieve that.

Indeed I have another condition for the coefficients:
K[p,l][\[omega]]=K[p,-l]][-\[omega]]
In this case, I have completely no idea how to model this in Mathematica.

Can anybody give me a hint (on the problem or on literature) how to
solve this problem?