%Enter Directory
load('DIRECTORY/input.m')

boxplot(S_no_0,'Labels',{'1964','1965','1966','1967','1968',...
    '1969','1970','1971','1972','1973','1974','1975','1976',...
    '1977','1978','1979','1980','1981','1982','1983','1984',...
    '1985','1986','1987','1988','1989','1990','1991','1992',...
    '1993','1994','1995','1996','1997','1998','1999','2000',...
    '2001','2002','2003','2004','2005','2006','2007','2008','2009','2010'})
set(gca,'XTickLabelRotation',45)

title('The boxplot of the annual sparse values','fontsize',25)
filename='boxplot_S';
set(gcf,'units','normalized','outerposition',[0 0 1 1])

%Enter Directory
fname = 'DIRECTORY';

saveas(gca, fullfile(fname, filename),'tif');
