[gdal-dev] Computing transformation matrix from set of points

[ogr user]

Sorry to butt in with a question that must have been answered a few
times before.

I am looking for a solver (ready made or an algorithm to implement) to
derive an affine transformation matrix for a given coordinate
transformation based on a set of points with coordinates in each system.

Any help or ideas would be very much appreciated!
--Mike