On numerical solution of stochastic partial differential equations of elliptic type