Python - Luminosity of a Star

print ("Calculating Luminosity of a Star")
import math

print ("The formula is: ")
print ("Luminosity = 4 * Pi * (Stellar Radius**2) * Stefan-Boltzmann Constant * (Surface Temperature**4)")
# Stefan-Boltmann Constant is: 5.67 * (10**−8)
# Built in pi: math.pi

#  Input the data
EnteredStellarRadius = input ("Enter the Stellar Radius: ")
EnteredSurfaceTemperature = input ("Enter the Surface Temperature: ")

StefanBoltzmannConstant = ( 5.67*(10**(-8)) )
Pi = math.pi

# Convert entered numbers to float
StellarRadius = float(EnteredStellarRadius)
SurfaceTemperature = float(EnteredSurfaceTemperature)
SurfaceTemperature = float(EnteredSurfaceTemperature)

Luminosity = 4* Pi * (StellarRadius**2) * StefanBoltzmannConstant * (SurfaceTemperature**4) 

print ("The Luminosity of the star is: ", Luminosity)

print ("Thank you to astronomy.osu.edu and fxsolver.com for assistance with this equation.")