So I have a system of ODEs that I am solving through NDSolve, and I can plot them in mathematica, but I am curious how to get the x and y values used in the plotsfor manipulation in programs like excel. Here is the code I am currently working with. Clear["Global'*"];
ka = .004; kd = 1.2; km = 1500; ns = 25; sens = 1; ro2 = 11;
odecella = ka nut[t] cella[t] - (kd cella[t])/(1 + (nut[t]/ns)^sens);
odecelld = (kd cella[t])/(1 + (nut[t]/ns)^sens);
odeo2 = ro2 (1256 - o2[t]) - km cella[t];
odenut = -ka nut[t] cella[t];
vars = {cella[t], celld[t], o2[t], nut[t]};
solution =
NDSolve[{Derivative[1][cella][t] == odecella,
Derivative[1][celld][t] == odecelld,
Derivative[1][o2][t] == odeo2, Derivative[1][nut][t] == odenut,
cella[0] == .01, celld[0] == 0, o2[0] == 1110, nut[0] == 100},
vars, {t, 0, 120}];
vars;
Plot[cella[t] /. solution, {t, 0, 120}, PlotRange -> Full]
Plot[celld[t] /. solution, {t, 0, 120}, PlotRange -> Full]
Plot[o2[t] /. solution, {t, 0, 120}, PlotRange -> Full]
Plot[nut[t] /. solution, {t, 0, 120}, PlotRange -> Full]