function [xval]=gridx(x) % % maps uniform grid into nonuniform grid % xval = .5/sqrt(.5)*sqrt(x);