同核分子优化算法 Homonuclear Molecules Optimization HMO【matlab代码】

function [Nelectron_Layer]= Electron_arrangement(x) if x==1 % H Nelectron_Layer(1,1) = 1; end if x==2 % He Nelectron_Layer(1,1)=2; end if x==3 % Li Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=1; end if x==4 % Be Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=2; end if x==5 % B Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=3; end if x==6 % C Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=4; end if x==7 % N Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=5; end if x==8 % O Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=6; end if x==9 % F Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=7; end if x==10 % Ne Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; end if x==11 % Na Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=1; end if x==12 % Mg Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=2; end if x==13 % Al Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=3; end if x==14 % Si Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=4; end if x==15 % P Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=5; end if x==16 % S Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=6; end if x==17 % Cl Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=7; end if x==18 % Ar Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=8; end if x==19 % K Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=8; Nelectron_Layer(4,1)=1; end if x==20 % Ca Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=8; Nelectron_Layer(4,1)=2; end if x==21 % Sc Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=9; Nelectron_Layer(4,1)=2; end if x==22 % Ti Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=10; Nelectron_Layer(4,1)=2; end if x==23 % V Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=11; Nelectron_Layer(4,1)=2; end if x==24 % Cr Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=13; Nelectron_Layer(4,1)=1; end if x==25 % Mn Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=13; Nelectron_Layer(4,1)=2; end if x==26 % Fe Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=14; Nelectron_Layer(4,1)=2; end if x==27 % Co Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=15; Nelectron_Layer(4,1)=2; end if x==28 % Ni Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=16; Nelectron_Layer(4,1)=2; end if x==29 % Cu Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=1; end if x==30 % Zn Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=2; end if x==31 % Ga Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=3; end if x==32 % Ge Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=4; end if x==33 % As Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=5; end if x==34 % Se Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=6; end if x==35 % Br Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=7; end if x==36 % Kr Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=8; end if x==37 % Rb Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=8; Nelectron_Layer(5,1)=1; end if x==38 % Sr Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=8; Nelectron_Layer(5,1)=2; end if x==39 % Y Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=9; Nelectron_Layer(5,1)=2; end if x==40 % Zr Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=10; Nelectron_Layer(5,1)=2; end if x==41 % Nb Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=12; Nelectron_Layer(5,1)=1; end if x==42 % Mo Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=13; Nelectron_Layer(5,1)=1; end if x==43 % Tc Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=13; Nelectron_Layer(5,1)=2; end if x==44 % Ru Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=15; Nelectron_Layer(5,1)=1; end if x==45 % Rh Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=16; Nelectron_Layer(5,1)=1; end if x==46 % pd Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; end if x==47 % Ag Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=1; end if x==48 % Cd Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=2; end if x==49 % In Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=3; end if x==50 % Sn Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=4; end if x==51 % Sb Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=5; end if x==52 % Te Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=6; end if x==53 % I Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=7; end if x==54 % Xe Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=8; end if x==55 % Cs Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=8; Nelectron_Layer(6,1)=1; end if x==56 % Ba Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=8; Nelectron_Layer(6,1)=2; end if x==57 % La Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=18; Nelectron_Layer(5,1)=9; Nelectron_Layer(6,1)=2; end if x==58 % Ce Nelectron_Layer(1,1)=2; Nelectron_Layer(2,1)=8; Nelectron_Layer(3,1)=18; Nelectron_Layer(4,1)=20; Nelectron_Layer(5,1)=8; Nelectron_Layer(6,1)=2; end if x==59 % Pr Nelectron_Layer(1,1)=2; Nelectron_Lay