In[]:=
​​duncap=3.5;w=1;numnodes=10;​​For[i=1,i<=numnodes,i++,​​Dunbarexport[i]=0;​​Dunbarimport[i]=0;​​For[j=1,j<=numnodes,j++,​​If[i!=j,Dunbarexport[i]+=adj[i,j]];​​If[i!=j,Dunbarimport[i]+=adj[j,i]]]​​]​​Dunbarcons[x_]:=(For[i=1;Dunbarconstraints={},i<=numnodes,i++,Dunbarconstraints=Append[Dunbarconstraints,Dunbarexport[i]+Dunbarimport[i]+adj[i,i]<=2duncap/(1+Tanh[x])]];​​Dunbarconstraints​​);​​(*constraints[x_]:=(Join[Dunbarcons[x],levelcoherenceconstraints,nodeconstraints]);*)​​env[i_,t_]:=Cos[i-wt];​​totaladj[adj_]:=(​​For[i=1;predict=0,i<=numnodes,i++,For[t=1,t<=10,t++,predict+=(adj[i,t]-env[i,t])^2]];​​predict​​);​​Adj=Table[adj[i,j],{i,1,numnodes},{j,1,numnodes}];​​vars=Flatten[Adj];​​sol=NMinimize[{totaladj[adj],Dunbarcons[totaladj[adj]]},vars];
In[]:=
sol
Out[]=
{5.55111×
-15
10
,{adj[1,1]1.,adj[1,2]0.540302,adj[1,3]-0.416147,adj[1,4]-0.989993,adj[1,5]-0.653644,adj[1,6]0.283662,adj[1,7]0.96017,adj[1,8]0.753902,adj[1,9]-0.1455,adj[1,10]-0.91113,adj[2,1]0.540302,adj[2,2]1.,adj[2,3]0.540302,adj[2,4]-0.416147,adj[2,5]-0.989993,adj[2,6]-0.653644,adj[2,7]0.283662,adj[2,8]0.96017,adj[2,9]0.753902,adj[2,10]-0.1455,adj[3,1]-0.416147,adj[3,2]0.540302,adj[3,3]1.,adj[3,4]0.540302,adj[3,5]-0.416147,adj[3,6]-0.989993,adj[3,7]-0.653644,adj[3,8]0.283662,adj[3,9]0.96017,adj[3,10]0.753902,adj[4,1]-0.989993,adj[4,2]-0.416147,adj[4,3]0.540302,adj[4,4]1.,adj[4,5]0.540302,adj[4,6]-0.416147,adj[4,7]-0.989993,adj[4,8]-0.653644,adj[4,9]0.283662,adj[4,10]0.96017,adj[5,1]-0.653644,adj[5,2]-0.989993,adj[5,3]-0.416147,adj[5,4]0.540302,adj[5,5]1.,adj[5,6]0.540302,adj[5,7]-0.416147,adj[5,8]-0.989993,adj[5,9]-0.653644,adj[5,10]0.283662,adj[6,1]0.283662,adj[6,2]-0.653644,adj[6,3]-0.989993,adj[6,4]-0.416147,adj[6,5]0.540302,adj[6,6]1.,adj[6,7]0.540302,adj[6,8]-0.416147,adj[6,9]-0.989993,adj[6,10]-0.653644,adj[7,1]0.96017,adj[7,2]0.283662,adj[7,3]-0.653644,adj[7,4]-0.989993,adj[7,5]-0.416147,adj[7,6]0.540302,adj[7,7]1.,adj[7,8]0.540302,adj[7,9]-0.416147,adj[7,10]-0.989993,adj[8,1]0.753902,adj[8,2]0.96017,adj[8,3]0.283662,adj[8,4]-0.653644,adj[8,5]-0.989993,adj[8,6]-0.416147,adj[8,7]0.540302,adj[8,8]1.,adj[8,9]0.540302,adj[8,10]-0.416147,adj[9,1]-0.1455,adj[9,2]0.753902,adj[9,3]0.96017,adj[9,4]0.283662,adj[9,5]-0.653644,adj[9,6]-0.989993,adj[9,7]-0.416147,adj[9,8]0.540302,adj[9,9]1.,adj[9,10]0.540302,adj[10,1]-0.91113,adj[10,2]-0.1455,adj[10,3]0.753902,adj[10,4]0.96017,adj[10,5]0.283662,adj[10,6]-0.653644,adj[10,7]-0.989993,adj[10,8]-0.416147,adj[10,9]0.540302,adj[10,10]1.}}