call_AMPARODE_linear_CaCaMInput_scaleC1C2.m

function [t,y]= call_AMPARODE_linear_CaCaMInput_scaleC1C2(CIC,c1var, c2var, CaCaMinput)
%close all
%clear


tspan = (0:.01:500);
%initial conditions
names = ["CaMKII", .../
    "CaMKIIp", "CaNp", "CaN", "I1p", "I1", "PP1p", "PP1", .../
    "Aint", "Amem" ,"scaff", "Abound"];

CaMKII_IC = CIC;  %uM
CaMKIIp_IC = 0; %uM
CaNp_IC = 1;    %uM
CaN_IC = 0;     %uM
I1p_IC = 1.8;    %uM
I1_IC = 0;     %uM
PP1p_IC = 0.25;%0.27;   %uM
PP1_IC = 0;    %uM

Aint_IC = 0.1597;%0.2628;     %uM
%%%%=
Amem_IC = (7.288/(6.022e23))*1e12;      %moles/m2 12
%%%%
scaff_IC = (4000/(6.022e23))*1e12;     %moles/m2
Abound_IC = (1000/(6.022e23))*1e12;    %moles/m2

original_cond = [CaMKII_IC CaMKIIp_IC CaNp_IC CaN_IC I1p_IC I1_IC PP1p_IC PP1_IC.../
    Aint_IC Amem_IC scaff_IC Abound_IC];
[t,y] = ode23s(@(t,y) AMPARODE_linear_CaCaMInput_scaleC1C2(t,y,c1var, c2var, CaCaMinput,tspan), tspan, original_cond);


end