ThetaPlot = [40,50];   % Put incident angle as input

Count1 = 0;
Count2 = 0;
NewColor = 0;

for m = 1:length(ThetaPlot)
    NewColor = 0.1;
    Title = sprintf('Angle(%d)',ThetaPlot(m));
    figure('name',Title)
    Count1 = 0;
    for i = 1:length(xrange)-1             % X Coordinate of Start Node
        for j = 1:length(yrange)-1         % Y Coordinate of Start Node
            Count1 = Count1+1;
            for k = 2:length(xrange)       % X Coordinate of End Node
                for l = 2:length(yrange)   % Y Coordinate of End Node
                    Count2 = Count2+1;
                    if ReflRayData(Count2,ThetaPlot(m) - ThetaIncident(1) + 1,Count1) == 1
                        x = [xrange(i),xrange(k)];
                        y = [yrange(j),yrange(l)];
                        plot(x,y,['-',[0,NewColor,0],'o']);
                        hold on;
                    end
                end
                Count2 = 0;
            end
        end
    end
    NewColor = NewColor + 0.02;
end