this->P[0] = (-1*Cp*ro*dR/dTau)*((NT1(tS, tS)*N1_1()*C1(w1.r, w2.r)) + (NT2(tS, tS)*N2_1()*C2(w1.r, w2.r)));

    if (w2.isLast == true) {
        this->P[1] = ((-1*Cp2*ro2*dR/dTau)*((NT1(tS, tS)*N2_1()*C1(w1.r, w2.r)) + (NT2(tS, tS)*N2_2()*C2(w1.r, w2.r)))) - (2 * alfa*rmax*tO);
    }
    else {
        this->P[1] = (-1*Cp2*ro2*dR/dTau)*((NT1(tS, tS)*N2_1()*C1(w1.r, w2.r)) + (NT2(tS, tS)*N2_2()*C2(w1.r, w2.r)));
    }