function x= mesh1D(a,b,nq,C)
% x discretisation non régulière avec nq points
% obtenu en perturbant légèrement la discrétisation régulière à nq points
% 
if nargin==3,C=1/4;end % default value
assert(C<1/2 && C>=0)
h=(b-a)/(nq-1);
x=a:h:b; % discretisation régulière
per=C*h*(-1+2*rand(1,nq-2)); % perturbation entre -C*h et C*h pour les pts interieurs
x(2:end-1)=x(2:end-1)+per;
 
