This report presents a method for recovery of shape from shading in the case of a non-Lambertian surface illuminated by only a single light source. In the method, we derive simultaneous equations that should be satisfied by parameters determining reflectance map. After identifying the parameters by solving these equations, the surface normal is estimated by the iterative method using the surface normal information on the occluding boundary, and 3-D shape is recovered. The experimental results show the effectiveness of our method for shape from shading.