Программа вычисление золотого сечения, Циклы Опции

[Casio]

Здраствуйте. Помогите решить след. задачу:
Составить программу вычисления золотого сечения q с точностью Е=0.00001, используя представление q в виде:
Программа строится на вычислении предыдущего z1 и нового значения z2. Итерационный

Код

q=1 + ______1_________
      2 + ____1_______
           2 + ___1___
               2 + ...

цикл заканчивается, если |z2-z1|<=E
Z1 начальное =2+2/5
Z2=2+1/z1

[volvo]

Так пойдет?

program sistema;
uses crt;
const
  da=0.3;
  dx=2;

function max(a, b: real): real;
begin
  max := a;
  if b > a then max := b;
end;

var
  a, max_a,
  p, w1, w2, w3, maximum, next: real;
  i, x, max_x: word;
begin
  clrscr;
  a:=1.5;

  maximum := -maxint;
  repeat

    x:=8;
    repeat

      w1:=x*sqrt(x-a);
      w2:=x*sin(a*x);
      w3:= Exp(x*Ln(a))/cos(sqr(a*x-1))*cos(sqr(a*x-1))*cos(sqr(a*x-1));

      next := max(max(w1, w2), w3);

      if maximum < next then begin
        maximum := next;
        max_a := a; max_x := x
      end;

      x:=x+dx;

    until(x>10);

    a:=a+da;
  until(a>3);

  writeln('W max = ',maximum:2:1);
  writeln(' a = ', max_a:2:0, '; x = ',max_x);
  readln;
end.