4 Algorithm implementation based on modified Saunders and Smith algorithm

4.1 introduction
4.2 Worked examples