Mathematics and computer science are interested in methods of 2D curve interpolation and extrapolation using the set of key points (knots or nodes). Proposed method, called by author Probabilistic Nodes Combination (PNC), is such a method. This novel PNC method is introduced in the case of Hurwitz- Radon Matrices (MHR). MHR method is based on the family of Hurwitz-Radon (HR) matrices which possess columns composed of orthogonal vectors. Two-dimensional curve is modeled and interpolated via different functions as probability distribution functions: polynomial, sinus, cosine, tangent, cotangent, logarithm, exponent, arcsin, arccos, arctan, arcctg or power function, also inverse functions. It is shown how to build the orthogonal matrix operator and how to use it in a process of curve reconstruction.