Can anybody help me to code boundary conditions in MATLAB for Keller Box Method?

f^’=1,f=S,θ^’=-r_1 [1+θ],ϕ^’=-r_2 [1+ϕ] at η=0
f^’=0,f^”=0,θ=0,ϕ=0 as η→∞Can anybody help me to code boundary conditions in MATLAB for Keller Box Method?

f^’=1,f=S,θ^’=-r_1 [1+θ],ϕ^’=-r_2 [1+ϕ] at η=0
f^’=0,f^”=0,θ=0,ϕ=0 as η→∞ Can anybody help me to code boundary conditions in MATLAB for Keller Box Method?

f^’=1,f=S,θ^’=-r_1 [1+θ],ϕ^’=-r_2 [1+ϕ] at η=0
f^’=0,f^”=0,θ=0,ϕ=0 as η→∞ keller box method, differential equations MATLAB Answers — New Questions

​