Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

To Maple support:

I was investigating this pde from a different forum.

I noticed that when using an expanded version of the pde, Maple hangs. Without expanding the PDE, Maple gives an answer in 2 seconds. 

Why does expanding the PDE makes a difference? I do not have an earlier version of Maple on my new PC to check if this is a new issue or not.
 

interface(version);

`Standard Worksheet Interface, Maple 2022.0, Windows 10, March 8 2022 Build ID 1599809`

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 1230 and is the same as the version installed in this computer, created 2022, April 21, 9:8 hours Pacific Time.`

restart;
pde1:=VectorCalculus:-Laplacian(u(r,theta),'polar'[r,theta]);
pde1_expanded:=expand(pde1);
bc  := u(1,theta)=sin(theta)^4,u(3,theta)=1;
pdsolve([pde1=0,bc],u(r,theta))
 

(diff(u(r, theta), r)+r*(diff(diff(u(r, theta), r), r))+(diff(diff(u(r, theta), theta), theta))/r)/r

(diff(u(r, theta), r))/r+diff(diff(u(r, theta), r), r)+(diff(diff(u(r, theta), theta), theta))/r^2

u(1, theta) = sin(theta)^4, u(3, theta) = 1

u(r, theta) = (1/52480)*((328*r^6-26568*r^2)*ln(3)*cos(2*theta)+(-r^8+6561)*ln(3)*cos(4*theta)+19680*(ln(3)+(5/3)*ln(r))*r^4)/(ln(3)*r^4)

pdsolve([pde1_expanded=0,bc],u(r,theta)); #HANGS, Waited more than 40 minutes.

 


 

Download hangs_pde.mw

Let 
                              "a"

 and 
                              "b"

 be real numbers and 

       
"A = Matrix(3, 3, [[a, a - 1, -b], [a - 1, a, -b], [b, b, 2*a - 

   1]])"


,  
                              "B="


 "Matrix(5, 5, [[0, a, 3, 0, a], [3, 0, 0, b, 0], [0, 1, b, 0, 

    1], [b, 0, 0, 1, 0], [0, a, 1, 0, b]])"


(a) Show that if 
                            "0 <= a"


                            "a <= 1"

 and 
                      "b^2 = 2*a*(1 - a)"

, then A is an orthogonal matrix with determinant equal to one. 
(b) For what values of a and b is the matrix B singular? Determine the inverse of B (for those values of a and b for which B is invertible).
 

Find all rational function solutions to the Kadomtsev-Petviashvili equation 
                         
(&PartialD;)/(&PartialD;x);

diff(u, t) + 6*u*diff(u, x) + diff(u, x, x, x) - diff(u, y, y) = 0;

by u = 2 
diff(ln, x, x)*f;

=(2 (((&PartialD;)^2)/(&PartialD;x^2) f) f-2 ((&PartialD;)/(&PartialD;x) f)^2)/(f^2);
 with 
f;
  =
(a[1 ]x+a[2] y+a[3] t+a[4])^2+(a[5] x+ a[6] y+a[7] t+a[8])^2+a[9], ;
where 
a[i], i=1..9, ;
are real constants.
 

I am having a problem in evaluating a fluid flow problem at a boundary that includes infinity using a semi-analytic method Adomian Decomposition Method.

Find the attached code. Thanks in advance.

Just wondering how I do a loop over the contents of the table. I need both the index variable and the stored value associated with it.

This one apparently doesn't work.

F_vRk := table()

table( [ ] )

(1)

F_vRk["a"] := 1

1

(2)

F_vRk["b"] := 2

2

(3)

"for ind, val in F_vRk  do print(ind,val)  end do"

1, F_vRk

(4)

NULL

Download table_and_loop_1.mw

Hii!!

I need your help in state space system.... kindly guide me how to solve State Space system in maplesoft.....i excecute the command but i didn't find the answer....can you plz help me?I have been trying for two weeks now but it is not working.Thank you!!

Hello, we need to do an integral based on two numerically approximated functions given by a second degree coupled differential equation.

The integral:
Where xs and xu are functions of the variable t.

We wish to do a half-sums to approximate the integral, that is:

Define an N = number of points, such that we get a delta_x = L/N.

We then wish to get two vectors x_s = [x_s(0),x_s(1),...,x_s(N)] and x_u = [x_u(0),x_u(1),...,x_u(N)] by running some kind of for loop over our approximation procedure.

Is there any function for this already inside of Maple? We need to approximate the integral using Riemann sums and no other method.

When trying to print a Maple Document as a PDF on Mac, I am now getting an "Error while printing" message. A 1-page pdf is created that's empty.   This is a new error, one that I've never seen when printing as PDF for any other Maple Document I've created.

I'm running Mac OS 12.3.1, the latest version of Mac OS and using Maple 2020.

I am trying to develop a recursive code to the above equations.  Note, U,X&Y are multivariate functions (in this case bivariate functions of (x,y)) i.e. U=U(x,y), X=X(x,y) etc.

Dear all

I have a nonlinear system of algebraic equations, I would like to solve it  without using fsolve and Newton's method. 
Maybe one can use, fixed point method or Broyden's method

Fixed_Broyden_method.mw

Can we solve the system of 4 equations, using the proposed methods

thank you 

Hello;

Hope you are fine. Can i apply numerical scheme on maple for the following problem. This in integro-differential equation i think. Waiting for kind response.

Thanks

 


could you please help me ,the maple code for this given series.

restart

U[0](x) := x;

x

(1)

"U[k+1](x):=solve((k+1)*U[k+1](x)+(x*(-1)^((k-1)/(2)))/(k!)-x^(2)*((e)^(x))/(10){6/(k!)-sum((2^(k[1]))/(k[1!])(5*delta[k[]-k[1]](x)+(2^(k-k[1])*(-1)^((k-k[1])/(2)))/((k-k[1])!)+(2^(k-k[1]+1)*(-1)^((k-k[1]-1)/(2)))/((k-k[1])!)),k[1]=0..k)}-(cos^(2)(x)+sin^()(x))*((&PartialD;)^2)/((&PartialD;)^( )x^2) [U[k](x)]-(e)^(x)[sum(1/(k[1]!){1/(k-k[1])(sum(sum(1/(k[3]!)*U[k[2]-k[3]](x)*U[k-k[1]-k[2]-1](x)},k[1]=0..k),k[2]=0..k-k[1]-1),k[3]=0..k[2])),U[k+1](x)];  od;"

Error, unable to match delimiters

"U[k+1](x):=solve((k+1)*U[k+1](x)+(x*(-1)^((k-1)/2))/(k!)-x^2*((e)^x)/10{6/(k!)-sum((2^(k[1]))/(k[1!])(5*delta[k-k[1]](x)+(2^(k-k[1])*(-1)^((k-k[1])/2))/((k-k[1])!)+(2^(k-k[1]+1)*(-1)^((k-k[1]-1)/2))/((k-k[1])!)),k[1]=0..k)}-(cos^2(x)+sin(x))*((&PartialD;)^2)/(&PartialD;x^2) [U[k](x)]-(e)^x[sum(1/(k[1]!){1/(k-k[1])(sum(sum(1/(k[3]!)*U[k[2]-k[3]](x)*U[k-k[1]-k[2]-1](x)},k[1]=0..k),k[2]=0..k-k[1]-1),k[3]=0..k[2])),U[k+1](x)];  od;"

 

``


 

Download Chapter_6-Example-6.5.4.mw

Hi! Do you know maybe how to solve equation with Laplace operator in Maple like BZ equation?

restart;
a := 0.75;
rho := u(t) + v(t) + w(t);
ode := diff(u(t), t) = 10*Delta(u(t)) + u*(-a*v - rho + 1);
ode1 := diff(v(t), t) = 10*Delta(v(t)) + v*(-a*w - rho + 1);
ode2 := diff(w(t), t) = 10*Delta(w(t)) + w*(-a*u - rho + 1);

 

Edit: Sorry I guess it should be function of three variables so u,v,w depends on (x,y,z) not strictle from time

I am wondering how to animate something like this from BZ equation:

Let the curves y=2x^3-x^2-5x   and 𝑦=−𝑥^2+3𝑥be given.
a) Plot both curves on the same 𝑥y −axis.
b) Shade the region (between the curves).
c) Calculate the area of the region enclosed by the curves given.

Hi everyone! I'd really appreciate if I could get pointed in the right direction as I am a brand new maple user.

So im trying to solve this constrainted optimization problem (See picture) using Maple symbollically. I believe I should have a closed form solution given I can substitute the one constraint into the objective function. Specifically closed form solutions for the three phi variables.

Can someone point me in the right direction as to how I should go about this? I've already taken first order conditions and tried to using the solve() function to no avail, realizing my sytem of equations weren't linear );. 

First 33 34 35 36 37 38 39 Last Page 35 of 2097