Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

I have this system with 3 equations and 3 unknowns:

eqs = {q*(q - a) + u*(-r + 1) (p^r) (j^g)*(b - q) + 
        r*s*u*(p^(-r - 1)) (j^g) (b - q) == 0, 
       p*(2*q - a) - u*(p^(-r - 1)) (j^g) + (b - q) (s - b) + 
        s*u*(p^(-r)) (j^g) - s (b - 2*q) - c*(b - a) == 0, 
       u*(p^(-r)) g (j^(g - 1))*(b - q) - 
        s*u*g*(p^(-r)) (j^(g - 1)) (b - q) - (b - a) == 0} 

The unknown variables are {p, q, j}.

How can I solve this system, given that

Solve[{q*(q - a) + u*(-r + 1)*(p^r)*(j^g)*(b - q) +
        r*s*u*(p^(-r - 1))*(j^g)*(b - q) == 0,
       p*(2*q - a) - u*(p^(-r - 1))*(j^g) + (b - q)*(s - b) +
        s*u*(p^(-r))*(j^g) - s*(b - 2*q) - c*(b - a) == 0,
       u*(p^(-r))*g*(j^(g - 1))*(b - q) -
        s*u*g*(p^(-r))*(j^(g - 1))*(b - q) - (b - a) == 0
      }, {p, q, j}]

doesn't seem to work?

share|improve this question

closed as too localized by belisarius, Sjoerd C. de Vries, m_goldberg, Silvia, Yves Klett May 6 '13 at 9:08

This question is unlikely to help any future visitors; it is only relevant to a small geographic area, a specific moment in time, or an extraordinarily narrow situation that is not generally applicable to the worldwide audience of the internet. For help making this question more broadly applicable, visit the help center.If this question can be reworded to fit the rules in the help center, please edit the question.

    
The syntax is Solve[eqs,{p,q,j}], but it seems too hard as it stands. Do you have further restrictions like g being integer or alike? –  Rolf Mertig May 3 '13 at 7:06
    
You appear to be using Matlab functions instead of Mathematica functions. Which system are you using? –  Jonathan Shock May 3 '13 at 7:39
    
no , unfortunately I haven't got any information about "r" and "g" and "u" ... –  Rozita May 3 '13 at 7:40
    
If Matlab can solve this system , please tell me , which function or code I should use . –  Rozita May 3 '13 at 7:42
    
Are you using Mathematica? –  cormullion May 3 '13 at 7:46
show 4 more comments

1 Answer 1

If you specify some integer value for g and r then you get the solution for your system

  LE[g_, r_] := Module[{eqs},

  eqs = {q*(q - a) + u*(-r + 1) (p^r) (j^g)*(b - q) + 
  r*s*u*(p^(-r - 1)) (j^g) (b - q) == 0, 
p*(2*q - a) - u*(p^(-r - 1)) (j^g) + (b - q) (s - b) + 
  s*u*(p^(-r)) (j^g) - s (b - 2*q) - c*(b - a) == 0, 
u*(p^(-r)) g (j^(g - 1))*(b - q) - 
  s*u*g*(p^(-r)) (j^(g - 1)) (b - q) - (b - a) == 0};
 Solve[eqs, {p, q, j}]
   ]

Like for example

 LE[1,0]

But for higher values of g and r it is taking a lot of time and $mma$ ends up saying "A very large output is generated.."

share|improve this answer
    
Thanks a lot for your attention ... If I change my model to the uniform distribution , I think the equations become easier. But I don't know which code I should use ... –  Rozita May 3 '13 at 13:50
add comment

Not the answer you're looking for? Browse other questions tagged or ask your own question.