Help with a loop that computes at each step

hi im doing an assignmenr where i have 5 different stepsizes that compose a row vector hs. i have 100/hs steps for which id like to compute the value th2 at each stepsize.

th1=1;
th2=1;
w1=0;
w2=0;
hs(1)=[0.05];
hs(2)=[0.05/2];
hs(3)=[0.05/4];
hs(4)=[0.05/8];
hs(5)=[1/1000];
th2s=[];

????????????
    for i=1:N
        k1= h*fpend(y);
        k2=h*fpend(y + k1/2);
        k3=h*fpend(y + k2/2);
        k4= h*fpend(y + k3);
        y = y + (k1 + 2*k2 + 2*k3 + k4) / 6;
    end
    th2s=[th2s,y(2)]

[Awais Saeed]

Your question is not clear. Be specific. Why do you have a 2nd end at the end?

[Lavorizia Vaughn]

i have removed the 2nd end. i still need something relating h to hs

[Awais Saeed]

You still did not made the question clear. what is h and y and ?????????????

[Lavorizia Vaughn]

y are approximations using runge kutta. h is the step size

[Lavorizia Vaughn]

the assignment says,"Run Case 1 in problem 2 with the five stepsizes h = 0.05/2(k−1), k = 1, 2, 3, 4, and h = 0.001. Compute the value of θ2(t = 100) for each stepsize. '

[Jan]

The instruction is nonsense: For k=1, the value of h = 0.05/2(k-1) is Inf.

Answer 1

[Jan]

hs = [0.05, 0.05/2, 0.05/4, 0.05/8, 1/1000];
for k = 1:numel(hs)
   h = hs(k);
   
   N = totalDuration / h; % ???
   y = y0;
   for i = 1:N
        k1 = h * fpend(y);
        k2 = h * fpend(y + k1/2);
        k3 = h * fpend(y + k2/2);
        k4 = h * fpend(y + k3);
        y  = y + (k1 + 2*k2 + 2*k3 + k4) / 6;
   end
   
   th2s(k) = y;  % Why y(2)?
end

[Lavorizia Vaughn]

You helped me with the other but very similar thread. thanks a bunch. your help is much appreciated.