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 <[email protected]> wrote: > > 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 + > 33*k - 0.33, > lambda s00, s01, s10, s11, k, p:0.6*p*s01 + 0.6*p*s10 + 0.55*p*s11 + > 33*p - 0.33, > lambda s00, s01, s10, s11, k, p:0.11*k*s00 - 0.06*p*s01 - 0.06*p*s10, > lambda s00, s01, s10, s11, k, p:(0.06*k + 0.06*p)*s01 - 0.055*k*s00 - > 0.055*p*s11, > lambda s00, s01, s10, s11, k, p:(0.06*k + 0.06*p)*s10 - 0.055*k*s00 - > 0.055*p*s11] > > > s= findroot(f,(0.1,0.1,0.1,0.1,0.1,0.1)) > print [round(s[k,0],10) for k in range(6)] > > -- > To post to this group, send email to [email protected] > To unsubscribe from this group, send email to > [email protected] > For more options, visit this group at > http://groups.google.com/group/sage-support > URL: http://www.sagemath.org > -- Craig Jolley Foreign Postdoctoral Researcher RIKEN Center for Developmental Biology Laboratory for Systems Biology http://www.cdb.riken.jp/lsb/index.html -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
