function f = theoint(n, region, itmax)

% Theodorsen's method for interior regions
global ratio

% Initial guess for boundary correspondence
th = 2*pi*[0:n-1]/ n;
phi = th;
phil = phi;
disp('Iteration no.  Error between successive iterates');
z=exp(i*th);
fex=2*ratio*z./(1+ratio-(1-ratio)*z.^2);
f = bdrytheo(region,phi);
% Start of Theodorsen loop
for it = 1 : itmax
    c = log(abs(f));
  % Conjugation 
    c = conjug(c);
  % Compute new boundary corespondence and calcuate errors
    phi = real(c) + th;
    error=max(abs(phi-phil));
    errex=max(abs(f-fex));
    phil=phi;
    fprintf('%6.0f%30.14e%30.14e\n',it,error,errex);
    f = bdrytheo(region, phi);
end
% End of Theodorsen loop