function [option] = make_plot(nvars,species,iterations,ninst,AllInst,trials)
color = ['k' 'b' 'r' 'g' 'm' 'c'];
line = {'-' '-.' '--' ':'};
nS = 1;
inp = 1;
while inp ~= nvars+2 && inp ~= nvars+3
    disp('Choose a species for the time plot');
    for ind = 1:nvars
        disp([int2str(ind) '. ' species{ind} ]);
    end
    disp([int2str(nvars+1) '. Finished Entering']);
    disp([int2str(nvars+2) '. Exit time plot']);
    disp([int2str(nvars+3) '. Exit program']);
    disp('select the corresponding numbers to plot');
    inp = 1;
    in = 1;
    clear nSpecStr
    clear y
    names = {};
    while inp ~= nvars + 1 && inp ~= nvars + 2 && inp ~= nvars + 3
        nSpecStr{in} = input('Enter number: ');
        inp = nSpecStr{in};
        in = in + 1;
    end
    if inp ~= nvars + 2 && inp ~= nvars + 3

        wr = input('Pick the row that the well is on: ');
        wc = input('Pick the column the column that the well is on: ');
        figure
        paint = 1;
        type = 1;
        for num = 1:length(nSpecStr)-1
            nS = nSpecStr{num};
            names{num} = species{nS};
            
            for trial = 1:trials
                Ttime = 1;
                for SDiff = 1:iterations
                    for time = 1:ninst
                        y(Ttime,trial) = AllInst{trial}{nS,time,SDiff}(wr,wc);
                        Ttime = Ttime + 1;
                    end
                end
            end

            %"if" statements are used to vary the colors and
            %line types on the plots
            if paint > 6
                paint = 1;
                type = type + 1;
            end
            if type > 4
                type = 1;
            end
            for trial = 1:trials
                plot(y(:,trial),[line{type} color(paint)])
                hold on
            end
            hold off
                paint = paint + 1;
            hold on
            

        end
        
        title(['The Effect on the Species in Well(' int2str(wr) ',' int2str(wc) ') for all Diffusions' ])
        xlabel('time(min)')
        ylabel('Concentration(density)');
        legend(names)
        
  
    end
    if inp == nvars + 3
        option = 5;
    else
        option = 1;
    end
end