ModellingChapter2 -- Sage

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)

ModellingChapter2