function rss = met_nw_rss(p,t_data,Na_data)

global Np_init
global Na_init

Ca = p


tstart = 0;
tend = 10;

v_init(1) = Np_init;
v_init(2) = Na_init;

[t,v] = ode45(@met_network2,[tstart t_data tend],v_init,[],Ca);

Na_pred = v(:,2);

rss = 0;

for i = 2:length(Na_pred)-1
    
    rss = rss + (Na_pred(i) - Na_data(i-1)).^2;
end

rss