[gf,gc] = meshgrid(0:0.001:1,2:0.001:10);
Lc = gc./((gc-1).*(gc.*(2-gf)-gf));

% surface in 3D
figure(1);
surf(gf,gc,Lc,'EdgeColor','None');
view(2)
colormap jet;
colorbar