...dgesl.f
Alternatively, you can do this excercise with Maple, Mathematica, Matlab, or a similar symbolic manipulation package