function [ C,C0 ] = basic_model( x_limits, nx, dx )

xs = x_limits(1) + [0:nx-1]*dx;
C0 = ones(nx,1)';
dC = -100.0*(xs-0.5) .* exp(-((xs-0.5).^2)/(1e-4));

for i = 1:size(dC,1)
    if abs(dC(i)) < 1e-7
        dC(i) = 0;
    end

C = C0 + dC;

end

