(a + i*b) - (c + i*d) = ((a-c) + i*(b-d))
528529530531532533534535536537538
err *= EPS; if (b.mod() <= err) return x; g=d.div(b); g2=g.mul(g); h=g2.sub(f.div(b).mul(2.0)); sq = (h.mul(m).sub(g2).mul(m-1)).sqrt(); //sq=sqrt((m-1)*(m*h-g2)); gp=g.add(sq); gm=g.sub(sq); double abp=gp.mod();