Services & ResourcesWolfram ForumsMathGroup Archive

[Date Index] [Thread Index] [Author Index]

Re: How to access Solve solutions of eq.system as functions

• To: mathgroup at smc.vnet.net
• Subject: [mg126063] Re: How to access Solve solutions of eq.system as functions
• From: Bob Hanlon <hanlonr357 at gmail.com>
• Date: Sat, 14 Apr 2012 03:10:49 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• References: <201204130846.EAA19184@smc.vnet.net>

Clear[n, ne];

sol = With[
    {tnr = 1, te = 1, b = 1, d = 1, f = 1, nemax = 1},
    Solve[{
      g - n/tnr - b n^2 - f n^2 (1 - ne/nemax) + d ne == 0,
      f n^2 (1 - ne/nemax) - ne/te - d ne == 0},
     {n, ne}]] // Simplify;

n[g_] = n /. sol;

ne[g_] = ne /. sol;

Plot[Evaluate[n[g]], {g, -0.5, .5}]

Plot[Evaluate[ne[g]], {g, -0.5, 0.5}]


Bob Hanlon


On Fri, Apr 13, 2012 at 4:46 AM, getjonas <getjonas at gmail.com> wrote:
> i am solving a two equation system with two variables (n, ne) and one
> parameter (g):
>
> Solve[{g - n/tnr - b n^2 - f n^2 (1 - ne / nemax) + d ne == 0,
>
>  f n^2 (1 - ne/nemax) - ne/te - d ne == 0}, {n, ne}]
>
> Here tnr, te, b, d and f are all previously defined numbers.
>
> i get four pairs of analytical solutions, as expected.
>
> The question is, how to access these solutions (n, ne) as functions of
> parameter g? I.e., how to plot n[g], ne[g]?
>

• References:
  • How to access Solve solutions of eq.system as functions
    • From: getjonas <getjonas@gmail.com>