%% Parameter Scanning
% A brief section of code that will call our simulation function and scan
% over a range for a given parameter. 
% See <http://2015.igem.org/Team:Exeter/Modelling_Parameter> for more information.

function parameterscan()
    GFPovertemps=[];
    T=273:5:343; %0 degrees to 69
    r=[20 100 200 400];
    t=[20];
    Doubleparameterscan=[];
    
    for j=1:size(r,2);
        for i=1:size(T,2);
            Doubleparameterscan(i,j)=BM_07_09_insilicobinding(t(1),r(j),1000,T(i)); 
    %         xtemps[i]=char(T(i)-273);
        end
    end

    
    savefile= 'Doubleparameterscan';
    save(savefile, 'Doubleparameterscan')
    Doubleparameterscan; %Use this saved array to plot a Surface
end
