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