%
% Script de test pour le dernier programme 2011-12
%


function matlab8()
close all;

alpha = 4;
n = 4*alpha+1;
h = 2/(n-1);

tic
[U,Umax,Xmax] = theBigFish(alpha);
toc


fprintf('\n Maximum = %14.7e at (%14.7e,%14.7e,%4.7e) ',Umax,Xmax);
Xi = -1:h:1;
Yi = -1:h:1;
Zi = -1:h:1;

figure
subplot(2,2,1); contourf(Xi,Zi,reshape(U(:,3*alpha+1,:),n,n)); caxis([0,Umax]); axis off; axis equal;
subplot(2,2,2); contourf(Xi,Yi,U(:,:,alpha+1)); caxis([0,Umax]); axis off; axis equal;
subplot(2,2,3); contourf(Xi,Yi,U(:,:,2*alpha+1)); caxis([0,Umax]); axis off; axis equal;
subplot(2,2,4); contourf(Xi,Yi,U(:,:,3*alpha+1)); caxis([0,Umax]); axis off; axis equal;


figure
[X Y Z] = meshgrid(Xi,Yi,Zi); 
slice(X,Y,Z,U,[0,0.5,1],[ 1],[0.2 -1],'cubic'); axis off; axis equal;



end