Plot 如何加快九次积分计算？

我想画一个定积分

myIntegral[x_] := NIntegrate[Sqrt[(a - b)^2 + (c - d)^2]/ (a - b), {a, 0, x}, {b, x, 1}, {c, 0, 1}, {d, 0, 1}]
Plot[myIntegral[x], {x, 0, 1}]

我使用的是Mathematica 11.3。当我评估代码时，Mathematica立即给出了一些关于数值积分收敛太慢的警告，然后它似乎陷入了计算中。 例如，我尝试过NIntegrate的一些选项，但没有成功。

---

是否有办法在可接受的时间（例如几分钟）内获得绘图结果？

可以通过

Integrate

获得

Assuming[0<x<1,
 Simplify[
   Integrate[Sqrt[(a-b)^2 + (c-d)^2]/(a-b), {a,0,x}, {b,x,1}, {c,0,1}, {d,0,1}]]]

---

这应该比在单个

NIntegrate

上绘制数千个问题要快得多。如果你将问题发布在“谢谢”，我复制了我的问题，你会得到更好的帮助

(-11*Sqrt[2] + 11*Sqrt[2+(-2+x)*x] + x*(-8-7*Sqrt[2+(-2+x)*x] + 13*Sqrt[1+x^2] +
 2*x*(6-3*Sqrt[2+(-2+x)*x] + x*(-4+2*x + Sqrt[2+(-2+x)*x] - Sqrt[1+x^2]))) +
 18*x^2*ArcCoth[Sqrt[2+(-2+x)*x]] + 9*ArcSinh[1] - 3*ArcSinh[1-x] - 3*ArcSinh[x] +
 6*Log[1-x] + 9*x*Log[-1+Sqrt[2+(-2+x)*x]] + 3*((-2+x)*Log[1+Sqrt[2+(-2+x)*x]] -
 4*x*Log[-(((-1+x)*(1 + Sqrt[1+x^2]))/x)] + x^3*(Log[2+x^2+2*Sqrt[1+x^2]] +
 2*Log[(1-x)/(x+x*Sqrt[2-2*x+x^2])])))/36