Although elastic properties and their anisotropy play a crucial role in materials science and technology, it is not straightforward to visualize how properties such as shear modulus vary with the direction of stress/strain. Here, written in GNU Octave, we have developed an easy-to-use, flexible, and user-friendly open-source toolbox with a graphical user interface, VELAS, which automatically visualizes and analyses the elastic anisotropy of arbitrary crystal systems using second-order elastic constants. Using Voigt-Reuss-Hill averaging scheme, VELAS allows the calculation of more than ten significant mechanical properties such as modulus, hardness, anisotropy index, ductility, and bond type. It also supports the determination of the mechanical stability of crystals at atmospheric and high pressures using the Born mechanical stability criterion. VELAS provides a new easy-to-use tool for detecting and identifying unusual mechanical properties such as negative Poisson's ratio, negative linear compressibility, and negative bulk modulus. VELAS offers several alternative visualization solutions for elastic properties, such as unit spherical projection, and map projection, and supports direct output of high-quality images. Furthermore, VELAS provides users an interface to read data from the Materials Project API in both online and offline modes. After introducing the basic theory, we detail the software framework, technical route, and installation specifications of VELAS. Moreover, the reliability and versatility of VELAS are confirmed by the analysis of the cases of negative linear compressibility, negative Poisson's ratio, hardness, and fracture toughness. Program summaryProgram Title: VELASCPC Library link to program files:https://doi.org/10.17632/bj27474v6f.1Developer's repository link:https://github.com/ranzhengcode/VELASLicensing provisions: GNU General Public License 3Programming language: GNU OctaveNature of problem: To determine the mechanical stability of any crystal at atmospheric and high pressures. To automatically visualize and analyze the elastic anisotropy of arbitrary crystals and identify unusual elastic properties.Solution method: Firstly, the determination of the mechanical stability of crystals at atmospheric and high pressures using the Born mechanical stability criterion. Secondly, more than ten significant mechanical properties such as elastic moduli, ratios, hardness, and anisotropy index are analyzed using Voigt-Reuss-Hill averaging scheme. Then, the mechanical behavior of crystals in the elastic regime is derived from a comprehensive tensor analysis of the second-order elastic constants. Finally, Pugh Ratio, Hardness, and Fracture Toughness in 3D space are calculated based on our approximate treatment of elastic properties.