The good option of the same formula,
dy1 = (concconstb * depth * (12.0 * concconsta * concconsta * concconsta * strainslope * strainslope * widthslope * horizcentslope
- 6.0 * concconsta * concconsta * concconstb * (-2.0 * concconstc * strainslope * widthslope * horizcentslope + 2 * strainslope * strain0 * widthslope * horizcentslope + strainslope * strainslope * (2.0 * width0 * horizcentslope
+ 2.0 * widthslope * horizcent0 + widthslope * horizcentslope * depth))
+ 2 * concconsta * concconstb * concconstb * (-6.0 * concconstc * strain0 * widthslope * horizcentslope - 3.0 * concconstc * strainslope * (2.0 * width0 * horizcentslope
+ 2.0 * widthslope * horizcent0 + widthslope * horizcentslope * depth) + 3.0 * strainslope * strain0 * (2 * width0 * horizcentslope + 2.0 * widthslope * horizcent0 + widthslope * horizcentslope * depth)
+ strainslope * strainslope * (3.0 * width0 * (2.0 * horizcent0 + horizcentslope * depth) + widthslope * depth * (3.0 * horizcent0 + 2.0 * horizcentslope * depth)))
+ concconstb * concconstb * concconstb * (6.0 * concconstc * strain0 * (2.0 * width0 * horizcentslope + 2.0 * widthslope * horizcent0 + widthslope * horizcentslope * depth)
+ 2.0 * concconstc * strainslope * (3.0 * width0 * (2.0 * horizcent0 + horizcentslope * depth) + widthslope * depth * (3.0 * horizcent0 + 2.0 * horizcentslope * depth))
- 2.0 * strainslope * strain0 * (3 * width0 * (2.0 * horizcent0 + horizcentslope * depth) + widthslope * depth * (3.0 * horizcent0 + 2.0 * horizcentslope * depth))
- strainslope * strainslope * depth * (width0 * (6.0 * horizcent0 + 4.0 * horizcentslope * depth) + widthslope * depth * (4.0 * horizcent0 + 3.0 * horizcentslope * depth))))
+ 12.0 * (concconsta * strainslope + concconstb * concconstc) * (concconstb * strain0 - concconsta * strainslope) * (concconstb * width0 - concconsta * widthslope) * (concconstb * horizcent0 - concconsta * horizcentslope)
* log(concconsta + concconstb * depth)) / (12.0 * concconstb * concconstb * concconstb * concconstb * concconstb);