Advanced Search

Question: how I can plot a contour like attached figure

Question:how I can plot a contour like attached figure

9009134 110 Maple 2022
contour plotting duplicate-question labeling
August 11 2023
0

how I can plot phi[2] as a contour like attached figure?

tez-1.mw


 

> 

restart
> 

``
> 

beta := 2.5; lambda := 0.1e-1; b := Pi; a := Pi; alpha := 1; y[1] := 1.5; y[2] := 1.5; x[1] := -1; x[2] := 1; Q[1] := 40; Q[2] := 35

2.5

  

0.1e-1

  

Pi

  

Pi

  

1

  

1.5

  

1.5

  

-1

  

1

  

40

  

35

 (1)
> 

NULL
> 

NULL
> 

v := (2*n-1)*Pi/(2*b)

n-1/2

 (2)
> 

Delta := exp(2*v*a)*(alpha*v+beta)*(1+lambda)-(1-lambda)*(alpha*v-beta)

1.01*exp(2*(n-1/2)*Pi)*(n+2.000000000)-.99*n+2.970000000

 (3)
> 

g[22] := ((alpha*v+beta)*((1+lambda)*exp(-v*abs(x-xi))+(-1+lambda)*exp(-v*(x+xi)))*exp(2*v*a)+(alpha*v-beta)*((1+lambda)*exp(-v*(x+xi))+(-1+lambda)*exp(-v*abs(x-xi))))/(2*v*Delta)
> 

g[21] := ((alpha*v+beta)*exp(v*(2*a+xi))+(alpha*v-beta)*exp(-v*xi))*exp(-v*x)/(v*Delta)
> 

NULL
> 

u[2] := int(2*g[21]*Q[1]*Dirac(xi-x[1])*sin(n*Pi*y[1]/b)/b, xi = -a .. 0)+int(2*g[22]*Q[2]*Dirac(xi-x[2])*sin(n*Pi*y[2]/b)/b, xi = 0 .. infinity)
> 

NULL
> 

phi[2] := sum(u[2](x)*sin(v*y), n = 1 .. 30)
> 

NULL
> 

``
> 

plot3d(phi[2], x = 0 .. 5, y = 0 .. b)

 
> 

NULL


 

Download tez-1.mw


 

> 

restart
> 

``
> 

beta := 2.5; lambda := 0.1e-1; b := Pi; a := Pi; alpha := 1; y[1] := 1.5; y[2] := 1.5; x[1] := -1; x[2] := 1; Q[1] := 40; Q[2] := 35

2.5

  

0.1e-1

  

Pi

  

Pi

  

1

  

1.5

  

1.5

  

-1

  

1

  

40

  

35

 (1)
> 

NULL
> 

NULL
> 

v := (2*n-1)*Pi/(2*b)

n-1/2

 (2)
> 

Delta := exp(2*v*a)*(alpha*v+beta)*(1+lambda)-(1-lambda)*(alpha*v-beta)

1.01*exp(2*(n-1/2)*Pi)*(n+2.000000000)-.99*n+2.970000000

 (3)
> 

g[22] := ((alpha*v+beta)*((1+lambda)*exp(-v*abs(x-xi))+(-1+lambda)*exp(-v*(x+xi)))*exp(2*v*a)+(alpha*v-beta)*((1+lambda)*exp(-v*(x+xi))+(-1+lambda)*exp(-v*abs(x-xi))))/(2*v*Delta)
> 

g[21] := ((alpha*v+beta)*exp(v*(2*a+xi))+(alpha*v-beta)*exp(-v*xi))*exp(-v*x)/(v*Delta)
> 

NULL
> 

u[2] := int(2*g[21]*Q[1]*Dirac(xi-x[1])*sin(n*Pi*y[1]/b)/b, xi = -a .. 0)+int(2*g[22]*Q[2]*Dirac(xi-x[2])*sin(n*Pi*y[2]/b)/b, xi = 0 .. infinity)
> 

NULL
> 

phi[2] := sum(u[2](x)*sin(v*y), n = 1 .. 30)
> 

NULL
> 

``
> 

plot3d(phi[2], x = 0 .. 5, y = 0 .. b)

 
> 

NULL


 

Download tez-1.mw

 

 
Please Wait...