function [SN,FN]=faceNouvelle(S,F,A,k)
% [SN,FN]=FACENOUVELLE(S,F,A,k) coupe un polyèdre régulier donné par S(ommets et F(aces
% K est le coefficient de découpage des arêtes et crée les nouvelles faces.

[nbs,co]=size(S);		% nombre de sommets des faces 
[li,co]=size(F);		% taille de F
[i,j]=find(F(:,2:co-1)==1);	% recherche du nbre d'occurrences du sommet 1
							% dans F (= nombre d'arêtes en un sommet)
pole=length(i);
SN=zeros(nbs*pole,3); 	% matrice des sommets
FN=zeros(nbs,pole+2);	% matrice des faces
AA=ones(pole,1);			% reformattage de matrice
[nba,co]=size(A);
for i=1:nbs
	d=pole*(i-1);
	[j,m]=find(A(:,1)==i);			% recherche des indices des sommets 
	ind=A(j,2);
	SN(d+1:d+pole,:)=AA*S(i,:)+k*(S(ind,:)-AA*S(i,:));
	FN(i,:)=[pole+1 d+1:d+pole d+1];
end

