This text aims to develop algorithms of shape-preserving spline approximation for curves/surfaces with automatic choice of the tension parameters. Detailed descriptions of algorithms are given, with a strong emphasis on their computer implementation.