I am trying to solve these equations:

**eq1:=cot(theta[n])=(omega[a]^2 - omega[n]^2)/omega[n];**

**eq2:=cot(theta[n]+ omega[n])=-(omega[b]^2 - omega[n]^2)/omega[n];**

for some values of omega[a] and omega[b]. I will in the end create two 3d plots for omega[n] and theta[n] as a funciton of ometga[a] and omega[b]

So for example for some particular values have:

**eqs:=subs(omega[a]=0.2, omega[b]=1, [eq1, eq2]);**

Then solving,

**solve(eqs, [omega[n],theta[n]]);**

gives:

I'm however only interested in solutions where omega[n]>=0, but doing this:

**solve(eqs, [omega[n],theta[n]], UseAssumptions) assuming omega[n]>=0.0;**

returns an empty list. Adding the condition omega[n]>=0.0 to the list of equations also does not work.

Is this a bug or am I missing something? Note, I realize that i can manually go through the entries myself and pick the right solutions, I am just asking how to force maple to do this automagically.

thanks

Here is the full code:

**restart:**

**eq1:=cot(theta[n])=(omega[a]^2 - omega[n]^2)/omega[n];**

**eq2:=cot(theta[n]+ omega[n])=-(omega[b]^2 - omega[n]^2)/omega[n];**

**eqs:=subs(omega[a]=0.2, omega[b]=1, [eq1, eq2]);**

**solve(eqs, [omega[n],theta[n]]);**

**solve(eqs, [omega[n],theta[n]], UseAssumptions) assuming omega[n]>=0.0;**