просто... Мда, видимо соображал тупо, а то бы догадался.
Проверил М[ax,by]=a*b*M[x,y] - сходится.

И D[ax,by]=(a*b)^2*D[x,y] - тоже сходится.
Попробовал разобрать в общем M[(ax+b)*c,(u*y+v)*l] и вот, что получилось (a, b, c, u, v, l - постоянные)
=cl*(p11*(ax1+b)*(uy1+v)+p12*(ax1+b)*(uy2+v)+...+p1m*(ax1+b)*(uym+v)+...pnm*(axn+b)*(uym+v))=cl*(p11*(aux1y1+avx1+buy1+bv)+...)=cl*(au*(p11*x1*y1+...+pnm*xn*ym)+avm*(x1*(p11+...+p1m)+...+xn*(pn1+...+pnm))+bun*(y1*(p11+...+pn1)+...+ym*(p1m+...+pnm))+bv)
мисс_граффити
А разве у программистов в программе нет теории вероятностей?