function colicin_char(GFP_in, OD_in, time, anz, ratio, graph)

for i=1:1:anz
    GFP(1:8, 1:12, i) = GFP_in((i-1)*8+1:i*8, 1:12);
    OD(1:8, 1:12, i) = OD_in((i-1)*8+1:i*8, 1:12);
end

if (graph == 1)
for l=1:1:6
 figure

    hold all
    for k=1:1:7
        for j=1:1:anz
            GFP_plot(k,j)=GFP(k,l,j);
        end
    plot(time(:,1),GFP_plot(k,:), 'LineWidth', 2);
    end
    legend('0 pM', '250 pM', '500 pM', '750 pM', '1 nM', '10 nM', '25 nM', 'Location', 'EastOutside');
    xlabel('time [h]');
    ylabel('GFP intensity');
    title(['Prey:Killer ' int2str(ratio(l,1)) ':1'])

 hold off
end

for t=7:1:12
 figure

    hold all
    for k=1:1:7
        for j=1:1:anz
            GFP_plot(k,j)=GFP(k,t,j);
        end
    plot(time(:,1),GFP_plot(k,:), 'LineWidth', 2);
    end
    legend('0 pM', '250 pM', '500 pM', '750 pM', '1 nM', '10 nM', '25 nM', 'Location', 'EastOutside');
    xlabel('time [h]');
    ylabel('GFP intensity');
    title(['Prey:Reference ' int2str(ratio(t-6,1)) ':1'])

 hold off
end

for f=1:1:6
 figure

    hold all
    for k=1:1:7
        for j=1:1:anz
            OD_plot(k,j)=OD(k,f,j);
        end
    plot(time(:,1),OD_plot(k,:), 'LineWidth', 2);
    end
    legend('0 pM', '250 pM', '500 pM', '750 pM', '1 nM', '10 nM', '25 nM', 'Location', 'EastOutside');
    xlabel('time [h]');
    ylabel('OD');
    title(['Prey:Killer ' int2str(ratio(f,1)) ':1'])

 hold off
end

for z=7:1:12
 figure

    hold all
    for k=1:1:7
        for j=1:1:anz
            OD_plot(k,j)=OD(k,z,j);
        end
    plot(time(:,1),OD_plot(k,:), 'LineWidth', 2);
    end
    legend('0 pM', '250 pM', '500 pM', '750 pM', '1 nM', '10 nM', '25 nM', 'Location', 'EastOutside');
    xlabel('time [h]');
    ylabel('OD');
    title(['Prey:Reference ' int2str(ratio(z-6,1)) ':1'])

 hold off
end

end

GFP_start_killer = GFP(1:8,1:12,1);
GFP_12_killer = GFP(1:8,1:12,anz);

save gfp_start_killer.txt GFP_start_killer -ASCII;
save gfp_12_killer.txt GFP_12_killer -ASCII;
        
   GFP_final(:,:)= GFP_12_killer(1:7,1:6);
   [x,y]=meshgrid([1,10,25,50,100,500],[0,250,500,750,1000,10000, 25000]);

    figure
    hold all
    surf(x,y,GFP_final(:,:));
    zlabel('GFP intensity');
    xlabel('c(AHL) [pM]');
    ylabel('Prey/Killer ratio')
    colormap(jet);
    colorbar
    shading interp        


end