out qq;
off nat;

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

FI2 := (E**(O(1)*S)*G(3)*G(2)*O(3) - E**(O(1)*S)*G(3)*G(2)*O(2) + E**
(O(2)*S)*G(3)*G(1)*O(3) - E**(O(2)*S)*G(3)*G(1)*O(1) + E**(O(
3)*S)*G(2)*G(1)*O(2) - E**(O(3)*S)*G(2)*G(1)*O(1))/(E**
(O(3)*S + O(2)*S + O(1)*S)*(O(3)**2*O(2) - O(3)**2*O(1) - O(3)*O(2
)**2 + O(3)*O(1)**2 + O(2)**2*O(1) - O(2)*O(1)**2))$


%========= ok solutions ======================================

q(1):= -log(fi1);
q(2):=  log(fi1) - log(fi2) ;
q(3):= +log(fi2);

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

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

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


O(3):=-(O(1)+O(2));
g(3):=1/(g(1)*g(2));


eqq(1);
eqq(2);
eqq(3);

%%@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
%@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

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

FI2 := (E**(O(1)*S)*G(3)*G(2)*O(3) - E**(O(1)*S)*G(3)*G(2)*O(2) + E**
(O(2)*S)*G(3)*G(1)*O(3) - E**(O(2)*S)*G(3)*G(1)*O(1) + E**(O(
3)*S)*G(2)*G(1)*O(2) - E**(O(3)*S)*G(2)*G(1)*O(1))/(E**
(O(3)*S + O(2)*S + O(1)*S)*(O(3)**2*O(2) - O(3)**2*O(1) - O(3)*O(2
)**2 + O(3)*O(1)**2 + O(2)**2*O(1) - O(2)*O(1)**2))$

q(1):= -log(fi1);
q(2):=  log(fi1) - log(fi2) ;
q(3):= +log(fi2);

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

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

O(1):=l1*O(1);
O(2):=l1*O(2);
O(3):=l1*O(3);

g1:=l1**2*exp(O1*(l2+2l1));
g2:=l1**2*exp(O2*(l2+2l1));
g3:=l1**2*exp(O3*(l2+2l1));


%===========================================
O(3):=-(O(1)+O(2));
g(3):=1/(g(1)*g(2));

eqq(1);
eqq(2);
eqq(3);



shut qq;
bye;
end;