figure('Name','Insurance claims as a function of time')
t=59
for i=1:100
    vec=0
    x=0
    while x<t
        if rand<=(2+cos(pi*x/182.5))/3
            vec=[vec x];
        end
        x=x-log(rand)/3
    end
    plot(vec,1:1:numel(vec))
    hold on
end
title('Insurance claims as a function of time')
xlabel('days')
ylabel('total claims')