MATLAB: Integrate for a specific period of time

euler

Please help me. I am trying to use Euler integration to integrate for 10 seconds with a step size of .01 seconds. Plot x versus time.
x(0) = 1
t=0:.01:10;
x0=1;
xdot=-2*(x^3)+sin(0.5*t)*x;
for t=0:0.01:10
    x=integrate(xdot,t,x0);
end
plot(t,x)

Best Answer

t=0:.01:10;
x = zeros(numel(t));
x(1) = 1;
fun_xdot = @(t,x) -2*(x^3) + sin(0.5*t)*x;
for i = 1:numel(t)-1
    x(i+1) = x(i) + (t(i+1)-t(i))*fun_xdot(t(i),x(i));  
end
plot(t,x)
Related Question