In this problem,how to plot various values of Nc. That is,

 Sh = 0.5, R[d] = 0.7, alpha = (1/2)*Pi, Nc = 0.5, H = 0.4, `ε` =0.4;

 Sh = 0.5, R[d] = 0.7, alpha = (1/2)*Pi, Nc = 1, H = 0.4, `ε` = 0.4;

 Sh = 0.5, R[d] = 0.7, alpha = (1/2)*Pi, Nc = 1.5, H = 0.4, `ε` = 0.4;

 Sh = 0.5, R[d] = 0.7, alpha = (1/2)*Pi, Nc = 2, H = 0.4, `ε` = 0.4

MSNL.mw

HTR.mw

In above problem, Additionally How to  plot  heat transfer rate  Q versus L^2  for distinct porosity parmeters(Sh) , using  heat transfer rate formula, Q = (q*L)/(k*A*T[b])=theta'(1).

using  [Sh = 0.1, L^2 = 0.1, Nr =0 .1, Ha =0 .1, Pe = 0.1],  [Sh = 0.3, L^2 = 0.3, Nr = 0.1, Ha = 0.1, Pe =0 .1],   [Sh = 0.5, L^2 =0 .5, Nr =0 .1, Ha = 0.1, Pe =0 .1].

Boundary condition of my problem is   f(0) =1  and  d(f)(1) = -(k[f]/k[nf])*Bi*f(1).

But i am getting error when solve this.How to clear this .

my attached file is,

HNF_1.mw

How to rectify this error.

NF-1.mw

Good evening,could you please help me to find the values of Phi[2],Phi[3],Phi[4],....

The right hand side of the for loop is the equation i need.Is this maple code correct.

rps1.mw