In[]:=
numnodes=2;numfamilies=1;Adj=Table[adj[i,j],{i,1,numnodes},{j,1,numnodes}];varsnodes=Flatten[Adj];Adjfam=Table[famadj[i,j],{i,1,numfamilies},{j,1,numfamilies}];varsfam=Flatten[Adjfam];vars=Flatten[{varsnodes,varsfam}];(*totaladj[args__]*)totaldj[adj__,famadj__]:=((*adj[1,1]=args[[1]];adj[1,2]=args[[2]];adj[2,1]=args[[3]];adj[2,2]=args[[4]];famadj[1,1]=args[[5]];*)For[k=0;sum=0,k<=1,k++,cost[k]=0;imax=Which[k==0,numnodes,k==1,numfamilies];For[i=1,i<=imax,i++,For[j=1,j<=imax,j++,If[Which[k==0,adj[i,j],k==1,famadj[i,j]]>0.1/imax^2,Print["hi"];cost[k]+=1]]];sum+=cost[k]];sum);vars2={1,1,1,1,1};totaladj[vars2]sol=NMinimize[totaladj[vars],vars];
Out[]=
totaladj[{1,1,1,1,1}]
