1.  解析解
據(jù) Acos(wx)cosh(wy)=-B；-Asin(wx)sinh(wy)=C
有
cosh(wy)cos(wx)=-B/A、；sinh(wy)sin(wx)=-C/A
cosh(wy)cos(wx) + i*sinh(wy)sin(wx) = -B/A + i* -C/A
cosh(wy+i*wx)=-B/A + i* -C/A=z
wy+i*wx = Arccosh(z)=Ln(z+sqrt(z^2-1))
=Ln|z+sqrt(z^2-1)|+i*(2k*pi + Arg(z+sqrt(z^2-1)))
2.  
1stop code
parameters  wx, wy;
Constan  A=…, B=…, C=…;   // 自己賦值
Function  A*cos(wx)*cosh(wy)=-B;
        -A*sin(wx)*sinh(wy)=C;