% plot likelihood

n = length(rho);

lik_1 = zeros(n,1);

i = 1;

while i <= n
   
   r = rho(i);
   
   lik_1(i) = sar_lik(r,low_1,high_1,val_1,W1,y,X);
   
   i = i + 1;
   
end