What happen with my code????, please help me fix this!!

When I test f=1.991, i have result Fkq=1.9908

But when I use 'for' i have result Fkq=1.9523 (f=1.925) and Fkq=1.9524 (f=1.9530) then it stop increase f.

Please someone help me resolve this problem, why f do not increase 1.954,1.955..... This is my code:

w= [40000 120000 200000 220000 100000]; gama= [0.0995 0.0995 0.0995 0.1972 0.7854];
n=length(w); c=500;
b=40; phi=0.314; SFm=0; SFt=0;
%f=1.991;
for f=1:0.001:2
  for i=1:n
      Fm(i)=w(i)*tan(gama(i));
      T(i)=(c*b*sec(gama(i))+w(i)*sec(gama(i))*tan(phi))/(f+tan(gama(i))*tan(phi));
      Ft(i)=(((c*b*sec(gama(i))+(w(i)*sec(gama(i))-T(i)*tan(gama(i)))*tan(phi))*sec(gama(i))));
  end
  for i=1:n
      SFm = SFm + Fm(i);
      SFt = SFt + Ft(i);
  end
  Fkq = SFt/SFm;  
  if abs(f-Fkq) <= 0.001
      disp(Fkq)     
  end
end

[Stephen23]

Floating-point numbers have been explained many times on this forum:

https://www.mathworks.com/matlabcentral/answers/57444-faq-why-is-0-3-0-2-0-1-not-equal-to-zero

https://www.mathworks.com/matlabcentral/answers/69-why-does-1-2-3-1-3-not-equal-zero#answer_271

https://www.mathworks.com/matlabcentral/answers/316889-why-two-equal-numbers-are-not-equal

https://www.mathworks.com/matlabcentral/answers/41536-why-is-2-24-10-22-4-not-equal-to-0

https://www.mathworks.com/matlabcentral/answers/20830-1-0010-does-not-equal-1-0010

etc, etc, etc

Answer 1

[James Tursa]

See if this link helps you understand what is going on with your floating point calculations:

http://matlab.wikia.com/wiki/FAQ#Why_is_0.3_-_0.2_-_0.1_.28or_similar.29_not_equal_to_zero.3F

[nam tran]

Thank you. But i so bad english and i just begin use matlab. I can't understand that post in your link. Can you help me directly on my code.... :(. Thank you again.