clear
clc
 
time=[0 10];
y=[0 1]';
[t1 y1]=ode45('finode3',time,y);


sum=0;
b=0;


for i=1:1:5

    for k=0.0000000000001:1:i
    a=1;
    b=0;
    b =b+((y1(max(find(t1<k)))+y1(max(find(t1<k+a))))*a)/2;
    
   
 
    end
    
     sum=sum+b

     plot(t1,sum);
end
    
---------------------------------------------------------------------
副程式


function dy =finode3 (t,y)
dy(1)=y(2);
dy(2)=-y(1);
dy=dy';
end

-----------------------------------------------------------------------

我想繪製出時間和迴圈跑出SUM的那些值  0.3992,1.2677,1.8367  等等.......

但是我在workspace裡看到的sum只有一個值,而不是一個向量,不知道該怎麼會出圖出來

懇請各位大大可以教小弟一下,謝謝