We have the first order ODE
[tex] y'=4t \sqrt y,~y(0)=1, [/tex]
for which i have found the exact solution, namely a fourth order polynomial.
I want a numerical method to solve the problem exactly. This method has to be a fourth order method, since this implies that the local error vanishes.
Now we change the problem so it becomes
[tex] y'=4t \sqrt y - \lambda(y-(1+t^2)^2),~y(0)=a, [/tex]
and the question is: for which values of [itex]\lambda[/itex] and [itex]a[/itex] does a method that has the above mentioned property solve the new problem exactly.
Of course, the obvious case is [itex]\lambda=0[/itex] and [itex]a=1[/itex], because in this case the new problem reduces to the first problem.
My idea is that the solution must be a fourth order polynomial, since a fourth order numerical method has to solve the new problem exactly.
Although I want your view on this and a strategy to find the values of [itex]\lambda[/itex] and [itex]a[/itex] for which the new problem is solved exactly by a fourth order numerical method.