figure('Name','fig')
for i=1:100
    vec=[0]
    t=59
    x=0
    claim=0
    sum=0
    k=0
    while x<t
        k=k+1
        if rand<=(2+cos(pi*x/182.5))/3
            vec=[vec x];
            sum=sum+lognrnd(-2,1)
            claim=[claim sum]
        end
        x=x-log(rand)/3
    end
    plot(vec,(10^6)*claim)
    hold on
end
title('100 realizations claims amounts (59 days each)')
xlabel('time')
ylabel('total claim amount'