void obliczanieLFiLK(MACIERZ* _mac){
   size_t i;
   //bo dziwne wyniki wyskakuja dlatego tworze i przypisuje Alfe
   double Alfaa=0.;

   for(i=0;i<liczba_elementow;i++){
      dR=_mac->element[i].wezel2->x - _mac->element[i].wezel1->x;
      size_t j;
      if(i == (liczba_elementow-1))
          Alfaa=Alfa;

      _mac->element[i].K[0][0] = 0.;
      _mac->element[i].K[0][1] = 0.;
      _mac->element[i].K[1][0] = 0.;
      _mac->element[i].K[1][1] = 0.;

      _mac->element[i].F[0] = 0.;
      _mac->element[i].F[1] = 0.;

      for(j=0;j<Np;j++){
          double Rp=N1[j]*_mac->element[i].wezel1->x + N2[j]*_mac->element[i].wezel2->x;
          double TpTau= N1[j]*_mac->element[i].wezel1->temp + N2[j]*_mac->element[i].wezel2->temp;

          _mac->element[i].K[0][0] += K*Rp/dR + C*Ro*dR*Rp*N1[j]*N1[j]/dTau;
          _mac->element[i].K[0][1] += -(K*Rp/dR) + C*Ro*dR*Rp*N1[j]*N2[j]/dTau;
          _mac->element[i].K[1][0] = _mac->element[i].K[0][1];
          _mac->element[i].K[1][1] += K*Rp/dR + C*Ro*dR*Rp*N2[j]*N2[j]/dTau ;

          _mac->element[i].F[0] += -(C*Ro*dR*TpTau*Rp*N1[j]/dTau);
          _mac->element[i].F[1] += -(C*Ro*dR*TpTau*Rp*N2[j]/dTau) - 2*Alfaa*Rmax*TempPow;
      }
	  if(i == (liczba_elementow-1)){
		  _mac->element[i].K[1][1] += 2*Alfaa*Rmax;
	  }
   }
}