function [lvsp] = lvdemo_ode(t,lvs,circs)
  % odefile for the multiple
  
  lvsp = zeros(size(lvs));
  

  for J = 1:length(lvs)
    sum = 0;
    for j = 1:length(lvs)
      if (j ~= J)
	sum = sum - i*circs(j)/2*pi * 1/(lvs(J)-lvs(j));
      end
      lvsp(J) = conj(sum);
    end
  end
  
  