disp(' ');
disp(' RD_N_point.m  ver 1.0  March 12, 2013 ');
disp(' by Tom Irvine ');
disp(' ');
%
disp(' Enter type ');
disp('  1=calculate ratio for N');
disp('  2=calculate N for ratio ');
ic=input(' ');
%
disp(' ');
%
if(ic==1)
    N=input('Enter N ');
    a=(6.05-log10(N))/6.4;
    rd_Z_ratio=10^a;
%
    out1=sprintf(' rd_Z_ratio=%8.4g ',rd_Z_ratio);
%
else
    rd_Z_ratio=input('Enter rd_Z_ratio ');
    a=-(6.4*log10(rd_Z_ratio)-6.05);
    N=10^a;
%
    out1=sprintf(' N=%8.4g ',N);
%
end
%
disp(' ');
disp(out1);