(*This file contains routines which make use of Nminimize*) (*In order to use these functions, an exepac
must be loaded, and a set of ps values must be defined*) (*The "f" equations and yy21 are used*)

yrules={0<=y1<=1,0<=y2<=1,y1+y2<=1,0<=y3<=1,0<=y4<=1,0<=y5<=1,0<=y6<=1,0<=y7<=1,0<=y8<=1,0<=y9<=1,0<=y10<=1,
0<=y11<=1,0<=y12<=1,0<=y13<=1,0<=y14<=1,0<=y15<=1,0<=y16<=1,0<=y17<=1,0<=y18<=1,0<=y19<=1,0<=y20<=1,0<=y21<=1};

nMin1[rules_]:= Module[{nmw1, nmw2, r, f2},

Clear[f2];
prep[];

f2=f[[1]]+f[[2]]+f[[3]]+f[[4]]+f[[5]]+f[[6]]+f[[7]]+f[[8]]+f[[9]]+f[[10]]+f[[11]]+f[[12]]+f[[13]]+f[[14]]+
f[[15]]+f[[16]]+f[[17]]+f[[18]]+f[[19]]+f[[20]]+f[[21]];

nmw1={f2};
nmw2=Append[nmw1,rules];
r=NMinimize[nmw2,yy21];

Return[r]; ]

end=lastline