from scipy.stats import norm
sigma=24.6
avg=56.6
N=48
point=77
rank=N-norm.cdf(x=point,loc=avg,scale=sigma)*N
print(rank)
9.76686260423945
rank=10
point=norm.ppf(q=(1-rank/N),loc=avg,scale=sigma)
print(point)
76.58055791689786
9.76686260423945
76.58055791689786