I believe your data gives orbital radii in meters and time in seconds. This means that your result for the slope will come out in SI and you will be able to get a value for Saturn's mass in kilograms (which you can then compare with references in texts or on the 'Net).
You will want to plot T^2 on the y-axis and r^3 on the x-axis. These will then function as effective variables Y = T^2 and X = r^3 , so that your data will fall on a line
Y = mX + b .
Given the form of Newton's extension of Kepler's Third Law (which you quoted), you should find a rather small value for b, the intercept of the line (ideally it would be zero, but you are working with real data) and a tiny value for m (on the order of 10^-15), which is your experimental estimate for [ 4·(pi^2) / GM ] .
You can make the values a bit easier to deal with by dividing your values for T^2 by 10^10 and those for r^3 by 10^25 , which will give numbers that are easier to plot or use in graphing software. This will keep you from getting extremely small output for m and b. You would then have to undo this rescaling by dividing your slope value by 10^15 before extracting the result for M. (The correct value for the line intercept would then be your result for b, multiplied by 10^10.)
EDIT: I fixed some of the exponents, as I was remembering planetary masses in grams, not kilograms. I tried out the plot I described a bit crudely and found that the slope is quite close to 1, making [ 4·(pi^2) / GM ] close to 1·10^-15 . Working carefully, you should get a rather satisfactory result for Saturn's mass in kilograms.