complexzeros: v=erfs(z,N) - File Exchange

Code covered by the BSD License

Highlights from
complexzeros

[Zeros,Fzeros]=muller(fn,npre... [Zeros,Fzeros]=muller(fn,nprev,maxit,ep1,ep2,fnreal,Zeros)
[rts,frts,rtsm,frtsm,func]=in... [rts,frts,rtsm,frtsm,func]=interact(func,xlim,ylim)
f=zetac(x) f=zetac(x) computes the Riemann zeta function
v=erfs(z,N) v=erfs(z,N) evaluates the series expansion for the
v=fresnelc(z,N) v=fresnelc(z,N) sums the Bessel series expansion of
v=fresnels(z,N) v=fresnels(z,N) sums the Bessel series expansion of
v=gammac(z) v=gammac(z) gives the gamma function for complex argument.
readme.m
runexamples.m
View all files

from complexzeros by Howard Wilson
Complex zeros of a function of z=x+i*y are computed using graphics and Muller's method.

v=erfs(z,N)

function v=erfs(z,N)
% v=erfs(z,N) evaluates the series expansion for the
% error function of  complex argument scaled by 
% exp(-abs(z)). The default value of N is 50.
if nargin<2, N=50; end
s=size(z); z2=z(:).^2; k=0:N;
v=(besseli(2*k+.5,z2,1)-besseli(2*k+1.5,z2,1))*(-1).^k';
v=reshape(v,s);

Highlights from complexzeros

Highlights from
complexzeros