Hellow, help required to remove the errors

 How to  obtain the pressure drop i am unable to get the output. I am uploading the file  and the equations

help_dp.mw

Help required to get the derire output  for the differential equation . I am attacing the codes and the sample output.

restart:
with(DETools):
a1:=m^2*((1-N)/(2-N))*(r/2)*Dp:
a2:=((1-N)*m^2/(2-N))*(((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)):

how to evaulate the value of R(z) for different values of z=0 to 1 with an interval of 0.1 and print ten values  in one column

R(z):= 1-cos^2*(Pi*z);
 

Greetings, any one help  me to remove the error  in plots display,here is the codes 

i want 2 rows and 2 columns graphs with 9 lines in each graphs. 

i am attaching the maple questions also  like this i want 

 https://www.mapleprimes.com/questions/224443-How-Do-I-Get-A-Legend-Only-Three-Curves--In-Maple

https://www.mapleprimes.com/questions/225078-How-To-Remove-Error-In-Plot-Display

restart;
 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:
 
 alpha:=Pi/6:
 with(plots):
 display
  ( Vector[row]
    ([seq [ seq
        ( plot
          ( [ seq
              ( eval(lambda(Gamma2), Nb=j),
                j in [0.1,0.2,0.3]
              )
           d1 in [0.1,0.2,0.5])] ],
            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,40]
        )
      ]
    )
  );

Greetings,

help 

i want to add text with greek letters and english letters on top of the graph with out using title in plot

restart:
with(plots):
p1:=plot(sin(2*x), x = -Pi .. Pi):
p2:= textplot(effect of alpha =3,gamma=5 in the flow):
display(p1,p2):