Maple 2022 Questions and Posts

These are Posts and Questions associated with the product, Maple 2022

I wanted to solve a pde of two variables in spherical coordinates with a boundary condition at a specific radius. I am doing something wrong because the pdsolve just evaluates for several minutes and does not actually give anything.

Could someone have a look and let me know what I am doing wrong.

restart

``

with(Physics[Vectors])

Setup(mathematicalnotation = true)

[mathematicalnotation = true]

(1)

eq3 := Laplacian(Phi_out(r, theta)) = 0

((r^2*(diff(diff(Phi_out(r, theta), r), r))+2*r*(diff(Phi_out(r, theta), r))+diff(diff(Phi_out(r, theta), theta), theta))*sin(theta)+(diff(Phi_out(r, theta), theta))*cos(theta))/(r^2*sin(theta)) = 0

(2)

bc1 := eval(Phi_out(r, theta)-Omega*mu*(-1+cos(2*theta))/(2*r), r = R) = 0

Phi_out(R, theta)-(1/2)*Omega*mu*(-1+cos(2*theta))/R = 0

(3)

``

pdsolve([eq3, bc1])

NULL

Download laplacian.mw

how to make a phaseportrait curves with animation? so i can see how the curve going on. Thank you

caring_phase_pemanenan_predator_dgn_parametersesuai_jurnal.mw

Hi everyone, I am trying to draw horizontal and vertical line on the graph. Command is working but vertical line is not displayed. Looking for help in this regrads.

vertical_and_horizonatl_line_help.mw

Hi everyone, I am trying to draw bar graph for my expression, but received error. Could anyone please help me in this regard.

Help_Bar_graph.mw

Hi

On the latest Maple 2022 version, using dsolve with 1st order ODE and Dirac function returns an incorrect solution

dsolve({D(y)(t) + y(t)/tau = Dirac(t)/tau, y(0) = 0})

See attached image. The correct impulse response should be :

y(t) = exp(-t/tau)*Heaviside(t)/tau

The returned solution is 0 for t = 0 and scaled by a 1/2 for t > 0. I never had such an issue with older Maple versions. Is this a bug or am I doing something wrong

Here is the expected solution followed by the Maple solution:

acer helped me out last week with converting strings to dates which was very helpful. However now I stumble on 2 other problems:

- I do not succeed to plot DS6 (time difference versus a numeric value) and to not understand why

- in the cvs data sheet the length of the columns is not equal. The empty cells are automatical filled with "". This prevents me to use numelems as it will also take into consideration the cells filled with "". Is there a way to prevent this

As usual I would like to thank you all, power users, for your patience and help

Test.mw

Hi Everyone; I have plotted two different expressions separately and then shown them in a single graph. From the final display, we can observe that the peak of one expression, A is larger than that of another expression B. Now, the question is, how do I calculate the decrease or increase in the percentage of both curves? This means the peak of expression A is 20%, 8%, or 7% decreases or increases as of expression A.

Help_percentage_inccrease_or_decrease.mw

more explanation, each curve has its maximum point at some value of x, need to calculate that point for each curve. Then, we combined all curves we observed that the peak of each curve is smaller than others, so I am interested to observing that difference in percentage at the maximum value of x.

for reference, I uploading a file of which idea I want to implement for my problem. This is not my problem, but I want to implement it like this.

For_reference.mw

This is most likely a simple question for the power users of this forum, but I do not manage to find a solution. I have date in an Excel file. The first column consists of dates (mm/dd/yyyy) and the second of time (hh:mm). I can easily concatenate both in Maple using cat("9/7/2023", "10:22") but how can I convert the obtained string into a date+time that Maple understands? 

Thank you in advance for your help.

I create a system of equations (with 9 linear equations and 7 variables).
I get 7 equations from the multiplication of a matrix M with a transposed vector of the 7 variables equals the transposed vector of the 7 variables. Other two simple equations are necessary because they are restrictions. Those two equations a re very simple and thave the 7 variable sin it.
To start I have been trying with fsolve but i haven't been able to solve it yet, as I also get the error: "Error, (in fsolve) number of equations, 9, does not match number of variables, 7"
Have you andy idea to solve this problem? Thanks in advance.

