
function grid = UNEQspace(Gmin,Gmax,nG)

Del     = (Gmax-Gmin)/(nG-1)^2;
grid    = zeros(nG,1);
grid(1) = Gmin;
for i=2:nG
    grid(i) = grid(1) + Del*(i-1)^2;
end