I tried to create a program, referring to Minoru Hirose’s Multiple Zeta Values online tools, that expresses multiple zeta values in terms of shuffle-symmetric multiple zeta values and products of multiple zeta values. However, the computation becomes very heavy and times out even for small weights. I recommend downloading the Mathematica file from the my website.
Example: mzv(2,3)