on matrix approach was adopted to calculate the basic reproduction number (R0) from the RP model to