CONFIG

USERHOMOTOPY: 1;
TRACKTOLBEFOREEG: 1e-3;
TRACKTOLDURINGEG: 1e-4;
FINALTOL: 1e-5;
IMAGTHRESHOLD: 1e-4;
END;

INPUT

variable u1, u2, lambda1, lambda2, lambda3, lambda4, mu1, mu2, mu3, mu4, mu5, mu6, iL2, iL3, vo2, vo3;
constant iL11, vo11, ro1, uold1, iLmax1, vref1;
constant iL12, vo12, ro2, uold2, iLmax2, vref2;
parameter iL1, vo1, ro, uold, iLmax, vref;
pathvariable t;
function f1, f2, f3, f4, f5, f6, f7, f8, f9, f10, f11, f12, f13, f14, f15, f16;

iL11 =  -5.7623+(28.763)*I;
vo11 =  -333.88+(51.269)*I;
ro1 =  7.88+(-5.269)*I;
uold1 = 8.31+(-31.269)*I;
iLmax1= 3.4+(1.269)*I;
vref1 = 33.88+(51.269)*I;

iL12 =  0.5;
vo12 =  0.5;
ro2 =  10;
uold2 = 0.4;
iLmax2= 2;
vref2 = 1;

iL1 = t*iL11 + (1-t)*iL12;
vo1 = t*vo11 + (1-t)*vo12;
ro = t*ro1 + (1-t)*ro2;
uold = t*uold1 + (1-t)*uold2;
iLmax = t*iLmax1 + (1-t)*iLmax2;
vref = t*vref1 + (1-t)*vref2;

f1=0.4*u1-0.2*uold-0.2*u2+lambda1+mu1*(2-2*vo1)+mu2-mu3;
f2=-0.5*lambda1+3.52*lambda2*ro+0.48*lambda3-2.52*lambda4*ro+mu4;
f3=20*vo2-20*vref+lambda2*(-352*ro-3.52)-lambda3+lambda4*(352*ro+2.52)-2*mu4*u2;
f4=0.2*u2-0.2*u1+lambda3+mu4*(2-2*vo2)+mu5-mu6;
f5=-0.5*lambda3+3.52*lambda4*ro;
f6=20*vo3-20*vref+lambda4*(-352*ro-3.52);
f7=0.48*iL1-vo1+u1-0.5*iL2;
f8=-2.52*iL1*ro-352*ro*vo2+352*ro*vo1+3.52*ro*iL2+2.52*vo1-3.52*vo2;
f9=0.48*iL2-vo2+u2-0.5*iL3;
f10=-2.52*ro*iL2-352*ro*vo3+352*ro*vo2+3.52*ro*iL3+2.52*vo2-3.52*vo3;
f11=mu1*(iL1-iLmax+2*u1-2*vo1*u1);
f12=mu2*(u1-1);
f13=-mu3*u1;
f14=mu4*(iL2-iLmax+2*u2-2*vo2*u2);
f15=mu5*(u2-1);
f16=-mu6*u2;

END;