

save_set_movie1 = figure('Name', 'Interactive Movie', 'Position', [320 120 700 700]);

M_Interact = moviein(mframes, save_set_movie1);


 

%11111111111111111111111111111	

for j=1:mframes ;			
	xmdata_pos1 = ones(ndata_pos1,1) ;			
	xmplot_pos1 = ones(nplot_pos1, 1);	
	x2data_pos1 = xdata_pos1 - xc_pos1(j) ;		
	x2plot_pos1 = xplot_pos1 - xc_pos1(j);
	if p_pos1 > 0 ;		
        	for jj = 1:p_pos1 ;		
        	   xmdata_pos1 = [xmdata_pos1,(x2data_pos1.^jj)];
		   xmplot_pos1 = [xmplot_pos1,(x2plot_pos1.^jj)];
        	   jj = jj + 1;		
       		end ;
	end ;
		if can_pos1 == 0;
		if ikern_pos1 == 1 ;			
			kerndata_pos1 = (epan(x2data_pos1./(band_pos1)));
			kernplot_pos1 = (epan(x2plot_pos1./(band_pos1)));		
		elseif ikern_pos1 == 0 ;	
			kerndata_pos1 =sqrt(2*pi)*pdfn(x2data_pos1./(band_pos1)) ;
			kernplot_pos1 = sqrt(2*pi)*pdfn(x2plot_pos1./(band_pos1));
		elseif ikern_pos1 == 2;
			kerndata_pos1 =  uniform(x2data_pos1./band_pos1);
			kernplot_pos1 = uniform(x2plot_pos1./band_pos1);	
		elseif ikern_pos1 == 3;
			kerndata_pos1 =  biwt(x2data_pos1./band_pos1);
			kernplot_pos1 = biwt(x2plot_pos1./band_pos1);
		elseif ikern_pos1 == 4;
			kerndata_pos1 = triwt(x2data_pos1./band_pos1);
			kernplot_pos1 = triwt(x2plot_pos1./band_pos1);

		end ;
	else;
		if ikern_pos1 == 1 ;			
			delta_pos1 = 15^.2			
			kerndata_pos1 = (epan(x2data_pos1./(delta_pos1*band_pos1)));
			kernplot_pos1 = (epan(x2plot_pos1./(delta_pos1*band_pos1)));		
		elseif ikern_pos1 == 0 ;	
			delta_pos1 = (1/(4*pi))^.1;
			kerndata_pos1 = sqrt(2*pi)*pdfn(x2data_pos1./(delta_pos1*band_pos1)) ;
			kernplot_pos1 = sqrt(2*pi)*pdfn(x2plot_pos1./(delta_pos1*band_pos1));
		elseif ikern_pos1 == 2;
			delta_pos1 = (9/2)^.2;
			kerndata_pos1 = uniform(x2data_pos1./(delta_pos1*band_pos1));
			kernplot_pos1 = uniform(x2plot_pos1./(delta_pos1*band_pos1));	
		elseif ikern_pos1 == 3;
			delta_pos1 = 35^.2
			kerndata_pos1 =  biwt(x2data_pos1./(delta_pos1*band_pos1));
			kernplot_pos1 = biwt(x2plot_pos1./(delta_pos1*band_pos1));
		elseif ikern_pos1 == 4;
			delta_pos1 = (9450/143)^.2
			kerndata_pos1 = triwt(x2data_pos1./(delta_pos1*band_pos1));
			kernplot_pos1 = triwt(x2plot_pos1./(delta_pos1*band_pos1));
		end;
	end;

	xwdata_pos1 = xmdata_pos1.*(kerndata_pos1*ones(1,p_pos1+1));	
	hat2_pos1 = inv(xwdata_pos1'*xmdata_pos1)*xwdata_pos1';		
	hatplot_pos1 = xmplot_pos1*hat2_pos1 ;
	hat2_pos1= hat2_pos1(1,:) ;
	yhatplot_pos1 = hatplot_pos1*ydata_pos1 ;
	x_snake1_pos1 = xplot_pos1(abs(x2plot_pos1)<band_pos1) ;
	y_snake1_pos1 = yhatplot_pos1(abs(x2plot_pos1)<band_pos1) ;
	mhat_pos1(j) = hat2_pos1*ydata_pos1 ;

	subplot(2,2,1);

	if color_pos1 ==0			
		for k = 1:ndata_pos1;
			if kerndata_pos1 (k) < .0001;
				a_pos1 = 0;
				b_pos1 = 0;
				else;
				a_pos1 = kerndata_pos1(k);
				b_pos1 = kerndata_pos1(k);
			end;		
			col_pos1 = [a_pos1.^2/2, b_pos1.^2, 1];
			q_pos1 =  plot (xdata_pos1(k), ydata_pos1(k), '+');
			set (q_pos1, 'Color', col_pos1);
			k = k+1;
			hold on;
		end;	
	elseif color_pos1 == 1;
		for k = 1:ndata_pos1;
			if kerndata_pos1 (k) < .0001;
				c_pos1 = 6;
				else;
				c_pos1 = ((kerndata_pos1(k).^2).*15)+6;
			end;
			q_pos1 = plot (xdata_pos1(k), ydata_pos1(k), 'b+');
			set (q_pos1, 'MarkerSize', c_pos1);
			k = k+1;
			hold on
		end;
	
	elseif color_pos1 == 2;
		for k = 1:ndata_pos1;	
			if kerndata_pos1 (k)< .0001;
				a_pos1 = 0;
				b_pos1 = 0;
				c_pos1 = 6;				
				else;
				a_pos1 = kerndata_pos1(k);
				b_pos1 = kerndata_pos1(k);
				c_pos1 = ((kerndata_pos1(k).^2).*15)+6;			
			end;		
			col_pos1 = [a_pos1^2/2, b_pos1^2, 1];
			q_pos1 =  plot (xdata_pos1(k), ydata_pos1(k), '+');
			set (q_pos1, 'Color', col_pos1);	
			set (q_pos1, 'MarkerSize', c_pos1);	
			k =k+1;
			hold on;
		end;
		else;
		plot (xdata_pos1, ydata_pos1, 'b+');
	end;	
	hold on;

	plot(xplot_pos1,mplot_pos1,'w:',xplot_pos1,kernplot_pos1,'m-',x_snake1_pos1,y_snake1_pos1,'y-', xdata_pos1, ydata_pos1, 'b+') ;
	title(pos1_label);	
	hold on;	
	a_pos1 = plot(xc_pos1(1:j), mhat_pos1(1:j), '-');	
	set (a_pos1, 'Color', [.8 .8 0]);
	hold off;	
	set(gca,'Ylim',[0,2.5]);
	set(gca,'Xlim',[-2,2]);	
	h = line([xc(j),xc(j)],[0,2.5]);
	set(h, 'Color', [0 1 1]);
		

%2222222222222222222222222222222222

				
	xmdata_pos2 = ones(ndata_pos2,1) ;			
	xmplot_pos2 = ones(nplot_pos2, 1);	
	x2data_pos2 = xdata_pos2 - xc_pos2(j) ;		
	x2plot_pos2 = xplot_pos2 - xc_pos2(j);
	if p_pos2 > 0 ;		
        	for jj = 1:p_pos2 ;		
        	   xmdata_pos2 = [xmdata_pos2,(x2data_pos2.^jj)];
		   xmplot_pos2 = [xmplot_pos2,(x2plot_pos2.^jj)];
        	   jj = jj + 1;		
       		end ;
	end ;
		if can_pos2 == 0;
		if ikern_pos2 == 1 ;			
			kerndata_pos2 = (epan(x2data_pos2./(band_pos2)));
			kernplot_pos2 = (epan(x2plot_pos2./(band_pos2)));		
		elseif ikern_pos2 == 0 ;	
			kerndata_pos2 =sqrt(2*pi)*pdfn(x2data_pos2./(band_pos2)) ;
			kernplot_pos2 = sqrt(2*pi)*pdfn(x2plot_pos2./(band_pos2));
		elseif ikern_pos2 == 2;
			kerndata_pos2 =  uniform(x2data_pos2./band_pos2);
			kernplot_pos2 = uniform(x2plot_pos2./band_pos2);	
		elseif ikern_pos2 == 3;
			kerndata_pos2 =  biwt(x2data_pos2./band_pos2);
			kernplot_pos2 = biwt(x2plot_pos2./band_pos2);
		elseif ikern_pos2 == 4;
			kerndata_pos2 = triwt(x2data_pos2./band_pos2);
			kernplot_pos2 = triwt(x2plot_pos2./band_pos2);

		end ;
	else;
		if ikern_pos2 == 1 ;			
			delta_pos2 = 15^.2			
			kerndata_pos2 = (epan(x2data_pos2./(delta_pos2*band_pos2)));
			kernplot_pos2 = (epan(x2plot_pos2./(delta_pos2*band_pos2)));		
		elseif ikern_pos2 == 0 ;	
			delta_pos2 = (1/(4*pi))^.1;
			kerndata_pos2 = sqrt(2*pi)*pdfn(x2data_pos2./(delta_pos2*band_pos2)) ;
			kernplot_pos2 = sqrt(2*pi)*pdfn(x2plot_pos2./(delta_pos2*band_pos2));
		elseif ikern_pos2 == 2;
			delta_pos2 = (9/2)^.2;
			kerndata_pos2 = uniform(x2data_pos2./(delta_pos2*band_pos2));
			kernplot_pos2 = uniform(x2plot_pos2./(delta_pos2*band_pos2));	
		elseif ikern_pos2 == 3;
			delta_pos2 = 35^.2
			kerndata_pos2 =  biwt(x2data_pos2./(delta_pos2*band_pos2));
			kernplot_pos2 = biwt(x2plot_pos2./(delta_pos2*band_pos2));
		elseif ikern_pos2 == 4;
			delta_pos2 = (9450/143)^.2
			kerndata_pos2 = triwt(x2data_pos2./(delta_pos2*band_pos2));
			kernplot_pos2 = triwt(x2plot_pos2./(delta_pos2*band_pos2));
		end;
	end;


	
	xwdata_pos2 = xmdata_pos2.*(kerndata_pos2*ones(1,p_pos2+1));	
	hat2_pos2 = inv(xwdata_pos2'*xmdata_pos2)*xwdata_pos2';		
	hatplot_pos2 = xmplot_pos2*hat2_pos2 ;
	hat2_pos2= hat2_pos2(1,:) ;
	yhatplot_pos2 = hatplot_pos2*ydata_pos2 ;
	x_snake2_pos2 = xplot_pos2(abs(x2plot_pos2)<band_pos2) ;
	y_snake2_pos2 = yhatplot_pos2(abs(x2plot_pos2)<band_pos2) ;
	mhat_pos2(j) = hat2_pos2*ydata_pos2 ;
	
	subplot(2,2,2);
	
	if color_pos2 ==0			
		for k = 1:ndata_pos2;
			if kerndata_pos2 (k) < .0001;
				a_pos2 = 0;
				b_pos2 = 0;
				else;
				a_pos2 = kerndata_pos2(k);
				b_pos2 = kerndata_pos2(k);
			end;		
			col_pos2 = [a_pos2.^2/2, b_pos2.^2, 1];
			q_pos2 =  plot (xdata_pos2(k), ydata_pos2(k), '+');
			set (q_pos2, 'Color', col_pos2);
			k = k+1;
			hold on;
		end;	
	elseif color_pos2 == 1;
		for k = 1:ndata_pos2;
			if kerndata_pos2 (k) < .0001;
				c_pos2 = 6;
				else;
				c_pos2 = ((kerndata_pos2(k).^2).*15)+6;
			end;
			q_pos2 = plot (xdata_pos2(k), ydata_pos2(k), 'b+');
			set (q_pos2, 'MarkerSize', c_pos2);
			k = k+1;
			hold on
		end;
	
	elseif color_pos2 == 2;
		for k = 1:ndata_pos2;	
			if kerndata_pos2 (k)< .0001;
				a_pos2 = 0;
				b_pos2 = 0;
				c_pos2 = 6;				
				else;
				a_pos2 = kerndata_pos2(k);
				b_pos2 = kerndata_pos2(k);
				c_pos2 = ((kerndata_pos2(k).^2).*15)+6;			
			end;		
			col_pos2 = [a_pos2^2/2, b_pos2^2, 1];
			q_pos2 =  plot (xdata_pos2(k), ydata_pos2(k), '+');
			set (q_pos2, 'Color', col_pos2);	
			set (q_pos2, 'MarkerSize', c_pos2);	
			k =k+1;
			hold on;
		end;
	else;
	plot (xdata_pos2, ydata_pos2, 'b+');
	end;	
	hold on;
	plot(xplot_pos2,mplot_pos2,'w:',xplot_pos2,kernplot_pos2,'m-',x_snake2_pos2,y_snake2_pos2,'g-', xdata_pos2, ydata_pos2, 'b+') ;
	title(pos3_label);	
	hold on;	
	a_pos2 = plot(xc_pos2(1:j), mhat_pos2(1:j), '-');	
	set (a_pos2, 'Color', [.2 .6 0]);
	hold off;	
	set(gca,'Ylim',[0,2.5]);
	set(gca,'Xlim',[-2,2]);	
	h = line([xc(j),xc(j)],[0,2.5]);
	set(h, 'Color', [0 1 1]);


%3333333333333333333333333333333333333
	
	xmdata_pos3 = ones(ndata_pos3,1) ;			
	xmplot_pos3 = ones(nplot_pos3, 1);	
	x2data_pos3 = xdata_pos3 - xc_pos3(j) ;		
	x2plot_pos3 = xplot_pos3 - xc_pos3(j);
	if p_pos3 > 0 ;		
        	for jj = 1:p_pos3 ;		
        	   xmdata_pos3 = [xmdata_pos3,(x2data_pos3.^jj)];
		   xmplot_pos3 = [xmplot_pos3,(x2plot_pos3.^jj)];
        	   jj = jj + 1;		
       		end ;
	end ;
		if can_pos3 == 0;
		if ikern_pos3 == 1 ;			
			kerndata_pos3 = (epan(x2data_pos3./(band_pos3)));
			kernplot_pos3 = (epan(x2plot_pos3./(band_pos3)));		
		elseif ikern_pos3 == 0 ;	
			kerndata_pos3 =sqrt(2*pi)*pdfn(x2data_pos3./(band_pos3)) ;
			kernplot_pos3 = sqrt(2*pi)*pdfn(x2plot_pos3./(band_pos3));
		elseif ikern_pos3 == 2;
			kerndata_pos3 =  uniform(x2data_pos3./band_pos3);
			kernplot_pos3 = uniform(x2plot_pos3./band_pos3);	
		elseif ikern_pos3 == 3;
			kerndata_pos3 =  biwt(x2data_pos3./band_pos3);
			kernplot_pos3 = biwt(x2plot_pos3./band_pos3);
		elseif ikern_pos3 == 4;
			kerndata_pos3 = triwt(x2data_pos3./band_pos3);
			kernplot_pos3 = triwt(x2plot_pos3./band_pos3);

		end ;
	else;
		if ikern_pos3 == 1 ;			
			delta_pos3 = 15^.2			
			kerndata_pos3 = (epan(x2data_pos3./(delta_pos3*band_pos3)));
			kernplot_pos3 = (epan(x2plot_pos3./(delta_pos3*band_pos3)));		
		elseif ikern_pos3 == 0 ;	
			delta_pos3 = (1/(4*pi))^.1;
			kerndata_pos3 = sqrt(2*pi)*pdfn(x2data_pos3./(delta_pos3*band_pos3)) ;
			kernplot_pos3 = sqrt(2*pi)*pdfn(x2plot_pos3./(delta_pos3*band_pos3));
		elseif ikern_pos3 == 2;
			delta_pos3 = (9/2)^.2;
			kerndata_pos3 = uniform(x2data_pos3./(delta_pos3*band_pos3));
			kernplot_pos3 = uniform(x2plot_pos3./(delta_pos3*band_pos3));	
		elseif ikern_pos3 == 3;
			delta_pos3 = 35^.2
			kerndata_pos3 =  biwt(x2data_pos3./(delta_pos3*band_pos3));
			kernplot_pos3 = biwt(x2plot_pos3./(delta_pos3*band_pos3));
		elseif ikern_pos3 == 4;
			delta_pos3 = (9450/143)^.2
			kerndata_pos3 = triwt(x2data_pos3./(delta_pos3*band_pos3));
			kernplot_pos3 = triwt(x2plot_pos3./(delta_pos3*band_pos3));
		end;
	end;

	
	xwdata_pos3 = xmdata_pos3.*(kerndata_pos3*ones(1,p_pos3+1));	
	hat2_pos3 = inv(xwdata_pos3'*xmdata_pos3)*xwdata_pos3';		
	hatplot_pos3 = xmplot_pos3*hat2_pos3 ;
	hat2_pos3= hat2_pos3(1,:) ;
	yhatplot_pos3 = hatplot_pos3*ydata_pos3 ;
	x_snake3_pos3 = xplot_pos3(abs(x2plot_pos3)<band_pos3) ;
	y_snake3_pos3 = yhatplot_pos3(abs(x2plot_pos3)<band_pos3) ;
	mhat_pos3(j) = hat2_pos3*ydata_pos3 ;
	
	subplot(2,2,3);
	
	if color_pos3 ==0			
		for k = 1:ndata_pos3;
			if kerndata_pos3 (k) < .0001;
				a_pos3 = 0;
				b_pos3 = 0;
				else;
				a_pos3 = kerndata_pos3(k);
				b_pos3 = kerndata_pos3(k);
			end;		
			col_pos3 = [a_pos3.^2/2, b_pos3.^2, 1];
			q_pos3 =  plot (xdata_pos3(k), ydata_pos3(k), '+');
			set (q_pos3, 'Color', col_pos3);
			k = k+1;
			hold on;
		end;	
	elseif color_pos3 == 1;
		for k = 1:ndata_pos3;
			if kerndata_pos3 (k) < .0001;
				c_pos3 = 6;
				else;
				c_pos3 = ((kerndata_pos3(k).^2).*15)+6;
			end;
			q_pos3 = plot (xdata_pos3(k), ydata_pos3(k), 'b+');
			set (q_pos3, 'MarkerSize', c_pos3);
			k = k+1;
			hold on
		end;
	
	elseif color_pos3 == 2;
		for k = 1:ndata_pos3;	
			if kerndata_pos3 (k)< .0001;
				a_pos3 = 0;
				b_pos3 = 0;
				c_pos3 = 6;				
				else;
				a_pos3 = kerndata_pos3(k);
				b_pos3 = kerndata_pos3(k);
				c_pos3 = ((kerndata_pos3(k).^2).*15)+6;			
			end;		
			col_pos3 = [a_pos3^2/2, b_pos3^2, 1];
			q_pos3 =  plot (xdata_pos3(k), ydata_pos3(k), '+');
			set (q_pos3, 'Color', col_pos3);	
			set (q_pos3, 'MarkerSize', c_pos3);	
			k =k+1;
			hold on;
		end;
		else;
		plot (xdata_pos3, ydata_pos3, 'b+');
	end;	
	hold on;
	plot(xplot_pos3,mplot_pos3,'w:',xplot_pos3,kernplot_pos3,'m-',x_snake3_pos3,y_snake3_pos3,'r-', xdata_pos3, ydata_pos3, 'b+') ;
	title(pos3_label);	
	hold on;	
	a_pos3 = plot(xc_pos3(1:j), mhat_pos3(1:j), '-');	
	set (a_pos3, 'Color', [.6 0 0]);
	hold off;	
	set(gca,'Ylim',[0,2.5]);
	set(gca,'Xlim',[-2,2]);	
	h = line([xc(j),xc(j)],[0,2.5]);
	set(h, 'Color', [0 1 1]);

%4444444444444444444444444444444444
	subplot(2,2,4);
	plot(xplot_pos1,mplot_pos1, 'w:', xplot_pos2,mplot_pos2, 'w:', xplot_pos3, mplot_pos3, 'w:',x_snake1_pos1, y_snake1_pos1, 'y-', x_snake2_pos2, y_snake2_pos2,'g-', x_snake3_pos3, y_snake3_pos3, 'r-');
	title ('Comparison Window');	
	hold on;
	a1 = plot(xc(1:j), mhat_pos1(1:j), '-');	
	set (a1, 'Color', [.8 .8 0]);
	b2 = plot(xc(1:j), mhat_pos2(1:j), '-');
	set (b2, 'Color', [.2 .6 0]);
	c3 = plot(xc(1:j), mhat_pos3(1:j), '-');
	set (c3, 'Color', [.6 0 0]);
	hold off;


	set(gca, 'Ylim', [0, 2.5]);
	set(gca, 'Xlim', [-2,2]);
	h = line([xc(j),xc(j)],[0,2.5]);
	set(h, 'Color',[0 1 1]);
figure(save_set_movie1);	
M_Interact(:,j) = getframe(save_set_movie1);

end ;

       		
clf;
h = axes('Visible', 'off');
h = axes('Position', [0 0 1 1]);
movie(M_Interact, 6);

