How to write c code for density map