MATLAB: “Inputs must be floats, namely single or double”

integration

system of equations is:
F = @(t, r, y) [y(1), y(2), y(3), f_integral]
And then calling the ode solver:
[t, y] = ode45(F, [0 5], [1 1 1]).

Best Answer

function main
tspan=[0 5];
y0=[1; 1; 1; 0];
[t,y]=ode45(@fun,tspan,y0)
function dy=fun(t,y)
dy=zeros(4,1);
dy(1)=y(1);
dy(2)=y(2);
dy(3)=y(3);
funint=@(r)(y(3)^2+y(4)^2)*r;
dy(4)=integral(funint,-1,1);
Best wishes
Torsten.
Related Question