g = 9.8;
U195 = 20.6;
T1 = 10

f =  0.075;
w = 2 * pi * f;
alpha = 8.1e-3;
beta = 0.74;
w0 = g / U195;


PM_Sw = ((alpha * g^2) ./ (w.^5)) .* exp((-1*beta)*((w0./w).^4))
modPM_Sw = (0.11/(2*pi))*(((w*T1)/(2*pi))^-5)*exp(-.44*((w*T1)/(2*pi))^-4);
plot(f,PM_Sw);
figure
plot(f,modPM_Sw)