Help with knowing what kind of vector im using

1 view (last 30 days)
could someone tell me what th2s=[th2s,y(2)] means?
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 h=hs %for all the stepsizes
y=[th1,th2,w1,w2]; we have this vetir which has the values
N=100/h;
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)]
end

Accepted Answer

Cris LaPierre
Cris LaPierre on 13 Nov 2021
It is appending y(2) to the end of the vector th2s at the end of each h loop.
th2s=[];
for a=1:3
a
th2s = [th2s,a]
end
a = 1
th2s = 1
a = 2
th2s = 1×2
1 2
a = 3
th2s = 1×3
1 2 3
  3 Comments
Lavorizia Vaughn
Lavorizia Vaughn on 13 Nov 2021
i am doing a rk4. and i wish to find a better way. this way seems a bit hard.

Sign in to comment.

More Answers (0)

Products


Release

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by