% Recall using Lukasiewicz IFAM...

clear;
load faces;

W = minI_P(X,X'); theta=min(X,[],2);

% Para efeitos computacionais!
Theta=[]; 
for i=1:size(Xin,2)
	Theta=[Theta,theta];
end
	
yIFAM_P=max(maxC_P(W,Xin),Theta);
disp('Goguen IFAM test - concluded');

figure(3); clf;
for i=1:size(yIFAM_P,2)
	subplot(4,6,i); imshow(reshape(yIFAM_P(:,i),s1,s2));
	disp(sprintf('%.4f',norm(yIFAM_P(:,i)-X(:,1))/norm(X(:,1))));
end
print -deps facesIFAM_P.eps
