Mp = 56.5/50;
fp = 7.5;

dr1 = sqrt((1-sqrt(1-1/Mp^2))/2);
wn = 2*pi*fp/sqrt(1-2*dr1^2);