diff --git a/GUI/descriptive_part.m b/GUI/descriptive_part.m
index 98e0b4c..97bcee6 100644
--- a/GUI/descriptive_part.m
+++ b/GUI/descriptive_part.m
@@ -80,9 +80,9 @@ function descriptive_part(varargin)
     histfit(MeasureM.meanPathLength,nbins,'gamma');
     title('Path Length(MEAN)')
     if isunix
-        print(strcat('Descriptive_STAT/Eye_Mvt_Measure_Dist'),'-dpng')
+        print(h1,strcat('Descriptive_STAT/Eye_Mvt_Measure_Dist'),'-dpng')
     else
-        print(strcat('Descriptive_STAT\Eye_Mvt_Measure_Dist'),'-dpng')
+        print(h1,strcat('Descriptive_STAT\Eye_Mvt_Measure_Dist'),'-dpng')
     end
     
     % Mean of Fixation Map
@@ -93,9 +93,9 @@ function descriptive_part(varargin)
     imagesc(FixationMean);
     axis off,axis equal
     if isunix
-        print(strcat('Descriptive_STAT/Mean_Intensity_Map'),'-dpng')
+        print(h2,strcat('Descriptive_STAT/Mean_Intensity_Map'),'-dpng')
     else
-        print(strcat('Descriptive_STAT\Mean_Intensity_Map'),'-dpng')
+        print(h2,strcat('Descriptive_STAT\Mean_Intensity_Map'),'-dpng')
     end
     close(h1,h2)
     
@@ -136,9 +136,9 @@ function descriptive_part(varargin)
         boxplot(MeasureM.meanPathLength,Conditiontmp)
         title('Path Length(MEAN)')
         if isunix
-            print(strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Boxplot'),'-dpng')
+            print(h1,strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Boxplot'),'-dpng')
         else
-            print(strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Boxplot'),'-dpng')
+            print(h1,strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Boxplot'),'-dpng')
         end
         close(h1)
         
@@ -181,9 +181,9 @@ function descriptive_part(varargin)
                 
                 if pp==1
                     if isunix
-                        print(strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
+                        print(h(k-1),strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
                     else
-                        print(strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
+                        print(h(k-1),strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
                     end
                     h(k)=figure('NumberTitle','off','Name',['Mean_FixaionIntensity_Map - ' CName{ic} ' continue'],...
                         'position',[xpos,ypos,screensize(3),screensize(4)]);
@@ -192,9 +192,9 @@ function descriptive_part(varargin)
                 end
             end
             if isunix
-                print(strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map',num2str(k)),'-dpng')
+                print(h(k),strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map',num2str(k)),'-dpng')
             else
-                print(strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map',num2str(k)),'-dpng')
+                print(h(k),strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map',num2str(k)),'-dpng')
             end
             close(gcf)
         else
@@ -212,9 +212,9 @@ function descriptive_part(varargin)
             end
             % save figures
             if isunix
-                print(strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map'),'-dpng')
+                print(h(1),strcat('Descriptive_STAT/',CName{ic},'/',CName{ic},'_Mean_Intensity_Map'),'-dpng')
             else
-                print(strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map'),'-dpng')
+                print(h(1),strcat('Descriptive_STAT\',CName{ic},'\',CName{ic},'_Mean_Intensity_Map'),'-dpng')
             end
             %close figure
             close(gcf);
@@ -282,9 +282,9 @@ function descriptive_part(varargin)
                     
                     if pp==1
                         if isunix
-                            print(strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
+                            print(h(k-1),strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
                         else
-                            print(strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
+                            print(h(k-1),strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map',num2str(k-1)),'-dpng')
                         end
                         h(k)=figure('NumberTitle','off','Name',['Mean Fixation Intensity - JointCondion continue'],...
                             'position',[xpos,ypos,screensize(3),screensize(4)]);
@@ -293,9 +293,9 @@ function descriptive_part(varargin)
                     end
                 end
                 if isunix
-                    print(strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map',num2str(k)),'-dpng')
+                    print(h(k),strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map',num2str(k)),'-dpng')
                 else
-                    print(strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map',num2str(k)),'-dpng')
+                    print(h(k),strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map',num2str(k)),'-dpng')
                 end
                 close(gcf)
             else
@@ -308,9 +308,9 @@ function descriptive_part(varargin)
                     axis off,axis equal
                 end
                 if isunix
-                    print(strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map'),'-dpng')
+                    print(h(1),strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Mean_Intensity_Map'),'-dpng')
                 else
-                    print(strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map'),'-dpng')
+                    print(h(1),strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Mean_Intensity_Map'),'-dpng')
                 end
                 close(gcf);
             end
@@ -320,7 +320,7 @@ function descriptive_part(varargin)
             % %hfigs(hfigs ==handles.figure1)=[];
             % close(hfigs(1:end-1)) % keep the figure of the linear model open
             
-            %% box plot for  joint condition
+            %% box plot for joint condition
             h(1) = figure('NumberTitle','off','Name','Eye movement measurement Boxplot',...
                 'position',[xpos,ypos,screensize(3)/2,screensize(4)]);
             
@@ -344,9 +344,9 @@ function descriptive_part(varargin)
             boxplot(MeasureM.meanPathLength,conditiontmp);
             title('Path Length(MEAN)')
             if isunix
-                print(strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Boxplot'),'-dpng');
+                print(h(1),strcat('Descriptive_STAT/',CNameall,'/',CNameall,'_Boxplot'),'-dpng');
             else
-                print(strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Boxplot'),'-dpng');
+                print(h(1),strcat('Descriptive_STAT\',CNameall,'\',CNameall,'_Boxplot'),'-dpng');
             end
             close(gcf)
             
diff --git a/Matlab_Installation_Package/iMap4.mlappinstall b/Matlab_Installation_Package/iMap4.mlappinstall
index 53c6471..2b41484 100644
Binary files a/Matlab_Installation_Package/iMap4.mlappinstall and b/Matlab_Installation_Package/iMap4.mlappinstall differ