complexzeros: v=fresnels(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=fresnels(z,N)

function v=fresnels(z,N)
% v=fresnels(z,N) sums the Bessel series expansion of 
% the Fresnel sine integral. N+1 terms are used in
% the series with a default value of N=100. 
if nargin<2, N=100; end
s=size(z); z=pi/2*z(:).^2; n=1.5+2*(0:N);
v=sum(besselj(n,z),2);
v=reshape(v,s);

Highlights from complexzeros

Highlights from
complexzeros