%PLOT_F yields a plot of functional values
%
%INPUTS; must specify domain [a,b] and number of points, k .  

a = -.99 ;

b = .99 ;

k = 50 ;

r = (b - a)/k ;

Xval= [a:r:b]' ;

Yval = zeros(k+1,1) ;

for i = 1:k+1
   
   Yval(i) = car_lik(Xval(i),WW,y,X);  
  
end

plot(Xval,Yval)