％rng（14，'twister'）％取消注释前一行获取相同的输出ms = multiStart;opts = Optimoptions（@ Fmincon，'算法'，'内部点'）;六分= @（x）（4 * x（1）^ 2  -  2.1 * x（1）^ 4 + x（1）^ 6/3 ... + x（1）* x（2） -  4 * x（2）^ 2 + 4 * x（2）^ 4）;问题= createOptimprobrom（'fmincon'，'x0'，[ -  1,2]，...'目标'，六分钟，'lb'，[ -  3，-3]，'UB'，[3,3]，...'选项'，选择）;[Xminm，Fminm，Flagm，Outptm，Manyminsm] =运行（MS，问题，50）;posscolors ='kbgcrm';保持i = 1：这个行的大小（minyminsm，2）％cidx = REM（I-1，长度（POSSCOLORS））+ 1;Color = POSSCOLORS（CIDX）;％绘图开始点U = monyminsm（i）.x0;x0thismin =重塑（[u {：}]，2，长度（u））; plot(x0ThisMin(1, :), x0ThisMin(2, :), '.', ... 'Color',color,'MarkerSize',25); % Plot the basin with color i plot(manyminsm(i).X(1), manyminsm(i).X(2), '*', ... 'Color', color, 'MarkerSize',25); end % basin center marked with a *, start points with dots hold off