Home > SPlaC v1_0 > Mie > DIP > DipIncEefn1fn0.m

DipIncEefn1fn0

PURPOSE ^

Calculates e_{n,1}, f_{n,1}, and f_{n,0} for n=1:nNmax for the incident dipole field

SYNOPSIS ^

function stDipEefn1fn0=DipIncEefn1fn0(nNmax,xp)

DESCRIPTION ^

 Calculates e_{n,1}, f_{n,1}, and f_{n,0} for n=1:nNmax for the incident dipole field
 using Eqs. H.88 and H.90

 Parameters:
 - nNmax:  integer scalar
 - xp:     column vector [L x 1]
           wavelength-dependent xp=kM*Rp (Rp is distance from origin)

 Returns: structure with 3 fields
          each a matrix [L x nNmax]
 - stDipEefn1fn0.en1: e_{n,1}
 - stDipEefn1fn0.fn1: f_{n,1}
 - stDipEefn1fn0.fn0: f_{n,0}

 This file is part of the SPlaC v1.0 package (copyright 2008)
 Check the README file for further information

EXAMPLE OF OUTPUT ^

Example figure output

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 
0002 stZnAll=GenZnAll(nNmax,xp,'j');
0003 % Z0 is psi(xp)/xp, Z2 is (x psi)'(xp)/xp
0004 % Z1 is psi(xp)/xp^2
0005 n=1:nNmax;
0006 
0007 % Parallel dipole
0008 coeff1=3/8*(2*n+1); % [1 x nNmax]
0009 nCoeffmat=repmat(sqrt(coeff1),length(xp),1); % [L x nNmax]
0010 
0011 stDipEefn1fn0.en1 = i * nCoeffmat .* stZnAll.Z0;
0012 stDipEefn1fn0.fn1 = - nCoeffmat .* stZnAll.Z2;
0013 
0014 % Perpendicular dipole
0015 coeff2=4*coeff1.*n.*(n+1); % [1 x nNmax]
0016 nCoeffmat=repmat(sqrt(coeff2),length(xp),1); % [L x nNmax]
0017 
0018 stDipEefn1fn0.fn0 = nCoeffmat .* stZnAll.Z1;

This web page is part of the SPlaC package © 2008. Contact: Eric Le Ru
Generated on Wed 03-Dec-2008 11:10:14 by m2html © 2003 (adapted)