%Modifié le 03/02/2004 par Claire
% fonction d'ecriture de equations

function FONC=equat(x,y)
global nc REPAR X Vrup BI1 BI2 ;

TERME=zeros(nc,1);
TERME2=zeros(nc,1);
TERME1=zeros(nc,1);

for i=1:nc
  j=1;
  while REPAR(i,j,1)>0
     delta=1;
     if REPAR(i,j,1)==REPAR(i,j,2)
       delta=0.5;
     end
     TERME2(i)=TERME2(i)+delta*REPAR(i,j,3)*q(X(REPAR(i,j,1)),X(REPAR(i,j,2)))*y(REPAR(i,j,1))*y(REPAR(i,j,2));
     j=j+1;
  end
  for j=1:nc
     TERME(i)=TERME(i)-y(i)*q(X(i),X(j))*y(j);
  end
end
FONC=TERME2+TERME;

VrupY=y.*Vrup; % y matrice de dimension [x,nc]; Vrup vecteur colonne (nc,1)
SOMME= BI1 + BI2; % somme de 2 matrice de meme dimension [nc,nc]
TERME1=SOMME*VrupY; % VrupY vecteur colonne de (nc,1)
%FONC= FONC + TERME1-VrupY; % FONC vecteur colonne (nc,1), doit etre de la mme dimension que vecteur initial N
 
%   for k=i:nc
%      TERME1(i)=TERME1(i)+(BI1(i,k)+BI2(i,k))*Vrup(k)*y(k);
%   end
  %FONC(i)=TERME2(i)+TERME(i)+TERME1(i)-Vrup(i)*y(i);

%FONC=FONC';