basic model working. Adding Dirichlet impulses