Before revision:  Our main method is to use the solving method of a partial differential equation for bond pricing with constant default intensity in every subinterval and to take expectation to remove the random constants.
After revision:  Our main method is to use the solving method of a partial differential equation with a random constant in every subinterval and to take expectation to remove the random constants.