function rt = rt400(jaar,beta)
    % Geef een LS waarde terug van een jaartal (waarbij 1972 van wordt
    % afgetrokken) en parameter vector voor dat jaartal
    rt = 0;
    for i=1:size(beta)
        rt = rt + beta(i) .* (jaar ^ (i - 1));
    end
