Finding solution to nonlinear equations numerically in a range

asked 2012-10-15 19:15:30 +0100

Anonymous

Hi,

this question is probably quite easy to solve. I've got two nonlinear equations with a solution in the range of (0,1). As described in the sage tutorial I tried the following:

var('p t')
eq2 = p==1 - (1 - (1 - (1 - 0.01)^1036)^20)*(1 - t)^(9)
eq1 = t==(2*(1 - 2*p))/((1 - 2*p)*16 + 15*p(1 - (2*p)^1023))
solns = solve([eq1,eq2],t, solution_dict=True)
[[s[t].n(10)] for s in solns]

Unfortunately this doesn't work. In addition I tried find_root, but failed miserably. Any ideas?

edit retag flag offensive close merge delete

add a comment

2 Answers

Sort by » oldest newest most voted

answered 2012-10-16 03:24:20 +0100

achrzesz
2803 ●5 ●25 ●61

You can also use mpmath

from mpmath import * 
mp.dps = 30; mp.pretty = True
f1 = lambda t,p:1 - (1 - (1 - (1 - 0.01)**1036)**20)*(1 - t)**(9)-p  
f2 = lambda t,p:(2*(1 - 2*p))/((1 - 2*p)*16 + 15*p*(1 - (2*p)**1023))-t
findroot([f1,f2],(1,1))
[2.39470409963011871605945052144e-309]
[    0.999398871164408619804930822283]

edit flag offensive delete link

add a comment

answered 2012-10-16 01:33:04 +0100

robert.marik
15 ●1 ●7

Hi, the system is special in some sense, I think that you can do the following

var('t')
p=1 - (1 - (1 - (1 - 0.01)^1036)^20)*(1 - t)^(9)
eq1 = t==(2*(1 - 2*p))/((1 - 2*p)*16 + 15*p(1 - (2*p)^1023))
eq1.find_root(-100,100)

and then

p.subs(t=_)

edit flag offensive delete link

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Finding solution to nonlinear equations numerically in a range

2 Answers

Your Answer

Question Tools

Stats

Related questions

Finding solution to nonlinear equations numerically in a range edit

2 Answers

Your Answer

Question Tools

Stats

Related questions

Finding solution to nonlinear equations numerically in a range