out uu;
off nat;

%FI1 := ( E**(-O(1)*(exp(-l1/2)*S-l2))*G(1)+ E**(-O(2)*(exp(-l1/2)*S-l2) )*G(2) )
%/( O(2)-O(1))$

%FI1 := ( E**(-O(1)*l1*S)*G(1)+ E**(-O(2)*l1*S )*G(2) )
%/( O(2)-O(1))$

FI1 := ( E**(-O(1)*(l1*S-l2))*G(1)+ E**(-O(2)*(l1*S-l2) )*G(2) )
/( O(2)-O(1))$



%------------------------------------

q(1):= -log(fi1)+ l4-log(l1**2) -l3*(exp(-l1/2)*S-l2);
q(2):=  log(fi1)+ l4- 2*log(l1**2)-l3*(exp(-l1/2)*S-l2);

%q(1):= -log(fi1);
%q(2):=  log(fi1);


 

%----------------------------------

eqq(1):= df(q(1),S,2) + exp(q(1) -q(2))  ;

eqq(2):=df(q(2),S,2) - exp(q(1) -q(2))   ;

%eqq(1):= df(q(1),S,2) +2*exp(q(1)) -2*exp(q(2))  ;
%eqq(2):=df(q(2),S,2) - 2*exp(q(1)) +2*exp(q(2))   ;






%--------------------------------



O(2) :=  -  O(1)$

g(1):=1/g(2);



eqq(1);

eqq(2);

shut uu;
bye;
end;