This is C++ code written to simulate the viscoelastic behavior of materials. I always try to improve it in different aspects, including computational efficiency and accuracy. Please feel free to leave comments if you think there is something I could do to improve the code.