function [EPSPnmdaT] = EPSPnmda(T,PreTrain,deltaT,V,TauF,TauS,Vr,I_f,P0,Gnmda,BP,EPSP,lastone)


T = T/10;
deltaT = deltaT/10;
%TauF --> 50 and TauS --> 200
%Vr --> 0 mV
%I_f --> 0.5
%I_f = 0.5;
I_s = 1 - I_f;
Theta = 0;
EPSPnmdaT = 0;
%index = find((preTrain/1000) <= T,1,'last');
%indexs = find(((PostTrain/10) <= T+deltaT)&((PostTrain/10) > T));
%indexs = find(((PreTrain/10) <= T+deltaT));

%if((T+deltaT) - PreTrain(find((PreTrain/10) < T+deltaT, 1, 'last'))) > 10;
%Thetan = 0;    
%else
%Thetan = ((I_f*exp((-deltaT)/TauF))+ (I_s*exp(-(deltaT)/TauS)));
%end
%Vn = BP(1) + EPSP(1);
%InmdaT = (P0 * Gnmda*Thetan*BV(Vn)*(Vn-Vr));


    indexs = find((PreTrain/10) <= (T+deltaT),1,'last');

%for i = 1:length(indexs)
if isempty(indexs) == false
Theta = Theta + ((I_f*exp((-((T+deltaT)-(PreTrain(indexs)/10)))/TauF))+ (I_s*exp(-((T+deltaT)-(PreTrain(indexs)/10))/TauS)));
end
    %Theta = ((I_f*exp((-deltaT)/TauF))+ (I_s*exp(-(deltaT)/TauS)));

    



        %Vt = BP(i+1)+EPSP(i+1);
        
    %end

EPSPnmdaT = EPSPnmdaT + (P0 * Gnmda*Theta*BV(V)*(V-Vr));
if (PreTrain/10) == (T+deltaT)
    EPSPnmda = EPSPnmda + lastone;
end