function [ PHISOL ] = calculgramraideur(NPT,dt,z0M,dz0M,alpha0M,dalpha0M,xeqM,Rad,B0,F0,a0)
%==================
% construction de la matrice de Gram
%==================
GRAMM=zeros(4,4);
%==================
phi0=0.*linspace(1,4,4);
deltaphi0=0.*linspace(1,4,4);
Lsec=linspace(1,2,4)*0.;
%=================
for k1=1:4
  for i1=1:4
        if(i1 == k1)
                phi0(i1)=1;
        else
                phi0(i1)=0;
        end
  end
    [P1,P2]=soluedo(Rad,NPT,dt,phi0);
    P1dotT=(P1(NPT)-P1(NPT-1))/dt;
    P2dotT=(P2(NPT)-P2(NPT-1))/dt;
     Lsec(k1)=dz0M*phi0(1)+dalpha0M*phi0(3)-z0M*phi0(2)-alpha0M*phi0(4)+xeqM(1)*P1dotT+xeqM(2)*P2dotT;
     squatro=0;
     for jkl=1:NPT
         squatro=squatro+F0(1)*P1(jkl)+F0(2)*P2(jkl);
     end
     Lsec(k1)=Lsec(k1)+squatro*dt;
  %==============
   for k2=k1:4
       for i2=1:4
          if(i2 == k2)
               deltaphi0(i2)=1;
          else
               deltaphi0(i2)=0;
          end
        end
%==================
% Remplissage des solution elementaire
%==================
    [Q1,Q2]=soluedo(Rad,NPT,dt,deltaphi0);
%========   
% calcul des termes de GRAMM(Monsieur Gram)
%========
        sbis=0;
        %
        for jj=1:NPT
        sbis=sbis+(B0(1)*P1(jj)+B0(2)*P2(jj))*(B0(1)*Q1(jj)+B0(2)*Q2(jj));
        end
        
  
    GRAMM(k1,k2)=sbis*dt;
    GRAMM(k2,k1)=GRAMM(k1,k2);
   end
end
GRAMM=GRAMM/a0;
%======
%  calcul de PHISOL
%======
PHISOL=linsolve(GRAMM,Lsec');
end




