Question:Apply Keller Box method for ODE

Question:Apply Keller Box method for ODE

Maple 2018

Dear sir ,

I have implemented Dsolve method the code was executed, but i need to apply Kellor Box method to solve the ODES

Please can any one help how to implement?

because there is no post regarding the Kellor box method.

 >
 >
 >
 >
 >
 >
 >
 >
 > with(plots):  cols := [red, blue, black,green]:  plotA:= display   ( [ seq       ( odeplot         ( Ans[k],[eta,(f(eta))],           eta=0..1,           color=cols[k]         ),         k=1..numelems(KpVals)       )     ],linestyle = "solid",     'axes'= 'boxed',labels=[eta,'f(eta)'],labelfont=[TIMES,BOLD,16]   );
 > with(plots):   cols := [red, blue, black,green]: plotB:= display( [ seq( odeplot         ( Ans[k],[eta,Theta(eta)],           eta=0..1,           color=cols[k]         ),         k=1..numelems(KpVals)       )     ],linestyle = "solid",     'axes'= 'boxed',labels=[eta,'Phi(eta)'],labelfont=[TIMES,BOLD,16]   );
 > with(plots):   cols := [red, blue, black,green]: plotC:= display( [ seq( odeplot         ( Ans[k],[eta,Phi(eta)],           eta=0..1,           color=cols[k]         ),         k=1..numelems(KpVals)       )     ],linestyle = "solid",     'axes'= 'boxed',labels=[eta,'Phi(eta)'],labelfont=[TIMES,BOLD,16]   );
 > with(plots):  cols := [red, blue, black,green]:  plotA:= display   ( [ seq       ( odeplot         ( Ans[k],[eta,(diff(f(eta),eta))],           eta=0..1,           color=cols[k]         ),         k=1..numelems(KpVals)       )     ],linestyle = "solid",     'axes'= 'boxed',labels=[eta,"f '(eta)"],labelfont=[TIMES,BOLD,16]   );
 >
 >