@robjohn @quid @DanielFischer I have to calculate approximations of the solution with the method $y^{n+1}=y^n+h \cdot [\rho \cdot f(t^n,y^n)+(1-\rho) \cdot f(t^{n+1},y^{n+1})] , n=0,....,N-1 \\ y^0=y0$ for various values of $\rho$ and the errors for uniform partitions with N=64, 128, ...., 4096,8192 subintervals. Determine the value of the parameter $\rho \in [0,1)$ such that the method has maximum order, let $\rho=\rho_m$.
PS: The order accuracy can be computed numerically as follows.
Let $E(N_1)$ (respectively $E(N_2)$) the error of the numerical method for $N_1$ ( respectively $N_2$) su…