ModellingChapter2

483 days ago by tazzalenghe

def a(n): return (1.+1/(n+1)) def b(n): return 100. 
       
(0.0, 100.0, 250.0, 17945.9472255)
(0.0, 100.0, 250.0, 17945.9472255)
pp = vector(RDF, 100) pp[0]= 0. for i in range(1, 51): pp[i] = a(i-1)*pp[i-1] + b(i-1) pp[0], pp[1], pp[2], pp[50] 
       
(0.0, 100.0, 250.0, 17945.9472255)
(0.0, 100.0, 250.0, 17945.9472255)
def p(n): if n==0: return RDF(0) else: return RDF(a(n-1)*p(n-1) + b(n-1)) p(0), p(1), p(2), p(50) 
       
(0.0, 0.0)
(0.0, 0.0)
p(100) 
       
42392.5129282
42392.5129282
i = .06 j = .03 def a(n): return RDF(1+i) def b(n): if n==0: return RDF(100) else: return RDF(b(n-1)*(1+j)) 
       
p(0), p(1), p(2), p(50) 
       
(0.0, 100.0, 209.0, 46787.4941876)
(0.0, 100.0, 209.0, 46787.4941876)