>

restart;

 

>	intr:=3.5; Dnull:=12999; lambda:= 1 + intr/1200;

 

 

 

	(1)

 

>	f:= n-> lambda^n *Dnull-475*(1-lambda^n)/(1-lambda);

 

(2)

 

>	plot(f(n), n=1..60);

 

>

 

>	intr2:=4.5; Dnull2:=12050; lambda2:= 1 + intr2/1200;

 

 

 

	(3)

 

>	f2:= n-> lambda2^n *Dnull2-475*(1-lambda2^n)/(1-lambda2);

 

(4)

 

>	plot(f2(n), n=1..60);

 

>

 

>

 

>	intr3:=6.5; Dnull3:=11900; lambda3:= 1 + intr3/1200;

 

 

 

	(5)

 

>	f3:= n-> lambda3^n *Dnull3-475*(1-lambda3^n)/(1-lambda3);

 

(6)

 

>	plot([f(n), f2(n),f3(n)], n=1..60, color=[blue,black,red]);

 

>

 

>	nterm:=28; nterm2:=26; nterm3:=26; P := nterm*475 + 1000 + f(nterm); P2:= nterm2*475 + 1500 + f2(nterm2); P3 := nterm3*475 + 500 + f3(nterm3);

 

 

 

 

 

 

	(7)

 

>