This function fits any series of data to the Generalized Extreme Value (GEV) function using the Lmoments method (Hosking, 1990) . The time series should be stored in a text file with the values to be fitted in the first column. A typical application is fitting annual extreme values of streamflow to the GEV function, similar to how such time series are typically fitted to the Log Pearson Type III distribution. The program computes the parameters of the equation x, a, and k.
Q = x + a/k x (1  (ln(P))^k )

A 
B 
C 
D 
E 
F 
1 
x 
a 
k 

T 
Q 
2 
4.886 
0.7608 
0.347 

100 
=$A$2+$B$2/$C$2*(1(1*LN(11/E2))^$C$2) 
3 




50 
11.18 