Y_bar = (OMG.^2)/(1-OMG.^2)
That overwrites all of Y_bar each time. You need
Y_bar(1,i) = (OMG.^2)/(1-OMG.^2)
But your code could be simplified:
Y_bar = omega_bar.^2 ./ (1-omega_bar.^2);
No loop is needed.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!