matlabpool open 2
comsolPort1 = 2036;
system( ['D:\COMSOL43b\bin\win64\comsolserver.exe -np 1 -port ' num2str(comsolPort1) ' &'] );
pause(1)
comsolPort2 = 2037;
system( ['D:\COMSOL43b\bin\win64\comsolserver.exe -np 1 -port ' num2str(comsolPort2) ' &'] );
pause(1)
PopulationSize=100;
Generations=50;
EliteCount=2;
CrossoverFraction=0.5;
nvars=2;
A=[];
b=[];
Aeq=[];
beq=[];
LB=[0.1 10];
UB=[1 30];
nonlcon=[];

options = gaoptimset(@ga);
options = gaoptimset(options,'PopulationSize', PopulationSize,'Generations',Generations,'EliteCount',EliteCount,'MutationFcn', {@mutationadaptfeasible},'CrossoverFraction',CrossoverFraction,'PlotFcns', {@gaplotbestf,@gaplotbestindiv,@gaplotdistance,@gaplotexpectation,@gaplotstopping},'UseParallel', 'always', 'Vectorized', 'off');
[x fval exitflag output population scores]=ga(@inverse2Dparga,nvars,A,b,Aeq,beq,LB,UB,nonlcon,options);