function view_HMM(P,Q,pi_i,gamma,spk_trains,time_vec),

	%Plot raster of response    
	Rasterize_Mult(spk_trains(:,250:4000,:),time_vec(250:4000),[1 64],[0 100]);
	
	%Color vector
    col = 'brgcmykb';

    %Plot states
	figure(2);
	for i = 1:5,
        subplot(5,1,i);
        tmp = sprintf('%s-',col(i));
        plot(time_vec(250:2000),gamma(i,:),tmp);
        hold on;
        title(sprintf('State %d',i));
	end;
	
	%Plot activation in state
	figure(3);
    top = max(max(Q(:,2:65)));
	for i = 1:5,
        subplot(5,1,i);
        tmp = sprintf('%s-',col(i));
        bar([1:64],Q(i,2:65),tmp);
        axis([1 64 0 0.03]);
        title(sprintf('Fire Prob. State %d',i));
	end;

return;