Hi, I'm trying to write an iteration program that runs a while loop until the variable 'theta' is almost the same as the output x (within a certain accuracy). I keep getting this error, but when I try to change the code to run, it only runs for the first loop. Can anyone help? My code is below (i've been using 0.01 as thetaInitial because i know this is pretty close to the correct answer)
function x=iteration(thetaInitial)
theta(1)=thetaInitial;
x(1)=0.0248;
error=1e-8;
Re=10000;
i=1;
while abs(x(i)-theta(i))>error
x(i)=(((6.25*log(2))+1.5)*log(Re*(theta(i)^0.5))+0.09)^-1; i=i+1;
end
x=x(1:i+1);
plot(x)
[x(end) i]
Best Answer