age <- c(0,0.25,0.5,0.75,1:96)
live <- c(11650,10310,9756,9203,
               8650,7283,6781,6446,6249,6065,5925,5815,5735,
          5675,5623,5573,5523,5473,5423,5373,5320,5262,5199,
          5132,5060,4985,4910,4835,4760,4685,4610,4535,4460,
          4385,4310,4235,4160,4085,4010,3935,3860,3785,3710,
          3635,3559,3482,3404,3326,3248,3170,3092,3014,2936,
          2857,2776,2694,2612,2530,2448,2366,2284,2202,2120,
          2038,1956,1874,1793,1712,1632,1552,1472,1392,1312,
          1232,1152,1072, 992, 912, 832, 752, 675, 602, 534,
           469, 406, 346, 289, 234, 186, 145, 111,  83,  62,
            46,  34,  24,  16,   9,   4,   1)
detr <- -c(diff(live),-1)
total <- sum(live[c(1,5:100)])
totaldetr <- sum(detr)