syms x y z; for i = 2:10 sol = solve(x/a(i)==y/b(i),y/b(i)==z/c(i),x*x+y*y==559*(z+300),x,y,z); aCopy(i) = sol.x(1); bCopy(i) = sol.y(1); cCopy(i) = sol.z(1); end