噴泉碼小程序

clc;clear;
K = 10000;
c = 0.2;
o = 0.05;
S = fix(c*sqrt(K)*log(K/o));
for d = 1:50  
    %*******************************P(d)開始************************************
    index_d = 2:K;
    if find(index_d==d)
        P = 1/(d*(d-1));
    elseif d==1
        P = 1/K;
    end
    P;
    PP(d) = P;
    %*******************************P(d)結束*******************************
    %*******************************t(d)開始*******************************
    a = 1:round(K/S)-1;
    if find(a==d)
        t = (S/K)/d;
    elseif d==round(K/S)
        t = (S/K)*log(S/o);
    else
        t = 0;
    end
    t;
    tt(d) = t;
end
%*******************************t(d)結束************************************
stem(1:50,PP)
hold on;
stem(1:50,tt)