Vo = [30 15 13 10];
MaxIter=1e10;
Cntr=0;
while true
Cntr=Cntr+1;
Re1 = (D(1).*Vo(1)) / nu;
Re2 = (D(2).*Vo(2)) / nu;
Re3 = (D(3).*Vo(3)) / nu;
Re4 = (D(4).*Vo(4)) / nu;
A = [Matrix];
b = [Matrix];
Vn = A\b;
for i = 1:4
Va(i) = (alpha * Vn(i)) + ((1-alpha) * Vo(i));
end
if all(abs(Va-Vn(1:4))<1e-10) && all(abs(Va-Vo)<1e-10)
break;
end
if Cntr==MaxIter
break;
end
Vo=Va;
end
Best Answer