Re: Numerical simulation
- To: mathgroup at smc.vnet.net
- Subject: [mg113703] Re: Numerical simulation
- From: Peter Pein <petsie at dordos.net>
- Date: Mon, 8 Nov 2010 03:39:45 -0500 (EST)
- References: <iaonn6$jm9$1@smc.vnet.net>
Am Tue, 2 Nov 2010 10:04:22 +0000 (UTC)
schrieb Equilib <petterbe at gmail.com>:
> Hi,
> Does anybody have a good idea on how to set up and run simple
> numerical simulations in Mathematica?
>
> I have an equation that consists of the variables p, g that I want to
> solve for p for different values of g=[0,1] and c=[0,1]. The function
> could be for example 1/2+c+g(1-p)+root(c+p*g-1)=0. (the real equation
> is to complex to solve algebraicly using Solve[XX,p]).
>
> I guess I need to construct a matrix for different values [0,1] of g
> and c and then use this to numerically solve for p. but how do i set
> this up in a smart way? And how do I use the results to plot p as a
> function of for example g?
>
> Best,
> Petter
>
>
Hi,
when you got (or will get) an equation which can be handled by Reduce,
you can try:
func[c_, g_] = Block[{listlist=
LogicalExpand[Reduce[1/2 + c + g*(1 - p) + Sqrt[c + p*g - 1] == 0
&& 0 <= c <= 1 >= g >=0 == Im[p], p]] /. Or -> List /.
HoldPattern[(c1___) && p == (pp_) && (c2___)] :> {pp, c1 && c2}},
Piecewise[{FullSimplify[#1, #2], FullSimplify[#2]}& @@@ listlist]
]
Peter