  function conductor
% function conductor
%
% Plots the current distribution in a circular wire with a skin depth of
% one tenth of the radius, both the current in phase with the total current
% and the current in quadrature with the total current.
%
r   =  (0:0.01:1); % 101 sampling points radially
tau =  (1 - i).*10.;
J_z =  tau.*BesselJ(0,tau.*r)./(2.*pi.*BesselJ(1,tau));
plot (r, real(J_z)) % current in phase with total current
hold on
plot (r, imag(J_z), 'k:') % quarter period later
title  ('Current Distribution in Circular Wire')
xlabel ('Radial Position')
ylabel ('Current Density')
hold off
return
