f<-function(X){ P0<-300000 P<-c() P[1]<-P0+8000+18000-X for (i in 2:10){ P[i]<-P[i-1]+0.08*min(10^5,P[i-1])+0.09*max(0,P[i-1]-10^5)-X } return(P) } f(10^4) f(45797)