for k = 0 : 1

    t_kT = [-A*T : Ts : A*T + k*T];
    t_kt1 = [A*T : Ts : A*T + k*T];
    ph1_B = [ph1 zeros(1, length(t_kt1))];
    ph1_kT = [zeros(1, length(t_kt1)) ph1_B];

    figure;
    plot(t_kT,ph1_B);
    hold on;
    plot(t_kT,ph1_kT);

    figure;
    plot(t_kT,(ph1_B.*ph1_kT));
end