function [SE,FE]=polyedreEtoile(S,F,k)
% [SE,FE]=POLYEDREETOILE(S,F,K) crée un polyèdre étoilé à partir
% du polyèdre donné par S et F. K est le rapport entre le sommet
% ajouté pour le rendre étoilé et la sphère.

[nbf,co]=size(F);		% nombre de faces
[nbs,co]=size(S);		% nombre de sommets

SE=zeros(nbs+nbf,3);	% matrice des sommets
SE(1:nbs,:)=S;			% les premiers sommets sont les mêmes
nbsf=F(1,1)-1;			% nombre de sommets par face
FE=zeros(nbsf*nbf,5);	% matrice des faces

for i=1:nbf
	d=nbsf*(i-1);
	ind=F(i,2:nbsf+1);
	p=[S(ind,1) S(ind,2) S(ind,3)];
	g=barycentre(p);
	SE(i+nbs,:)=k*g;
	ind=F(i,2:nbsf+2);
	m=1;
	for j=d+1:d+nbsf
		FE(j,:)=[4 i+nbs ind(m) ind(m+1) i+nbs];
		m=m+1;
	end
end
