@tomleslie

sorry, 

I have not mentioned  G[r] in my codes .

I am attaching my codes 

restart:
with(DETools):
a1:=m^2*((1-N)/(2-N))*(r/2)*Dp:
a2:=((1-N)*m^2/(2-N))*((G[r]*(Nb-Nt)/64)*(r^5/6-h^2*r/2)-(B[r]/4)*(r^3/4-h^2*r/2)*(Nt/Nb)):
a3:=a1-a2:
a4:=r^2*a3:
DE1:=r^2*diff(v(r),r,r)+r*diff(v(r),r)-(m^2*r^2+1)*v(r)=a4:
b1:=dsolve(DE1,v(r));
simplify(collect(coeff(rhs(collect(expand(b1), [B[r], Nt])), B[r]),r))*B[r];

@vv 

Atleast (Nb-Nt) and (Nb/Nt)  should be common

@Kitonum  Thanks  alot, one more thing can we  put legend inside the plot, i have seen  it is in the bottom , left , so on. I have seen your post " https://www.mapleprimes.com/questions/222181-Insering-Legend-Inside-The-Figure-Plot-",  whose codes are  given below. Is it possible  in text plot C  that can we put  line, Dash dot, and longdash  , circle, diamond

with(plots):
A:=plot([sin(x), cos(x)], x=-Pi..3*Pi, color=[red,blue]):
B:=textplot([[3.7,0.8,sin(x)], [-1.8,0.5,cos(x)]], font=[times,bold,
14]):
#C:=textplot([[3,1,_, color=red], [-2.5,0.7,__,color=blue]], font=[times,bold,24]):
display(A, B, C, scaling=constrained);
# Note  C:=textplot([[3,1,longdash, color=red], [-2.5,0.7,dash dot,color=blue]], font=[times,bold,24]):

@acer  apology  from next time i will mention version , curretly i am using maple13

@Kitonum  thanks but i am getting the out put like this, i,e unable to evaluate the value of Pi

1-cos(.1*Pi)^2

@Kitonum  thanks a lot, i have seen in others software line along with symbol can be drawn in plots(i.e line decorated with symbols {line +diamond, line with astrisk,... }), is it possible in maple or not?

@Kitonum  thanks for your suggestions, is it possible in the above maple code out of 9 curves in a single plot can we change the style of lines as (line,dash, long dash for three lines and remaining 6 lines can we put the symbols as circle , cross, diamond so on) 

i means lines and symbols both

@Kitonum

thanks for your help

@Kitonum thanks for your prompt reply. suppose if we have delta2=delta2 *beta and beta values are 0.1,0.2,0.3 then how we can plot the curves ,for beta=0.1,0.2,0.3  nine curves with in one graphs  for gamma=10, similarly for beta=0.1,0.2,0.3  nine curves with in one graphs for gamma=20, similary for gamma =30, and 40 

 
* i have changed h:=z->1-(beta*delta2/2)*(1 + cos(2*(Pi/L1)*(z - L1))):
restart;
 h:=z->1-(beta*delta2/2)*(1 + cos(2*(Pi/L1)*(z - L1))):
 K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
 lambda:=Int(K1,z=0..1):
 L1:=0.2:
 alpha:=Pi/6:
plots:-display(Matrix(2,2, [seq(plot([seq(evalf(eval(lambda,Nb=j,beta=m1)),j=[0.1,0.2,0.3]),m1=[0.1,0.2,0.3])],delta2=0.02..0.1, color=[red,blue,green], legend=[Nb=0.1,Nb=0.2,Nb=0.3], labels=[typeset(`δ1`), typeset(conjugate(`Δp`))], title=typeset("Effect of ", 'alpha', " when ", Gamma,"2=",Gamma2)), Gamma2=[10,20,30,40])]));

@acer 

thanks , but can it be possible for individual graphs subtitle can be write alongwith one title on top of all graph as your code is mention above

@acer  not like that please see the sample  , I want title should by effect of  so on.... which should cover all the three graphs and there should not be any change on sub title on each graphs.

Note: Only on the top of all three figures one title 

@acer  Thanks  for  the  codes given below  i am unable to get the output,  

can you help me to get output in such a way that  the text  should be display on all graph  .

 here is my codes  and sample display

restart:
with(plots):

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=unapply(Int(K1,z=0..1), Gamma2):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
P1:=( Vector[row]
( [ seq
   ( plot
         ( [ seq
              ( eval(lambda(Gamma2), Nb=j),
                j in [0.1,0.2,0.3]
              )
            ],
            delta2=0.02..0.1,
            legend=[Nb=0.1,Nb=0.2,Nb=0.3],
            labels=[typeset(`δ1`), typeset(conjugate(`Δp`))],
            title=typeset("Effect of ", ''alpha'', " when ", Gamma,"2=", Gamma2)
          ),
          Gamma2 in [10,20,30]
        )
      ]
    )
  );
P2:= textplot([0,1.2,typeset("effect of ",alpha=3,",",gamma=5," in the flow")]):
display(P1,P2);
 

@tomleslie 

sample graph

@tomleslie 

Thanks for the suggestions,

if we change the values of H=0.1, 0.2,0.3  we would get three graphs one after the other, i am attaching the maple codes,  and sample graphs also .

instead of applying individual values of H is it possible to apply all H values in some procedure such that if we run the programme all the graphs should be display one after the other.

my_codes.mw

@acer  thanks , for your prompt reply,  in the above graph i have taken t=Pi/3, if we take t:=0,Pi/4,Pi/3; then how we will  plot for different values of t