if ConnectElemPhase2(i)==1
%the point lies within the inclusion
lambda=Prop2(1);
mu=Prop2(2);
elseif ConnectElemPhase2(i)==2
%the point lies within the matrix
lambda=Prop2(3);
mu=Prop2(4);
end
%3 Matrix Phase Field (dans les noeuds d'element i)
%df=[dresult_old2(n1) dresult_old2(n2) dresult_old2(n3) dresult_old2(n4) dresult_old2(n5) dresult_old2(n6) dresult_old2(n7) dresult_old2(n8)];
%LOOP OVER ALL POINTS GAUSS
for kG=1:8
kkk=0;
Assemble=zeros(1,3*length(CubNodes));
for iii=1:length(CubNodes)
kkk=kkk+1;
Assemble(kkk)=CubNodes(iii)*3-2;
kkk=kkk+1;
Assemble(kkk)=CubNodes(iii)*3-1;
kkk=kkk+1;
Assemble(kkk)=CubNodes(iii)*3;
end
Ui=Uresult_old2(Assemble);
%Computer Strain
Epsilon=MatD2(:,:,kG)*Ui;
%Computer Part Positive of strain et Tensor Projection