Documentation of findzeros_of_dot

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

tau = 0.5; nt = length(tau);

Cross-Reference Information

Listing of script findzeros_of_dot

clean
cd /home/disk/tao/dvimont/matlab/BAT
x = 2*pi*[0.005:.01:.9];
y = sigifun(x, 3.9, 2.2, 0.5);
figure(1); fo(1); clf;
subplot(2,1,1);
plot(x, y);
b = .1:.1:6; nb = length(b);
c = -4:.1:6; nc = length(c);
tau = 0.02:0.02:1.26; nt = length(tau);

x0 = NaN*ones([nb nc 1]);

opts = struct('TolX', 1e-4);

for i = 1:nb;
  x0
  for j = 1:nc;
    for k = 1:nt;
      x0(i,j,k) = fsolve('sigifun', 1e-4, opts, ...
			 b(i), c(j), tau(k));
    end
  end
end
xi = x0;

x0 = xi;
%x0(x0 < 0.0013) = 0;

clean
cd /home/disk/tao/dvimont/matlab/BAT/data
load dot_sigimag.mat

x0 = xi;

nb = length(b);
nc = length(c);
nt = length(tau);

xr1 = NaN*ones(nb, nc);
xr2 = NaN*ones(nb, nc);
lind = find(tau == 0.5);
for i = 1:nb;
  for j = 1:nc;
    for k = lind:lind;
      xr1(i,j) = c(j)-x0(i,j,k)./tan(x0(i,j,k)*tau(k));
      xr2(i,j) = (1./tau(k))*log(b(i)*sin(x0(i,j,k)*tau(k))/...
				       x0(i,j,k));
      if x0(i,j,k)<0; xr1(i,j)=0; xr2(i,j)=0; end;
    end
  end
end


cd /home/disk/tao/dvimont/matlab/BAT/data
save dot_sigimagmat  xi b c tau
save dot_params.mat  xi xr1 xr2 b c tau

figure(1); fo(1); clf;
bind = find(b == 2.25);
subplot(4,2,1);
  pncont(b, c, x0'/(2*pi), [-1:0.1:5], 0, 'k');
subplot(4,2,3);
  pncont(b, c, xr2', [-5:.25:5], 0, 'k');
subplot(4,2,5);
  pncont(b, c, real(xr1)', [-5:.25:5], 0, 'k');
subplot(4,2,7);
  pncont(b, c, 0.5*(xr1'+xr2'), [-5:.25:5], 0, 'k');