clear all; close all; clf;
%Fmax=maximum flux for calcium pump 
%H=Hill coefficient 
%Km=Km is the [Ca]i at which Jpump is at half of Fmax
%Ccai=[ca]i---> Intracellular calcium concentration
%Jpump=flux for calcium pump in micro.moles/sec/cm2
%Jpump=Fmax/(1+(Km/[ca]i)^H)

%-----------------calcium pump Without loop-----------------%
Fmax=5e-6; % micromoles/sec/cm2
H=0.75;
Km=260;
Ccai=0.5:-0.001:0.001;
alpha=Km./Ccai;
Jpump=Fmax./(1+(alpha.^H));
plot(Ccai,Jpump);
set(gca,'XDir','rev')


%-----------------calcium pump With loop (For Check only- Do not use this in your code)-----------------%
% % dt=0.0002; t_final=0.3;
% % Fmax=5e-6; % micromoles/sec/cm2
% % H=0.75;
% % Km=260;
% % Ccai=0.5;
% % tt=0
% % for t=0:dt:t_final
% % tt=tt+1;    
% % alpha=Km/Ccai;
% % Jpump(tt)=Fmax/(1+(alpha^H));
% % Ccai=Ccai+(-1)*(1e-3*Jpump(tt));
% % end
% % t=0:dt:t_final;
% % plot(t,Jpump);xlabel('time(S)'); ylabel('Jpump(uM/S)');