function [age sur den]=survivalrate


death=[0.006713
0.000444
0.000300
0.000216
0.000179
0.000168
0.000156
0.000143
0.000125
0.000103
0.000086
0.000088
0.000125
0.000206
0.000317
0.000438
0.000552
0.000657
0.000747
0.000825
0.000905
0.000983
0.001033
0.001049
0.001038
0.001019
0.001006
0.000998
0.001002
0.001018
0.001042
0.001072
0.001113
0.001156
0.001212
0.001276
0.001355
0.001456
0.001585
0.001739
0.001903
0.002077
0.002268
0.002479
0.002706
0.002943
0.003190
0.003453
0.003741
0.004057
0.004405
0.004778
0.005166
0.005554
0.005939
0.006335
0.006760
0.007234
0.007796
0.008470
0.009282
0.010204
0.011178
0.012118
0.013024
0.013999
0.014995
0.016161
0.017527
0.019109
0.020890
0.022925
0.025280
0.027972
0.030997
0.034386
0.038027
0.042036
0.046447
0.051297
0.056623
0.062465
0.068867
0.075871
0.083524
0.091872
0.100962
0.110842
0.121558
0.133155
0.145675
0.159156
0.173631
0.189127
0.205661
0.223242
0.241869
0.261527
0.282188
0.303810
1.00000];

age=(1:size(death,1))-1;

sur(1)=1;
 
for i=2:size(death,1)
sur(i)=sur(i-1)*(1-death(i,1)); 


end

for i=2:size(death,1)
den(i)=sur(i-1)*death(i,1); 


end


 