function [beta,mu] = birth_mortality_rates()

n = 100; % max age
ages = (1:n)';

x = ages(1:n);
sigma_ = 6;
beta = 2/(sigma_*sqrt(2*pi)).*exp(-0.5*(x-32).^2./(sigma_^2));


% Mortility data from Sterbetafel 2010/2012 Statistisches Bundesamt
% average male/female
mu = [
    0.0034
    0.0003
    0.0002
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0001
    0.0002
    0.0003
    0.0003
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0004
    0.0005
    0.0005
    0.0005
    0.0005
    0.0006
    0.0006
    0.0006
    0.0007
    0.0007
    0.0008
    0.0009
    0.0010
    0.0011
    0.0012
    0.0014
    0.0015
    0.0017
    0.0019
    0.0022
    0.0024
    0.0028
    0.0031
    0.0034
    0.0038
    0.0043
    0.0047
    0.0051
    0.0056
    0.0061
    0.0066
    0.0072
    0.0078
    0.0085
    0.0092
    0.0100
    0.0108
    0.0117
    0.0126
    0.0136
    0.0146
    0.0157
    0.0171
    0.0186
    0.0205
    0.0228
    0.0256
    0.0288
    0.0325
    0.0368
    0.0416
    0.0469
    0.0529
    0.0596
    0.0671
    0.0756
    0.0852
    0.0961
    0.1083
    0.1219
    0.1369
    0.1532
    0.1708
    0.1896
    0.2093
    0.2298
    0.2508
    0.2720
    0.2932
    0.3140
    0.3341
    ];