How can I add the solutions of xd as a matrix , I tried with xd1 , but it still give an error. and how can I make the new value as a initial value for the next st. x0 known and kmax =10. While the solution should be a matrix of 10 rows .. Thank you !

for k=1:kmax t1=k*h; xd=modeu(Fdfpek,t1,x0,h); xd1=[xd'; x0]; x0=xd; end

Answer 1

[Mischa Kim]

Ammar, use

 for k = 1:kmax 
    t1    = k*h; 
    xd(k) = modeu(Fdfpek,t1,x0,h);
    x0    = xd(k); 
 end