%% Benefit reduction factor for early retirement
function f = benefit_factor_early(age, fra)
    if age < fra
        months_early = round((fra - age) * 12);
        if months_early <= 36
            f = months_early * ((5/9) * 0.01);
        else
            f = (5/9)/100 * 36 + (5/12)/100 * (months_early - 36);
        end
    else
        f = 0; % 0 percent reduction at FRA
    end
end