Optimization of Parameters in the Generalized Extreme-Value Distribution Type 1 for Three Populations Using Harmonic Search