  We consider two loop sunset diagrams with two mass scales m and M at the
threshold and pseudotreshold that cannot be treated by earlier published
formula. The complete reduction to master integrals is given. The master
integrals are evaluated as series in ratio m/M and in epsilon with the help of
differential equation method. The rules of asymptotic expansion in the case
when q^2 is at the (pseudo)threshold are given.
