So I am using ODE45 to solve a system of equations, and then I plot the solution. The problem I am having is that for each y value, the corresponding x value is a decimal in most cases. Therefore the datetime function is not working, and I cannot seem to edit the xtick values to account for it.
z=0;
Z=0;
i = [ 1 2 3 4 5];
it =[0 28 56 91 119];
j = [3 4 5 6 7];
jt = [0 28 56 91 119];
kt = [8 9 8 7 6];
k = [0 28 56 91 119];
TSPAN = [0 300];
ICC = 500;
for P_1 = [1]
for P_2 = [2]
for P_3 = [3]
[t,y] = ode45(@(t,y) ODE(t,y,P_1,P_2,P_3,it,i,jt,j,kt,k), TSPAN, ICC);
z=z+1;
Z=Z+1;
OutT{z,:} = [t];
OutY{Z,:} = [y];
end
end
end
plot(t,y);
with the ode file
function dydt = ODE(t,y,P_1,P_2,P_3,it,i,jt,j,k,kt)
i = interp1(it, i, t, 'linear', 'extrap');
j = interp1(jt, j, t, 'linear', 'extrap');
k = interp1(kt, k, t, 'linear', 'extrap');
dydt = P_1*i-P_2*j*y-P_3*k*y;
end
The dates that I want to plot correspond to [0 28 56 91 119], and they are;
07-Jul-2015
04-Aug-2015
01-Sep-2015
06-Oct-2015
03-Nov-2015
Any ideas how I can get this to work?
