I believe you need to read about for loop and subscripting. Try this:
As1=510; fy=420; h = 650; d = 590; Es = 200000; fc = 21; As2 = 520; b = 350; beta1 = .85;ecu=1;%couldn't find ecu, randomly given c = [0.0001:.1:1]*h; for i=1:numel(c) F = min(fy,max( -fy,((d-c(i))./c(i)*ecu).*Es)).*As1/1000; Cc = beta1*fc*b*.85*c(i)/1000; S(:,i) = sum(F)-Cc; end elementS = find(S>0,1) neededc = c(1,elementS)