Question: unable to get the desire plot

Dear Maple users Help me to  get the desire graph for this codes. 

restart:
with(plots):
with(IntegrationTools):
h:=z->piecewise( z<=d+1,   1,
                 z<=d+4,   1-(delta/(2))(1 + cos(2(Pi)*(z - 1 - 1/2))),                                                           z<=d+6,   1 ):
w0:=(-c*h(z)^2/4)+(3/64)(b*c-4*a)*h(z)^4+(19/2304)*b(b-4*a)*h(z)^6:
w1:=(c/4)+(1/16)*(4*a-b*c)*h(z)^2:
w2:=(1/256)(4(b*c-4*a)-b*h(z)^2):
w3:=(1/2304)b(b-4*a):
a:=(x4*S*Gr)*sin(alpha)/(4*x1*x5):
b:=(1/Da)+(x3*M/(x1*(1+m^2))):
c:=(1/x1)*Dp:
Dp:=96*x1/((6-b*h(z)^2)h(z)^4)(F+(a*h(z)/24)-((11/6144)b(b-4*a)*h(z)^8)):
x1:=1/((1-phi1)^2.5*(1-phi2)^2.5):
x2:=(1-phi2)((1-phi1)+phi1*Rs1/Rf)+phi2(Rs2/Rf):
x3:=(shnf)/(sf):
x4:=(1-phi2)((1-phi1)+phi1(RBs1)/(RBf))+phi2*((RBs2)/(RBf)):
x5:=khnf/kf:
shnf:=sbf*((ss2+2*sbf-2*phi2*(sbf-ss2))/(ss2+2*sbf+phi2*(sbf-ss2))):
sbf:=sf*((ss1+2*sf-2*phi1*(sf-ss1))/(ss1+2*sf+phi1*(sf-ss1))):
khnf:=kbf*((ks2+2*kbf-2*phi2*(kbf-ks2))/(ks2+2*kbf+phi2*(kbf-ks2))):
kbf:=kf*((ks1+2*kf-2*phi1*(kf-ks1))/(ks1+2*kf+phi1*(kf-ks1))):
RBs1:=(8933*16.7*10^6):
RBf:=(1063*1.8*10^6):
RBs2:=6320*18*10^6:
kf:=0.492:
sbf:=6.67*10^(-1): ss2:=2.7*10^(-8):
sf:=6.67*10^(-1):ss1:=59.6*10^(6):
ks2:=76.5:kf:=0.492: ks1:=401:
phi1:=0.01: phi2:=0.02:alpha:=Pi/4:m:=0.5:Da:=0.1:Gr:=5:delta:=1:S:=0.5:  d:=1:
                                      
W1:=w0+w1*r^2+w2*r^4+w3*r^6:

by varing M =2,5,7 and r varies from 0 to 1 i want this type of graphs.  please see the sample graphs