The purpose of this paper is to present a numerical method to solve partial stochastic differential equations. This concept remains the differential operator unchanged but discretizes the dimension of the problem. The response function will be decomposed by the Karhunen--Loeve expansion and...