function [ f ] = point_source( value, position, nx, dx )
f = zeros(nx,1);
i = int16(position/dx);
f(i) = value;
end

