clear all
close all

windSpeed-kts = 20.0
windSpeed-mps = windSpeed-kts * 1852.0 / 3600

g = 9.8;
U195 = windspeed-mps;
T1 = 10

f =  0.02:.01:0.5;
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,f,modPM_Sw);

figure
plot(f(30:end),PM_Sw(30:end),f(30:end),modPM_Sw(30:end));

