An Iterated Local Search Algorithm for Estimating the Parameters of the Gamma/Gompertz Distribution

<table><tr><td><table class="algorithm" id="alg1"><tr><td colspan="2"><i>s </i><sub>
                0</sub><i>←</i> Generate Initial Solution (<i>s</i>)</td></tr><tr><td colspan="2"><svg height="13.95" id="M50" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 14.1125 13.95" width="14.1125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,13.688)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,5.525)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg><i>←</i> Apply Local Search (<i>s </i><sub>0</sub>)</td></tr><tr><td colspan="2"><i>While</i> termination condition not met <i>do </i></td></tr><tr><td colspan="2"><i>  </i><svg height="17.9125" id="M51" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 10.65 17.9125" width="10.65" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.65)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,9.488)"><path d="M227 744l-123 -338l-31 15l73 368q12 3 41.5 -8t36.5 -20z" id="x2032"></path></g>
</svg>              <i>←</i>  Apply perturbation (<svg height="13.95" id="M52" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 14.1125 13.95" width="14.1125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,13.688)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,5.525)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg>)</td></tr><tr><td colspan="2"><i>  </i><svg height="17.9125" id="M53" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 17.637501 17.9125" width="17.637501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.65)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,9.488)"><path d="M227 744l-123 -338l-31 15l73 368q12 3 41.5 -8t36.5 -20z" id="x2032"></path></g><g transform="matrix(.012,-0,0,-.012,9.971,9.488)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg>  <i>←</i>  Apply Local Search (<svg height="17.9125" id="M54" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 10.65 17.9125" width="10.65" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.65)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,9.488)"><path d="M227 744l-123 -338l-31 15l73 368q12 3 41.5 -8t36.5 -20z" id="x2032"></path></g>
</svg>)</td></tr><tr><td colspan="2"><i>  </i><svg height="13.95" id="M55" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 14.1125 13.95" width="14.1125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,13.688)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,5.525)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg>         <i>←</i>  Apply acceptance criterion (<svg height="13.95" id="M56" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 14.1125 13.95" width="14.1125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,13.688)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,5.525)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg><i>,  </i><svg height="17.9125" id="M57" style="vertical-align:-0.13794pt" version="1.1" viewbox="0 0 17.637501 17.9125" width="17.637501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.65)"><path d="M352 391q0 -31 -27 -44q-14 -7 -24 6q-39 48 -84 48q-23 0 -39.5 -15t-16.5 -40q0 -43 73 -90q49 -32 70 -58t21 -57q0 -58 -62 -105.5t-129 -47.5q-40 0 -75.5 25t-35.5 52q0 28 32 46q7 4 15 3t11 -6q19 -31 48.5 -50.5t54.5 -19.5q34 0 54 19.5t20 42.5q0 43 -65 81
q-97 56 -97 123q0 50 51 96q19 17 58 32.5t62 15.5q37 0 61 -18t24 -39z" id="x1D460"></path></g>
<g transform="matrix(.012,-0,0,-.012,6.438,9.488)"><path d="M227 744l-123 -338l-31 15l73 368q12 3 41.5 -8t36.5 -20z" id="x2032"></path></g><g transform="matrix(.012,-0,0,-.012,9.971,9.488)"><path d="M471 153q-22 -15 -61 -13q-25 31 -45.5 49.5t-56.5 41.5q4 -71 28 -134q-17 -33 -42 -46q-24 12 -42 46q24 63 28 134q-36 -23 -56.5 -41.5t-45.5 -49.5q-36 -2 -61 13q0 28 19 59q65 10 130 43q-65 33 -130 43q-19 31 -19 59q22 15 61 13q25 -31 45.5 -49.5t56.5 -41.5
q-4 71 -28 134q17 33 42 46q24 -12 42 -46q-24 -63 -28 -134q36 23 56.5 41.5t45.5 49.5q36 2 61 -13q0 -28 -19 -59q-65 -10 -130 -43q65 -33 130 -43q19 -31 19 -59z" id="x2217"></path></g>
</svg>)</td></tr><tr><td colspan="2">   <i> </i>Memorize Best Found Solution</td></tr><tr><td colspan="2"><i>End While </i></td></tr></table></td></tr></table>

Modelling and Simulation in Engineering

alg1

Algorithm 1

Algorithm 1: An Iterated Local Search Algorithm for Estimating the Parameters of the Gamma/Gompertz Distribution 