theta_k = pi/2/G + (1:G)*pi/G [1 60] h =0.2 discarding the less reliable data points index = find(radius>=h) [1 50] radius = radius(index) [1 50] theta = theta(index) [1 50]-theta_k = pi/2/G+ (1:G)*pi/G [1 60] h =0.2 discarding the l