function [] = PlotdCadt(start,finish,interval,start2,finish2,interval2)
B1 = findobj(gcf,'Tag','TauCa')
TauCa = str2num(get(B1,'String'))
BPplot = zeros(floor((finish-start)/interval),floor((finish2-start2)/interval2));
BPplot = zeros(floor((finish-start)/interval),1);
BPplot(1) = 0;
h = waitbar(0,'Please wait...');
icount = 0;
jcount =0;
axisplot(1) = 0;
for i = start:interval:finish
    axisplot(length(axisplot)+1) = i;
    icount = icount +1;
    jcount = 0;
    waitbar(i/finish)
    
    for j = start2:interval2:finish2
        jcount = jcount +1;
    
    BPplot(icount,jcount) = dCaDt(i,j,TauCa)
    
    end
end
close(h)
contourf(BPplot)
xlabel('Inmda')
ylabel('Calcium concentration')
title('Rate of change of calcium concentration')
colorbar