import scipy.stats as stats
import math

x = 5

mean_wt = 15.7*0.5

sd_wt = math.sqrt(15.7*0.5)

prob=stats.norm.cdf((x-mean_wt)/sd_wt)

print('Probabilidad {:.2f}%'.format(prob*100))