Bug: Attempting to calculate a double sum crashes the CAS

[quinyu]

Calculating the following double sum causes a calculator crash:
   
While I understand that the closed form is nasty (involves the Psi function, plus the Euler-Mascheroni constant), the Prime could just say "no idea" instead of crashing. (Solution from Wolfram Alpha. TI nSpire CAS just returns the expression.)
   
This is Software Version 2015 6 17 (8151), HW ver C, CAS version 1.1.2-11, OS 0.037.526

[quinyu]

HP 50g delivers a half-solution, giving
   

[DrD]

Does this work on Wolfram, Matlab, et al?

[quinyu]

Wolfram Alpha gives identical result. (Sorry, earlier I thought that Mathematica would too - no, see below.) Maple 18 gives:
    (some little editing to keep style consistent, and factor things a little)

MuPAD gives the same half-solution as the HP 50g. Matlab R2015a also.

[quinyu]

Best sum that I could find with Mathematica 10.0 is:
    , where H_n is the nth harmonic number. This is the one I like the most, so far.