18 ## @deftypefn  {Function File} {@var{JM} = } jones_waveplate ()
19 ## @deftypefnx {Function File} {@var{JM} = } jones_waveplate (@var{p})
29 ## cell array. In the two latter cases, a cell array @var{JM} of
46 function JM = jones_waveplate(varargin)
60     JM = cell(size(phase_in_lambda_units));  variable
61     JM_subs = cell(1,ndims(JM));
62     for jmi=1:numel(JM)
63       [JM_subs{:}] = ind2sub(size(JM),jmi);
64       JM{JM_subs{:}} = s_waveplate(phase_in_lambda_units(JM_subs{:}));
69     JM = s_waveplate(phase_in_lambda_units);  variable
76 function JM = s_waveplate(phase_in_lambda_units)
78   JM = zeros(2,2);  variable
79   JM(1,1) = 1;
80   JM(2,2) = exp(-1i*phase_in_lambda_units*2*pi());