## 65 Reputation

6 years, 143 days

## @Kitonum  very nice , i want to kno...

@Kitonum  very nice , i want to know that what is the meaning of [4,2,1,1]

## thanks...

@tomleslie

sorry,

I have not mentioned  G[r] in my codes .

I am attaching my codes

restart:
with(DETools):
a1:=m^2*((1-N)/(2-N))*(r/2)*Dp:
a2:=((1-N)*m^2/(2-N))*((G[r]*(Nb-Nt)/64)*(r^5/6-h^2*r/2)-(B[r]/4)*(r^3/4-h^2*r/2)*(Nt/Nb)):
a3:=a1-a2:
a4:=r^2*a3:
DE1:=r^2*diff(v(r),r,r)+r*diff(v(r),r)-(m^2*r^2+1)*v(r)=a4:
b1:=dsolve(DE1,v(r));
simplify(collect(coeff(rhs(collect(expand(b1), [B[r], Nt])), B[r]),r))*B[r];

## @vv  Atleast (Nb-Nt) and (Nb/Nt)&n...

Atleast (Nb-Nt) and (Nb/Nt)  should be common

## @Kitonum  Thanks  alot, one mo...

@Kitonum  Thanks  alot, one more thing can we  put legend inside the plot, i have seen  it is in the bottom , left , so on. I have seen your post " https://www.mapleprimes.com/questions/222181-Insering-Legend-Inside-The-Figure-Plot-",  whose codes are  given below. Is it possible  in text plot C  that can we put  line, Dash dot, and longdash  , circle, diamond

with(plots):
A:=plot([sin(x), cos(x)], x=-Pi..3*Pi, color=[red,blue]):
B:=textplot([[3.7,0.8,sin(x)], [-1.8,0.5,cos(x)]], font=[times,bold,
14]):
#C:=textplot([[3,1,_, color=red], [-2.5,0.7,__,color=blue]], font=[times,bold,24]):
display(A, B, C, scaling=constrained);
# Note  C:=textplot([[3,1,longdash, color=red], [-2.5,0.7,dash dot,color=blue]], font=[times,bold,24]):

## @acer  apology  from next time...

@acer  apology  from next time i will mention version , curretly i am using maple13

## @Kitonum  thanks but i am getting t...

@Kitonum  thanks but i am getting the out put like this, i,e unable to evaluate the value of Pi

1-cos(.1*Pi)^2

## @Kitonum  thanks a lot, i have seen...

@Kitonum  thanks a lot, i have seen in others software line along with symbol can be drawn in plots(i.e line decorated with symbols {line +diamond, line with astrisk,... }), is it possible in maple or not?

## @Kitonum  thanks for your suggestio...

@Kitonum  thanks for your suggestions, is it possible in the above maple code out of 9 curves in a single plot can we change the style of lines as (line,dash, long dash for three lines and remaining 6 lines can we put the symbols as circle , cross, diamond so on)

i means lines and symbols both

## @Kitonum thanks for your help...

@Kitonum

thanks for your help

## @Kitonum thanks for your prompt rep...

@Kitonum thanks for your prompt reply. suppose if we have delta2=delta2 *beta and beta values are 0.1,0.2,0.3 then how we can plot the curves ,for beta=0.1,0.2,0.3  nine curves with in one graphs  for gamma=10, similarly for beta=0.1,0.2,0.3  nine curves with in one graphs for gamma=20, similary for gamma =30, and 40

```
* i have changed h:=z->1-(beta*delta2/2)*(1 + cos(2*(Pi/L1)*(z - L1))):
restart;
h:=z->1-(beta*delta2/2)*(1 + cos(2*(Pi/L1)*(z - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):
L1:=0.2:
alpha:=Pi/6:
plots:-display(Matrix(2,2, [seq(plot([seq(evalf(eval(lambda,Nb=j,beta=m1)),j=[0.1,0.2,0.3]),m1=[0.1,0.2,0.3])],delta2=0.02..0.1, color=[red,blue,green], legend=[Nb=0.1,Nb=0.2,Nb=0.3], labels=[typeset(`&delta;1`), typeset(conjugate(`&Delta;p`))], title=typeset("Effect of ", 'alpha', " when ", Gamma,"2=",Gamma2)), Gamma2=[10,20,30,40])]));```

## @acer  thanks , but can it be ...

thanks , but can it be possible for individual graphs subtitle can be write alongwith one title on top of all graph as your code is mention above

## @acer  not like that please see the...

@acer  not like that please see the sample  , I want title should by effect of  so on.... which should cover all the three graphs and there should not be any change on sub title on each graphs.

Note: Only on the top of all three figures one title

## @acer  Thanks  for  the&n...

@acer  Thanks  for  the  codes given below  i am unable to get the output,

can you help me to get output in such a way that  the text  should be display on all graph  .

here is my codes  and sample display

restart:
with(plots):

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=unapply(Int(K1,z=0..1), Gamma2):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
P1:=( Vector[row]
( [ seq
( plot
( [ seq
( eval(lambda(Gamma2), Nb=j),
j in [0.1,0.2,0.3]
)
],
delta2=0.02..0.1,
legend=[Nb=0.1,Nb=0.2,Nb=0.3],
labels=[typeset(`&delta;1`), typeset(conjugate(`&Delta;p`))],
title=typeset("Effect of ", ''alpha'', " when ", Gamma,"2=", Gamma2)
),
Gamma2 in [10,20,30]
)
]
)
);
P2:= textplot([0,1.2,typeset("effect of ",alpha=3,",",gamma=5," in the flow")]):
display(P1,P2);

sample graph

## @tomleslie  Thanks for the suggest...

Thanks for the suggestions,

if we change the values of H=0.1, 0.2,0.3  we would get three graphs one after the other, i am attaching the maple codes,  and sample graphs also .

instead of applying individual values of H is it possible to apply all H values in some procedure such that if we run the programme all the graphs should be display one after the other.

my_codes.mw

 1 2 3 4 5 6 7 Page 1 of 7
﻿