function [ax, co] = rock1(T)
% ROCK1 returns complex axial and coplanar field as a function of T
% for uniform infinite conducting medium - i.e. no boundary problem

% T is the ratio of distance to transition distance, i.e. T = kr

ax = 2*(1 + (1+i)*T);
co = 1 + (1+i)*T +2*i*T.^2;

scale = exp(1)/(2 * sqrt(5));

ax = scale * ax .* exp(-i*T) .* exp(-T) ./ T.^3;
co = scale * co .* exp(-i*T) .* exp(-T) ./ T.^3;

% the values are scaled so that the axial field has unity magnitude
% at T=1
