Dear Maple user, please help me to show two legends inside the plot. Here is the codes attached. 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):
T1:=Int(K1,z=0..1):
L1:=0.2: F:=10:
d1:=0.2:
T2:=evalf(T1):
T3:=unapply(T2,alpha,Nb,delta2):
plot([seq(seq(T3(alpha, Nb, delta2), Nb in [0.1, 0.2,0.3]), alpha in [Pi/6, Pi/4, Pi/3])], delta2 = 0.02.. 0.1,  titlefont = ["ROMAN", 15], labels = ["δ1", "T3"], labeldirections = ["horizontal", "vertical"], labelfont = ["SYMBOL", 11], linestyle = [solid, solid,solid,longdash, longdash,longdash, spacedash,spacedash, spacedash], symbol = [BOX, CROSS, CIRCLE, BOX, CROSS, CIRCLE, BOX, CROSS, CIRCLE], color = [blue,red, green, blue,red, green,blue,red, green], axes = boxed,size = [1000, 1000]);

Dear Maple user i want to extract the data from the given graph and store in excel file. where the first column contain the value of lambda in that substitude the values of delta2 ranging from 0.002 to 0.1 (atleast 20 values) and second column  contain Nb =0.1, third column Nb= 0.2 and third column Nb=0.3. Thanks in advance

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: F:=10:
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 Maple user

I want to add caption for all four graphs with a,b,c,d at the bottom centre of each graph that should be arrange in 2X2 array and the array bottom caption as Figure 6, also i want to save the output of the array graphs in desktop with jpeg format. sample graph i have uploaded.

restart:
with(plots):
A1:= plot(sin(x), x = -Pi .. Pi,title = "Graph at M^2=1",caption = "Fig a"):# This caption is at the centre bottom of first graph and the title of A1 is M*M (M^2)should be display
A2:= plot(cos(x), x = -Pi .. Pi,caption = "Fig b"):# This caption is at the centre bottom of second graph
A3 := plot3d(1.5^x*sin(y), x = -3 .. 3, y = -3 .. 3,caption = "Fig c"):# This caption is at the centre bottom of Third graph
A4 := plot3d(1.3^x*sin(y), x = -1 .. 2*Pi, y = 0 .. Pi,caption = "Fig d"):# This caption is at the centre bottom of forth graph
display(Array(1 .. 2, 1 .. 2, [[A1, A2], [A3, A4]]), aligncolumns = [1],caption = "Fig 6"):# This caption at the bottom centre of array graphs 

Dear Maple user I am facing error while running the codes  to plot the graph for two data sets .

I am attaching the files.

Error_in_Display_1.mw

Dear maple user i have a partial differential equation whose first approximation solution is w(r,z)=(r^2-h^2)/4 , i want to calculate the two approximations. I am attaching the file please see  

# I dont want to solve directly.

HP_Help.mw