Dear, I want to plot the graph for differerent values of F=10,20,30 one after the other  in the title automatic it should take the values of F and shown in the following sample here is my codes,  and sample graph

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);
 

Dear sir I want to plot the graph for different values of F=0.1, F=0.2,F=0.5 . My codes are given and i am attaching the sample graph. For F=0.1, Nb=0.1,0.2,0.3 three curves

Similarly For F=0.2,Nb=0.1,0.2,0.3 three curves and for F=0.5, Nb=0.1,0.2,0.3 three curves. Total nine curves 

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):

L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);

I want to plot contour for piecewise  function , here is my code

with(plots):
d1:=0.1:eta:=0.01:k:=2:L:=1:
mu:=0.01:delta1:=0.1:delta2:=0.2:
d1:=0.2:d2:=0.2:L1:=0.2: L2:=0.2:
h:=z->piecewise( z<=d1,    1,
                 z<=d1+L1,   1-(delta1/(2))*(1 + cos(2*(Pi/L1)*(z - d1 - L1/2))), 
                        z<=d1+L1+d2,  1 ,          
                    z<=d1+L1+d2+L2,  1-(delta2/(2))*(1 + cos(2*(Pi/L2)*(z - d1 - d2-L1-L2/2))),
                 z<=L,    1):
A1:=((k+2)*(k+3)*(1-mu)^(k+1)-2*(1-mu)^(k+2)*(k+2+mu))^(1/k):
A2:=(A1*r/h(z)^(1+3/k)):
contourplot(A2,z=0..1,r=-h(z)..h(z),colour=black,axes=boxed);

in plot line style should to like  and .so on....

line & diamond   ,   line  and circle   ,  line and cross,    line and cross,

 same should be displayed in legends

dear sir i want to plot a graph for different values of h(z) by applying do loop, hear is my codes

In K1 it should take first h(z)

In K2 it should take second h(z)

IN K3 it should take third h(z) 

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L2)*(z - d2 - L2))):
h:=z->1+(delta2/2):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
K2:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
K3:=((4/h(z)^4)+(cos(alpha)/F)+h(z)^2+Nb*h(z)^4):
lambda1:=Int(K1,z=0..0.2):
lambda2:=Int(K2,z=0.2..0.4):
lambda3:=Int(K3,z=0.4..0.6):
lambda:=(lmbda1+lambda2+lambda3):
> 
F:=0.3:
L1:=0.2:
d1:=0.2:
d2:=0.2:
L2:=0.3:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);