Petrozavodsk State University, St. Petersburg State University
Complex approach to the numerical studying of dynamical systems with two degrees of freedom is suggested. The program to study generalized Henon-Heiles model with the potential $U(x,y)=ax^2+by^2+cxy^2+dy^3$, where $a, b, c, d$ - are parameters of the model, is presented. The program allows to choose the method of numerical integration (Runge-Kutta-Fehlberg, Bulirsch-Stoer, symplectic and others), builds the trajectory, searches for the Lyapunov characteristic numbers, finds Poincare cross-section and finds points of the contours of orbit and folds of the velocity field.