home *** CD-ROM | disk | FTP | other *** search
- Xref: sparky sci.math:10468 sci.math.symbolic:2242
- Newsgroups: sci.math,sci.math.symbolic
- Path: sparky!uunet!munnari.oz.au!metro!news
- From: king@physics.su.OZ.AU (Peter King)
- Subject: Solving Rate Equations
- Message-ID: <1992Aug20.224623.7856@ucc.su.OZ.AU>
- Sender: news@ucc.su.OZ.AU
- Nntp-Posting-Host: langmuir.physics.su.oz.au
- Organization: School of Physics, University of Sydney, Australia
- Date: Thu, 20 Aug 1992 22:46:23 GMT
- Lines: 35
-
- I'm interested in solving the following equations (in Mathematica format)
- for a cw optically pumped laser:
-
- eqns = {
- wv N3 - wp (N0 - g0 / g2 N2) == 0
- B21 n (N2 - g2 / g1 N1) - wr N1 == 0
- wp (N0 - g0 / g2 N2) - wr N2 - B21 n (N2 - g2 / g1 N1) == 0
- N0 + N1 + N2 + N3 == N
- B21 n (N2 - g2 / g1 N1) - 2 Pi y n == 0
- }
-
- simultaneously. In Mathematica I naively tried
- Solve[eqns, {N0, N1, N2, N3, n}]
- and ran out of memory on a 40 Mb RAM workstation. I've tried eliminating n
- from the last equation to reduce the system to 4 equations but then I can't
- find the solution to n without running out of memory.
-
- Does anyone know of any tricks/techniques/simplications I can make
- to solve this system? Alternatively does anyone know of simple references
- to Groebner bases?
-
- thanks in advance
-
- Peter
-
- --
- ********************************************************************************
- Peter King * The home of TORTUS.
- Plasma Physics Department * T
- University of Sydney, NSW 2006 * O
- AUSTRALIA * Roid of
- voice: 61 2 692 3371 * The
- FAX: 61 2 660 2903 * University of
- email: king@physics.su.OZ.AU * Sydney
- ********************************************************************************
-