Dear All,
Can someone enlighten me by explaining why I'm getting following error message in code below?
Index exceeds matrix dimensions. Error in verysimple3 (line 14) [Tp]=ode45(@(l,Tp) myODE(l,Tp,Aap,Acp),Pool(k:k+1),y0);
If I add the parameters Aap,Acp as input parameters of my anonymous function then I get: Not enough input arguments. Error in verysimple3>@(l,Tp,Aap,Acp)myODE(l,Tp,Aap,Acp)
Sorry for this question with a probably very obvious answer! Originally I didn't want to include Aap and Acp as input parameters for myODE function as these are fixed for each interval (see loop) but I'm not sure to have another option in this case…
a =2;b =1;Pool = [0,a,b];y0=30;for k=1:numel(Pool) switch k case rem(k,2)==1 Aap=a; Acp=Aap; otherwise Aap=2*b; Acp=0; end [l,Tp]=ode45(@(l,Tp) myODE(l,Tp,Aap,Acp),Pool(k:k+1),y0);y0=y(end,1);endplot(l,Tp(:,:));function dTpdl=myODE(l,Tp,Aap,Acp)dTpdl=(140-Tp)*Acp+(100-Tp)*Aap-Aap/(Tp+273); end
Best Answer