This paper tries to implement the Multiphase Particle-In-Cell method in Matlab. The method was developed by Andrews and O'Rourke. It has the advantage of the Euler-Lagrange method, that different particle sizes can be simulated relatively simple. But at the same time this method reduces the computational effort for higher solid fractions because the solid pressure is not calculated for each individual particle, but for each cell. The created program solves one-dimensional problems and deals with different particle sizes and densities. It is possible to define the viscosity, density and flow velocity of the surrounding fluid. The particles are grouped into packages and the number of particles in each package can be specified. Wall collisions are simulated with a factor of restitution. It determines the fraction of the initial speed the particle has after the collision. The first test case deals with the free fall of a single particle. The difference in falling time between the analytical solution and the solution calculated by the program is only half a percent. The sedimentation of a pile of particles is the topic of the second test case. The result shows, that numerical diffusion appears. Particles with different diameters are used in the third test case. It can be concluded that the first results are positive, however a lot of effort is needed to implement correction algorithms for the numerical diffusion. The behaviour of the program when dealing with high solid fractions will be a subject to further research.