Hi all, any one help  me to find the values of constants by using given condition and then how to varify that the goiven condition varify the expression. I have found manually but want to varify through maple.

help.mw

HPM.mw

  • Please help me
  • How do I write Homotopy Perturbation method for a partial differential equation for this question in Maple?

Download HPM.mw

 

Using the SingularValues command with any output type specified causes this error

Error, (in sanity) invalid input: LinearAlgebra:-SingularValues expects value for keyword parameter output to be of type {list(identical(U,S,Vt,list)), identical(U,S,Vt,list)}, but received [U, S, Vt]

Below is a minimal example to trigger the error. Is there something wrong about how I am calling the function?

restart:
with(LinearAlgebra):

sanity := proc()
  local A,U,S,Vt;
  A :=RandomMatrix(3,10);
  U,S,Vt := SingularValues(A, output=['U','S','Vt']);
end proc:

sanity();

Good day all,

please I am new Maple user, i used the same code to solve boundary layer flow. Itried to use the same code in solving nanofluid but the solution is not coverging at some point. But the author i want to compared my results with had solution at those points.

Is there any way i can get solution of this model at this with these values?. Please i need your help. Thank you in advance

secod_code.mw

NULL

restart

with*plots; ode1 := (1+1/bet)*(diff(f(eta), eta, eta, eta))-(diff(f(eta), eta))^2+f(eta)*(diff(f(eta), eta, eta))+Gr*(theta(eta)+lambda1*theta(eta)^2)+Gc*(phi(eta)+lambda2*phi(eta)^2)-M*(diff(f(eta), eta)) = 0

(1+1/bet)*(diff(diff(diff(f(eta), eta), eta), eta))-(diff(f(eta), eta))^2+f(eta)*(diff(diff(f(eta), eta), eta))+Gr*(theta(eta)+lambda1*theta(eta)^2)+Gc*(phi(eta)+lambda2*phi(eta)^2)-M*(diff(f(eta), eta)) = 0

(1)

ode2 := (diff(theta(eta), eta, eta))*(1+(4*R*(1/3))*(theta(eta)+CT)^3)/Pr+4*R*(theta(eta)+CT)^2*(diff(theta(eta), eta))^2+f(eta)*(diff(theta(eta), eta))+Nb*(diff(phi(eta), eta))*(diff(theta(eta), eta))+Nt*(diff(theta(eta), eta))^2-de*(f(eta)*(diff(f(eta), eta))*(diff(theta(eta), eta))+f(eta)^2*(diff(theta(eta), eta, eta))) = 0

(diff(diff(theta(eta), eta), eta))*(1+(4/3)*R*(theta(eta)+CT)^3)/Pr+4*R*(theta(eta)+CT)^2*(diff(theta(eta), eta))^2+f(eta)*(diff(theta(eta), eta))+Nb*(diff(phi(eta), eta))*(diff(theta(eta), eta))+Nt*(diff(theta(eta), eta))^2-de*(f(eta)*(diff(f(eta), eta))*(diff(theta(eta), eta))+f(eta)^2*(diff(diff(theta(eta), eta), eta))) = 0

(2)

ode3 := diff(phi(eta), eta, eta)+Pr*Le*(diff(phi(eta), eta))*f(eta)+Nt*(diff(theta(eta), eta, eta))/Nb-Pr*Le*dc*(f(eta)*(diff(f(eta), eta))*(diff(phi(eta), eta))+f(eta)^2*(diff(phi(eta), eta, eta))) = 0

diff(diff(phi(eta), eta), eta)+Pr*Le*(diff(phi(eta), eta))*f(eta)+Nt*(diff(diff(theta(eta), eta), eta))/Nb-Pr*Le*dc*(f(eta)*(diff(f(eta), eta))*(diff(phi(eta), eta))+f(eta)^2*(diff(diff(phi(eta), eta), eta))) = 0

(3)

