Euler-Maclaurin Summation

Using EM summation formula estimate $$ \sum_{k=1}^n \sqrt k $$

up to the term involving $\frac{1}{\sqrt n}$

My attempt is $$ \sum_{k=1}^n \sqrt k = \frac{2 \sqrt{n^3}}{3} -\frac{2}{3} + \frac 1 2 (\sqrt n -1)+ \frac{1}{24} (\frac{1}{\sqrt n} -1) + \int_1^n P_{2k+1}(x)f^{(2k+1)}(x)dx $$ I am not sure what can be said about the integral. Please tell me if I have made a mistake and how I can solve the integral? Have I stopped the summation at the right point?

Using the Euler-Maclaurin Sum Formula, we get $$ \sum_{k=1}^n\sqrt{k}=\tfrac23n^{3/2}+\tfrac12n^{1/2}+\zeta(-\tfrac12)+\tfrac1{24}n^{-1/2}-\tfrac1{1920}n^{-5/2}+\tfrac1{9216}n^{-9/2}+O(n^{-13/2}) $$ where the constant $\zeta(-\frac12)=-0.20788622497735456602$ is explained below.

By the Euler-Maclaurin Sum Formula, we have for $\mathrm{Re}(s)\gt-1$ $$ \sum_{k=1}^nk^{-s}=\frac1{1-s}n^{1-s}+\frac12n^{-s}+\zeta_\ast(s)+O(n^{-s-1})\tag{1} $$ Define the sequence of functions $$ \zeta_n(s)=\sum_{k=1}^nk^{-s}-\frac1{1-s}n^{1-s}-\frac12n^{-s}\tag{2} $$ For all $n\ge1$, $\zeta_n(s)$ is analytic. For $\mathrm{Re}(s)\gt1$, $\lim\limits_{n\to\infty}\zeta_n(s)=\zeta(s)$. Estimate $(1)$ says that on compact subsets of $\{s\in\mathbb{C}\setminus\{1\}:\mathrm{Re(s)\gt-1}\}$, $\zeta_n(s)$ converges uniformly. Thus, for $s\in\mathbb{C}\setminus\{1\}$ and $\mathrm{Re}(s)\gt-1$, $\lim\limits_{n\to\infty}\zeta_n(s)$ is analytic, so by analytic continuation, for $\mathrm{Re}(s)\gt-1$, $$ \zeta_\ast(s)=\lim\limits_{n\to\infty}\zeta_n(s)=\zeta(s)\tag{3} $$

The correct expansion is given by

$$\begin{align} \sum_{i=1}^n f(i) &= \int_1^n f(x)dx + B_1 [f(n)+f(1)]\\\\ & + \sum_{k=1}^m \frac{B_{2k}}{(2k)!} \left(f^{(2k-1)}(n)-f^{(2k-1)}(1)\right)\\\\ & +\frac{1}{(2m+1)!}\int_1^n P_{2m+1}(x)f^{(2m+1)}(x) \end{align}$$

For $f(x)=x^{1/2}$ and $m=1$, we have

$$\begin{align} \sum_{i=1}^n\sqrt{i}&=\frac23 (n^{3/2}-1)+\frac12 (n^{1/2}+1)+\frac{1}{24}(n^{-1/2}-1)+\frac{1}{3!}\int_1^n P_3(x)f^{(3)}(x)dx\\\\ &=\frac23x^{3/2}+\frac12n^{1/2}+\left(-\frac23+\frac12-\frac{1}{24}\right)+\frac{1}{24}n^{-1/2}+\frac{1}{3!}\int_1^n P_3(x)f^{(3)}(x)dx \end{align}$$

The next task is to determine an estimate for the remainder $R$ where $R$ is the integral

$$R\equiv \frac{1}{3!}\int_1^n P_3(x)f^{(3)}(x)dx$$

We can find an estimate for the integral term

$$\int_1^n P_{3}(x)f^{(3)}(x)dx \tag1$$

using the well-known expression for the Bernoulli Polynomial

$$B_3(x)=x^3-\frac32 x^2+\frac12 x$$

We can easily verify that

$$-\frac{1}{12\sqrt{3}}<B_3(x)<\frac{1}{12\sqrt{3}}$$

for $x\in [0,1]$. Since $P_{2k+1}(x)=B_{2k+1}(x-\lfloor x\rfloor)$ then we see immediately that for $1<x<n$

$$-\frac{1}{12\sqrt{3}}< P_{3}(x) < \frac{1}{12\sqrt{3}}\tag 2$$

Now, we note that for $f(x)=x^{1/2}$, $f^{3}(x)=\frac38 x^{-5/2}>0$ for $x\in [1,n]$. Now, we use $(2)$ in the estimate of $(1)$ to reveal

$$\begin{align} -\frac{1}{48\sqrt{3}}\left(1-n^{-3/2}\right)<\int_1^n P_{3}(x)f^{(3)}(x)dx \le \frac{1}{48\sqrt{3}}\left(1-n^{-3/2}\right) \end{align}$$

Finally, we have the upper and lower bounds

$$\bbox[5px,border:2px solid #C0A000]{\sum_{k=1}^n\sqrt{k}\le \frac23n^{3/2}+\frac12n^{1/2}+\left(-\frac{5}{24}+\frac{1}{288\sqrt{3}}\right)+\frac{1}{24}n^{-1/2}-\frac{1}{288\sqrt{3}}n^{-3/2}}$$

$$\bbox[5px,border:2px solid #C0A000]{\sum_{k=1}^n\sqrt{k}\ge \frac23n^{3/2}+\frac12n^{1/2}+\left(-\frac{5}{24}-\frac{1}{288\sqrt{3}}\right)+\frac{1}{24}n^{-1/2}+\frac{1}{288\sqrt{3}}n^{-3/2}}$$