Вычислить значение величины, используя вычисление гиперболического котангенса по формуле - Pascal

uses
  sysutils, math;
 
function MyPower(base, exponent : double) : double;
{(c) volvo}
type
  TFrac =
    record
      num : longint; // chisl
      denom : longint; // znam
    end;
 
function Dec2Frac(f : double) : TFrac;
var
  df : Double;
  sign : boolean;
begin
    sign := f < 0;
    f := Abs(f);
    result.num := 1;
    result.denom := 1;
 
    df := result.num / result.denom;
    while (df <> f) do
    begin
      if (df < f) then
         inc(result.num)
      else
      begin
         inc(result.denom);
         result.num := trunc(f * result.denom);
      end;
      df := result.num / result.denom;
    end;
 
    if sign then result.num := - result.num;
  end;
 
var
  f : TFrac;
  sign : boolean;
  p : double;
begin
  f := dec2frac(exponent);
  if ((base < 0.0) and not odd(f.denom)) or
     ((base = 0.0) and (f.num < 0))
  then
    raise Exception.Create('invalid argument(s)')
  else
  begin
    sign := base < 0;
    p := power(abs(base), abs(f.num));
    p := power(p, 1 / f.denom);
    if sign then p := -p;
 
    if f.num < 0 then p := 1 / p;
    result := p;
  end;
end;
 
function cth(x:double):double;
begin
  cth:=(exp(x)+exp(-x))/(exp(x)-exp(-x));
end;
 
var
  x,y,chisl,znam:double;
begin
  readln(x);
  chisl:=4+cth(2*x-1)-cth(x+4);
  znam:=7.85*cth(4*x*x*x-2*x+1);
  znam:=MyPower(znam,1/3);
  y:=chisl/znam;
  writeln(y:10:3);
  readln;
end.