bcs1 := f(0) = 0, (D(f))(0) = 1+xi*((D^2)(f))(0), (D(f))(6) = 0, theta(6) = 0, phi(0) = 1+`ς`*(D(phi))(0), theta(0) = 1+Zeta*(D(theta))(0), phi(6) = 0

f(0) = 0, (D(f))(0) = 1+xi*(D(f))(0)^2, (D(f))(6) = 0, theta(6) = 0, phi(0) = 1+varsigma*(D(phi))(0), theta(0) = 1+Zeta*(D(theta))(0), phi(6) = 0

(4)

NULL

 

L := [0., .1, .4, 1]

[0., .1, .4, 1]

(5)

 
  

NULL


  for k to 10 do
      sol_All := dsolve
                 ( eval
                   ( {bcs1, ode1, ode2,ode3},
                     [xi= L[k],R=1, CT=0.2, Gr=0.5, Gc=0.5, lambda1=0.2, lambda2=0.2, Pr=0.72, Nb=0.5,  Le=1, dc=0.2, de=0.2, Nt=0.5, Zeta=0.5, varsigma=0.5, M=1.5, bet=1.5]
                   ),
                   [f(eta), theta(eta),phi(eta)],
                   numeric,
                   output = listprocedure
                 );
      Y_sol || k := rhs(sol_All[5]);
      YP_sol || k := -rhs(sol_All[6]);

feta || k := rhs(sol_All[4]);
      fpeta || k := rhs(sol_All[3]);

phieta || k := rhs(sol_All[8]);
      phipeta || k := rhs(sol_All[7])
  end do:

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

 

for k to 10 do L[k], [(Y_sol || k)(0), (YP_sol || k)(0)] end do

0., [Y_sol1(0), YP_sol1(0)]

 

.1, [Y_sol2(0), YP_sol2(0)]

 

.4, [Y_sol3(0), YP_sol3(0)]

 

1, [Y_sol4(0), YP_sol4(0)]

 

Error, invalid subscript selector

 

for k to 10 do L[k], [(feta || k)(0)] end do

0., [feta1(0)]

 

.1, [feta2(0)]

 

.4, [feta3(0)]

 

1, [feta4(0)]

 

Error, invalid subscript selector

 

NULL

Error, mismatched or missing bracket/operator

"for k  from 1 to  10 do L[k], ([phieta ||k(0),phipeta||k(0)]); od; "

 

NULL

  plot( [ seq((Y_sol||j)(eta), j = 1..6)],
         eta = 0 .. 6,
         labels = [eta, theta(eta)],
         axes = boxed
      );
plot( [ seq((YP_sol||j)(eta), j = 1..6)],
         eta = 0 .. 8,
         labels = [eta, thetaprime(eta)],
         axes = boxed
      );

 plot( [ seq((feta||j)(eta), j = 1..6)],
         eta = 0 .. 8,
         labels = [eta, f(eta)],
         axes = boxed
      );
  plot( [ seq((fpeta||j)(eta), j = 1..6)],
         eta = 0 .. 5,
         labels = [eta, fprime(eta)],
         axes = boxed
      );
plot( [ seq((phipeta||j)(eta), j = 1..6)],
         eta = 0 .. 5,
         labels = [eta, phiprime(eta)],
         axes = boxed
      );

plot( [ seq((phieta||j)(eta), j = 1..6)],
         eta = 0 .. 5,
         labels = [eta, phi(eta)],
         axes = boxed
      );

Warning, expecting only range variable eta in expression Y_sol1(eta) to be plotted but found name Y_sol1

 

 

Warning, expecting only range variable eta in expression YP_sol1(eta) to be plotted but found name YP_sol1

 

 

Warning, expecting only range variable eta in expression feta1(eta) to be plotted but found name feta1

 

 

Warning, expecting only range variable eta in expression fpeta1(eta) to be plotted but found name fpeta1

 

 

Warning, expecting only range variable eta in expression phipeta1(eta) to be plotted but found name phipeta1

 

 

Warning, expecting only range variable eta in expression phieta1(eta) to be plotted but found name phieta1