clear all; close all; clc;

gridesize=11;
% numbers are arbitrary
cube=zeros(11,11,11);
cube(3:9,3:9,3:9)=5; % Create a cube inside the region

% Boring: faces of the cube are a different color.
cube(3:9,3:9,3)=1;
cube(3:9,3:9,9)=1;
cube(3:9,3,3:9)=1;
cube(3:9,9,3:9)=1;
cube(3,3:9,3:9)=2;
cube(9,3:9,3:9)=1;

H=vol3d('Cdata',cube,'alpha',cube/5) %what is alpha?

figure;
[vol_handle]=VoxelPlotter(cube,0.5); %size of each voxel  
view(3);
daspect([1,1,1]);
set(gca,'xlim',[0 gridesize], 'ylim',[0 gridesize], 'zlim',[0 gridesize]);
xlabel('X');ylabel('Y');zlabel('Z');

%figure;
%plot3(cube(:,:,1),cube(:,1,:));

%connects surfaces with specified value

figure;
x = 1:11;
y = 1:11;
z = 1:11;
v = cube;
alpha=isosurface(x,y,z,v,1);
p = patch(isosurface(x,y,z,v,5));
isonormals(x,y,z,v,p)
p.FaceColor = 'red';
p.EdgeColor = 'none';

daspect([1,1,1])
view(3); axis tight
camlight 
lighting gouraud