for(el=1;el<=L;++el) {
  print(el);
  t=get_kex(el,1);
  plot2d(t(:,2),t(:,1));
  t=get_kex(el,2);
  plot2d(t(:,2),t(:,1));
}

;;;
