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);
 

Dear sir,

 i want to plot a graph with  tick marks on all sides , i am attaching codes and  a sample graph, see the attachmnent,

The graph range should be displace only on x and y axis and the vaues should not be on opposite to x and y axis

please see the 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):
F:=0.3:
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);

hai,

 i want to plot a graph here is my codes and i want the graph like this, can any body help me

restart:
with(DEtools):
with(plots):

h(z):=1-(delta2/2)*(1 + Cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h^4)-(sin(alpha)/F)-h^2+Nb*h^4):
lambda:=int(K1,z=0..1):
F:=0.3:
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
evalf(lambda);

Dear  i want to plot the graphs for different values of R[2,5,6.5]   here is the codes: Either ido loop or array i have to use or other procedure 

The graph display should be at R=2, one graph, R=5 another and R=6.5 

onemore graph

restart:
with(DEtools):
eta:=1+k*x+epsilon*sin(2*Pi*x):
A1(x):=-exp(-alpha*x)*J^2*R/(2*eta^3+6*xi*J^2*eta^2):

psi0:=A1(x)*y^3:
psi:=delta*psi0:
V:=-diff(psi,x):
delta:=0.1:
epsilon:=0.01:
alpha:=1:
xi:=0.001:
k:=0.1:
As:= [0, 2, 4]:
x:=0.2
plot([seq]([eval(V, J= A), y, y= 0..eval(eta, J= A)], A= As),title=" velocity at R=2.0 ",labels= ["v", "y"],color= [green, red, blue], linestyle= [solid, dash, dot],legend= [seq](J = A, A= As), axes=boxed);