We address the multi-view shape from shading problem, that is the recovery of 3-D shape, lighting configuration and surface albedo from multiple calibrated views. Previous approaches to this problem relied on physically impossible illumination models (negative light) and only work on constant albedo and resulted in biased estimates of shape and lighting positions. Furthermore, since the solution involves infinite-dimensional optimization, existing approaches were quite slow. We develop a new model that explicitly enforces positivity in the light sources with the assumption that the object is Lambertian and its albedo is piecewise constant and show that the new model significantly improves the accuracy and robustness relative to existing approaches. Furthermore, we show that the most computationally expensive step in the optimization can actually be solved in closed form. This significantly improves speed of convergence over existing schemes. We illustrate the behavior of our algorithm directly on the same data used by previous authors, so direct comparison is possible.