in this problem i used slip coundary condition.

the plot want to starts in 0.2, but the plot starts in 0

bc := f(0) = 0, (D(f))(0) = ak*((D^2)(f))(0), (D(f))(N) = 1, g(0) = -de*((D^2)(f))(0), g(N) = 0, T(0) = 1+be*(D(T))(0), T(N) = 0:

How to type double drivative in BC.

ψ =-0.09,-0.07,-0.04,-0.01,0,0.01,0.04,0.07,0.09

these are the ψ values.

then X=

here we can take eta =0..2 and X=-15..15
using this relation how to plot streamlines for eta against X.

Code:
restart; with(plots); fcns := {T(eta), f(eta)}; ep := .1; M := 1; kp := .5; n := 1; ec := .1; pr := 1; s := .1; N := 5; sys := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))*(diff(f(eta), eta))+ep*ep+(M+1/kp)*(ep-(diff(f(eta), eta))) = 0, diff(T(eta), eta, eta)+pr*(f(eta)*(diff(T(eta), eta))-n*(diff(f(eta), eta))*T(eta))+pr*(ec*(diff(f(eta), eta, eta))*(diff(f(eta), eta, eta))+ec*(M+1/kp)*(diff(f(eta), eta))^2+s*T(eta)) = 0; bc := f(0) = 0, (D(f))(0) = 1, (D(f))(N) = ep, T(0) = 1, T(N) = 0; R := dsolve(eval({bc, sys}), numeric, method = bvp[midrich], abserr = 0.1e-9, output = operator); psi = [-0.9e-1, -0.7e-1, -0.4e-1, -0.1e-1, 0, 0.1e-1, 0.4e-1, 0.7e-1, 0.9e-1]; for i to 9 do X[i] = psi[i]/f(eta); print(plots:-contourplot(X[i](X, eta), eta = 0 .. N, X = 0 .. 6)) end do

I am not interested in the mapleprime group, so I want to delete my account in this group.

Recently, my questions are deleted. last month I ask 3 questions but that all questions are deleted without my knowledge. yesterday also i ask one question "how to plot3d graph" that also deleted. "https://www.mapleprimes.com/questions/225618-How-To-Plot-The-3d-Graph?sq=225618". In the mapleprime, I am following the ethics and rules. 

Recently, my post is deleted without my knowledge.

https://www.mapleprimes.com/posts/209617-How-To-Plot-Stream-Line-contour?sp=209617