% kaustika, optiskt fenomen i en ring
  clear, clf
  dfi=2*pi/100;  fi=(0:dfi:2*pi)';  plot(cos(fi),sin(fi)), hold on
  v=(0:dfi:pi/2)';                             % vinkeln fram till pi/2
  xk=cos(v).^3; yk=sin(v).*(3/2-sin(v).^2);    % x och y-värden för kaustikan
  plot(xk,yk, -xk,yk), axis equal
% Fyll hela kaustikan med gult, lagra koordinater runt om!
  u=(0:dfi:pi)';               
  X=[cos(u); -xk; flipud(xk)];   
  Y=[sin(u);  yk; flipud(yk)];  fill(X,Y,'y')  % fyll kaustikan med gult


  

  
         

  
  
