1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
clc
clear
x=[0 0.1 0.2 0.3 0.4 0.5]
y(1,1)=-1
for i=1:5
B(1,i)=(-2*x(1,i)-y(1,i))
C(1,i)=0.1*B(1,i)
y(1,i+1)=y(1,i)+C(1,i)
D(1,i)=y(1,i+1)
E(1,i)=0.1*B(1,i)
F(1,i)=0.1*(-2*x(1,i+1)-y(1,i+1))
G(1,i)=(E(1,i)+F(1,i))/2
y(1,i+1)=y(1,i)+G(1,i)
end
for i=1:5
T=[x(1,i), y(1,i), C(1,i),D(1,i),F(1,i),G(1,i),y(1,i+1)]
disp(T)
end
|
