% Calculates wind speed extremes by Lieblein's BLUE for Mean Return 
% Intervals M from user's time series V of annual wind speed max
% speed = vector of extreme wind speeds with MRI M 
% M(i)does not have to be an integer 
% mu, sigma = location and scale parameters of Gumbel distribution
addpath('C:\windPRESSURE');
filename = 'WASP.csv'; % Name of Excel file of wind annual speeds
xlrange = 'A1:A22';             
V = xlsread(filename,xlrange);
M = [10 20 25 50 100 300 700 1700 3000]
[speed, mu, sigma] = blue8(V,M)

