clear
lambda=[700 550 400];
for i=1:length(lambda)
    disp(i)
   [X,Y,I(:,:,i),x0,y0]=singlescattering(lambda(i));
end    
%MAX=max(max(max(I)));
%I=I/MAX;
V=rgb2gray(I);
figure
caxis([0 2.5])
h=plot(x0,y0,'oy')
set(h,'MarkerFaceColor','y') 

[n,m]=size(X);
II=I;
for i=1:n
    for j=1:m
        MAX=max(max(I(i,j,:)));
        if MAX>0
         II(i,j,:)=I(i,j,:)/MAX;
        end 
    end
end
V=rgb2gray(II);
[V,map]=rgb2ind(II,65536);
V=double(V);
h=polar([0 2*pi], [0 90],'b');
delete(h);
hold on;
H=findall(gca);
for i=1:length(H);
    try
        set(H(i),'FaceColor',[0 1 0])
    end
end  
h=plot(x0,y0,'oy')
set(h,'MarkerFaceColor','y') 



for i=1:n
    for j=1:m
     dist=sqrt((X(i,j)-x0)^2        
    end
end




pcolor(X,Y,V); shading interp; 
colormap(map);
%h=plot(x0,y0,'oy')
%set(h,'MarkerFaceColor','y') 
%caxis([0 3.5]);