Delphi, ошибка деления на ноль, хотя знаменатель не ноль

[Эренка]

Здравствуйте, так вышло, что по учебе пишу программу на Delphi, но его не понимаю,
процедура:
s:=0.5; Ta:=1000; n2=100;

procedure calc_v0;
  var
    a,b,c,lam,i,a1,a2,a3:variant;
    j:integer;
  begin
  a:=VarComplexCreate;
  a1:=VarComplexCreate;
  a2:=VarComplexCreate;
  a3:=VarComplexCreate;
  b:=VarComplexCreate;
  c:=VarComplexCreate;
  i:=VarComplexCreate;
  i.Imaginary:=1;
  i.real:=0;
  lam:=VarComplexCreate;
  lam:=sqrt(sqrt(0.25*Ta))*(1+i);
  for j:=0 to n2 do
    begin
      a1:={Sech(2*lam)*(}lam*lam+2*exp(s*(1+j*dz))*s*lam*lam+exp(2*s)*(s*s-lam*lam){)*cosh(2*lam)};
      a2:=2*((1+2*s)*lam*lam+exp(2*s)*(s*s-lam*lam))*Cosh(lam-j*dz*lam)/(exp(2*lam)+exp(-2*lam));
      a3:=Sech(2*lam)*s*((exp(2*s)*s*cosh((1+j*dz)*lam)-(exp(s*(1+j*dz))-1)*lam*sinh(2*lam)-
         exp(2*s)*s*(1+2*lam*sinh((1+j*dz)*lam))));
      a:=(a1+a2+a3)/s/s;
       b:=a/lam/(s*s-lam*lam);
       c:=b/(2*lam-tanh(2*lam));
      u0[j]:=c.real;
      v0[j]:=c.imaginary;
    end;

Уже и разбила выражение на кусочки. Ругается на строке с a2: Project Project1.exe raised exception class EZeroDivide with message 'Division by zero'. Process stoped.
На строке с а1 тоже ругался, пока не убрала часть выражения в скобки {} (сокращается).
Delphi 6v. может стоит на 7 попробовать или все-таки дело в коде?

[Grag]

А для чего существует трассировка программы???

может стоит на 7 попробовать

Одна женщина сможет родить ребенка за девять месяцев, но девять женщин не смогут родить ребенка за один месяц...

[Эренка]

Сообщение от Grag

А для чего существует трассировка программы???

Плохо понимаю, что это и как ей пользоваться.
Другая похожая программа работает через раз, то ноль, то не ноль.
Прочитала по форуму, ошибка может быть из-за переполнения, но как ее решить с комплексными числами..

[Grag]

Обрабатывай исключительные ситуации блоками
try
.......
except
........
end;

[Эренка]

Так?

procedure calc_v0;
  var
    a,b,c,lam,i,a1,a2,a3:variant;
    j:integer;
  begin
  a:=VarComplexCreate; ...
  i.Imaginary:=1;
  i.real:=0;
  lam:=VarComplexCreate;
  lam:=sqrt(sqrt(0.25*Ta))*(1+i);
  for j:=0 to n2 do
    begin
      a1:={Sech(2*lam)*(}lam*lam+2*exp(s*(1+j*dz))*s*lam*lam+exp(2*s)*(s*s-lam*lam){)*cosh(2*lam)};
      try  
        try
      a2:=2*((1+2*s)*lam*lam+exp(2*s)*(s*s-lam*lam))*Cosh(lam-j*dz*lam)/(exp(2*lam)+exp(-2*lam));
      Except
         on E : Exception do
         ShowMessage(E.ClassName+' ошибка с сообщением : '+E.Message);
        end;
       a3:=Sech(2*lam)*s*((exp(2*s)*s*cosh((1+j*dz)*lam)-(exp(s*(1+j*dz))-1)*lam*sinh(2*lam)-exp(2*s)*s*(1+2*lam*sinh((1+j*dz)*lam))));
      a:=(a1+a2+a3)/s/s;
       b:=a/lam/(s*s-lam*lam);
       c:=b/(2*lam-tanh(2*lam));
      u0[j]:=c.real;
      v0[j]:=c.imaginary;    
      end;
end;

всё равно на a2 спотыкается и останавливает расчет

[Stilet]

В таком случае ставь точку прерывания на эту строку, когда до нее доползет выполнение проверь что из операндов дает ошибку.
Например выдели мышой (exp(2*lam)+exp(-2*lam)) и на выделенное наведи курсор - должно показать значение. Если 0 то проверяй формулу

[Anatole]

Наиболее вероятно, что функция exp, когда ей передать комплексный аргумент, возвращает 0

[Эренка]

Сообщение от Эренка

Delphi 6v. может стоит на 7 попробовать или все-таки дело в коде?

Проблема решается сносом Delphi 6 и установкой Delphi 7, может как-то и по другому, всем спасибо :-)