Thanks, this at least gives a reasonable answer and will work in the
short-term. I'm still confused about why solve() didn't work, though.
On Thu, Jul 21, 2011 at 11:44 PM, achrzesz wrote:
>
> Try something like this:
>
> from mpmath import *
> mp.dps = 30; mp.pretty = True
>
> f=[lambda s00, s
Try something like this:
from mpmath import *
mp.dps = 30; mp.pretty = True
f=[lambda s00, s01, s10, s11, k, p:0.55*k*s00 + 0.6*k*s01 + 0.6*k*s10
+ 0.6*p*s01 + 0.6*p*s10 +0.55*p*s11 + 33*s00 + 33*s01 + 33*s10 +
33*s11 - 33.0,
lambda s00, s01, s10, s11, k, p:0.55*k*s00 + 0.6*k*s01 + 0.6*k*s10 +
3
