[Help-gsl] how to get a matrix to transform point set P to point set Q?

[yaoyansi]

hi, all
I am new to GSL, and I have a question, but I don't know how to solve this 
problem by GSL.
My question is:
Given two point sets P={p0,...,pn} and Q={q0,...,qn},  I want to get a matrix M 
to transform P to Q'={q0',...,qn'}, that is:
      Q'=P*M
The constraint condition is let the SUM has the minimum value:
SUM=|q0-q0'|^2 + |q1-q1'|^2 + ... + |qn-qn'|^2;

Thank you.

Best regards
yaoyansi