1-D SPH Code, Matlab

Hi all,

I want to write a 1-d sph matlab code but i am lost on a number on issues and hoping someone can help.

How do i assign boundary conditions to particles? (i.e. i want to model a solid mechanics problem whereby a rod is impacted at one end and held at the other)

In regards to the initial conditions how to i assign that to the particles? (i.e for the particle that is intially impacted, do i assign that an initial velocity and keep all other particles at 0)

And time stepping i'm alittle lost as to this process in regards to updating the particles as time progresses.

Anyones help would be much appreciated.

Many thanks