Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Can anyone please assist with these errors?

#Clear memory.
restart;

#Initialise variables and arrays.
h:=Pi/10;
n:=10;
k:=0.1;  
m:=6;
t:=Array(0..m):
x:=Array(0..n):
u:=Array(0..n,0..m):

#Initialise the x array and the initial u(x,0) boundary.
for i from 0 to n do
    x[i]:=i*h;
    u[i,0]:=0;
end do:

#Initialise the t array and the u(x,t) side boundaries.
for j from 0 to m do
   t[j]:=j*k;
   u[0,j]:=0;
   u[n,j]:=3*sin*t[j];
end do:

#Use the 2D CTCS explicit wave method.
for i from 1 to n-1 do
   u[i,1]:=(u[i-1,0]/2*Pi^2+sin(x)/10*Pi^2+u[i+1,0]/2*Pi^2)-1;
end do:
 
for j from 1 to m-1 do
  for i from 1 to n-1 do
   u[i,j+1]:=(((u[i-1,j]+u[i+1,j]-u[i,j-1])/Pi^2)-2);
  end do;
end do:

#Display the u(x,t)values.
printf("2D CTCS Explicit Wave Method\n");
printf("----------------------------\n");
printf("x\t\t t\t\t u\n");
for i from 0 to n do
   printf("% f\t% f\t% f\n",x[i],t[m],u[i,m]);
end do;

(1/10)*Pi

 

10

 

.1

 

6

 

Error, invalid input: sin expects its 1st argument, x, to be of type algebraic, but received Array(0..10, [0,1/10*Pi,1/5*Pi,3/10*Pi,2/5*Pi,1/2*Pi,3/5*Pi,7/10*Pi,4/5*Pi,9/10*Pi,Pi])

 

2D CTCS Explicit Wave Method
----------------------------
x                 t                 u
 0.000000         0.600000         0.000000
 0.314159         0.600000        -2.015556
 0.628319         0.600000        -2.198509
 0.942478         0.600000        -2.215724
 1.256637         0.600000        -2.218015
 1.570796         0.600000        

 

Error, (in fprintf) number expected for floating point format

 

NULL


Download Asst_4_Question_2c.mw

Both, Maple and MapleSim in Versions from 2021 to 2023 do not start. A logfile is saved which starts as

# A fatal error has been detected by the Java Runtime Environment:

#

#  EXCEPTION_ACCESS_VIOLATION (0xc0000005) at pc=0x0000000000000000, pid=17256, tid=10520

#

# JRE version: OpenJDK Runtime Environment Temurin-19.0.1+10 (19.0.1+10) (build 19.0.1+10)

# Java VM: OpenJDK 64-Bit Server VM Temurin-19.0.1+10 (19.0.1+10, mixed mode, tiered, compressed oops, compressed class ptrs, g1 gc, windows-amd64)

# Problematic frame:

# C  0x0000000000000000

#

This only occurs when the notebook (on which the error occurs) is connected to the domain network. If the notebook is not connected, Maple and MapleSim can be started and will continue to run when the network is reconnected.

Has anybody seen this before or any ideas what this could be?

Edit: This all happened under Windows 10 on a notebook that has been running Maple for 3 years and only recently displayed this error.

Edit2: Connection via an USB doc which also connects to an external monitor.

Hello.

I would like to find the numerical solution of the Ginzburg-Landau equation for the1D, written in the form 

(diff(z(x), x, x))*kappa^2-(1-d)*kappa^2*(diff(z(x), x))/x+z(x)-z(x)^3 = 0

For d=1 and with the boundary conditions z(-infinity)=0 and z(infinity)=1 the solution is well-known, this is tanh(x) function.

My goal is to undertand the evolution of the numerical solution for the interval of the dimensionality 0<d<1. However, I faced with the problem that even for d=1 (the classical case) I get the error "Error, (in dsolve/numeric/bvp) Newton iteration is not converging".

Here is my simple code

restart;

with(plots):

R0 := (diff(z(x), x, x))*kappa^2-(1-d)*kappa^2*(diff(z(x), x))/x+z(x)-z(x)^3;

kappa := 2; d := 1;

p := dsolve({R0, z(-15) = 0, z(15) = 1}, type = numeric, range = -15 .. 15, maxmesh = 50000):

odeplot(p);

Thank you in advance.

Since I have updated from Maple 2021 to 2023 something is wrong with my stylesheets.

I rephrase my qeustion here:

I like to change the text-style of 2D-Input to courant 12, bold, color x, no italic.

Now it is courant 12, black, italic when opening a new worksheet.

I was going to format, style, 2D-input, modify and then the font has changed when opening a new woksheet but the color is going back to black and its also italic courant when pressing enter in that worksheet.

I When opening a new worksheet, i don't want ot looad an new style set by hand every time.

 

Thanks for support

 

 

How will I use maple 2016 to solve ODEs and showing the steps involved because this will increase my understanding in it. 

The worksheet below includes a sample use of the Mobius transformation which produces a hyperbolic reflection of a point in the Poincare disk, followed by my attempt to produce the same result from first principles in a procedure.

Do I misunderstand the Mobius transformation and its use and/or is my procedure incorrect? 

InversePoint.mw

Is there a way to directly code a complex conjugate such as z with overbar without using the verbose conjugate(z)?

Soit h : [1 ; + ∞[ℝ

1.a°) Complète le tableau ci-dessus et trace la courbe ( Ch) de h dans un repére orthonormé ( O,I,j )

 𝑥 1 2 3 4 5

h(𝑥 )

b°) Montre que h est une bijection

c°) Détermine la bijection réciproque h-¹ de h

d°) Calcule hoh-¹ et h-¹oh ( 𝑥 )

2.) Trace la courbe ( Ch-¹) dans le même repére que ( Ch)

I have solved a simple expression with different built-in methods.... I received differents answers ...What is the exact way to get the same answers by all methods

help_roots.mw.

Hi,

I am looking for a simple code to generate random magic squares (with the sum of integers) of dimensions 4x4 or 5x5.

Thank you

I have a very simple question:

How do I find the index of an list element.

For example to find 96 in the following list which is the 25th element.

R := [1, 2, 4, 8, 16, 32, 64, 31, 62, 27, 54, 11, 22, 44, 88, 79, 61, 25, 50, 3, 6, 12, 24, 48, 96, 95, 93, 89, 81, 65, 33, 66, 35, 70, 43, 86, 75, 53, 9, 18, 36, 72, 47, 94, 91, 85, 73, 49];
max(R);
                               96

So how to find 25 ?

I have the following procedure

export exportCsvMatrix := proc(m:: Matrix, fileName: string)

  ExportMatrix(fileName, m, target=csv):

end:

However, I would like to be able to pass a list of matrices as well.

How can I add a parameter type that is Matrix or list of Matrix?

Good day to all the members who read this question.
I ask for your help to find the right way to solve this two differential equations (is attached the file). Used the common code "dsolve" but is not working on this problem. Doesn't matter if the solution is numeric or algebraic.
I will thaks a lot your kind help.
Best regards to all of you

DIFFERENTIAL_EQUATION.mw

restart

with(PDEtools)

First Part

The differential equation to solve:

(Delta*LinearAlgebra:-Transpose(D[1/2])*D__0-I*`&mu;__e`*D__0/(lambda+I*`&mu;__e`*r)-2*(`&mu;__e`^2*r^2+lambda^2))*R(r)

Definitions

D__n = `&PartialD;__r`+I*omega(a^2+r^2)/Delta+I*a*m/Delta+2*n*(r-M)/Delta

LinearAlgebra:-Transpose(D[n]) = `&PartialD;__r`-I*omega(a^2+r^2)/Delta-I*a*m/Delta+2*n*(r-M)/Delta

a := 1; M := 1; omega := 1; m := 1; `&mu;__e` := 1; lambda := 1

Delta := -2*M*r+a^2+r^2

I divided the differential equation in 3 parts (A, B, C).

A := (-I*omega(a^2+r^2)-I*a*m+r-M)*(diff(R(r), r)+I*omega(a^2+r^2)*R(r)/Delta+I*a*m*R(r)/Delta)+Delta*(diff(R(r), r, r))+I*omega(a^2+r^2)*(diff(R(r), r))+I*a*m*(diff(R(r), r))

B := -I*`&mu;__e`*(diff(R(r), r)+I*omega(a^2+r^2)*R(r)/Delta+I*a*m*R(r)/Delta)/(lambda+I*`&mu;__e`*r)

C := -(2*(`&mu;__e`^2*r^2+lambda^2))*R(r)

DE := A+B+C

E := dsolve(DE)

R(r) = DESol({diff(diff(_Y(r), r), r)-(1-r+I/(1+I*r))*(diff(_Y(r), r))/(r^2-2*r+1)-(-(2*I)*((-1-2*I)+r)/(r^2-2*r+1)-2/((1+I*r)*(r^2-2*r+1))+2*r^2+2)*_Y(r)/(r^2-2*r+1)}, {_Y(r)})

(1)

dsolve({DE, DE(0) = 1}, numeric, range = 0 .. 20)

Error, (in dsolve/numeric/type_check) insufficient initial/boundary value information for procedure defined problem

 

Second Part

The differential equation to solve:

[`#msub(mi("L",fontweight = "bold"),mfrac(mn("1",fontweight = "bold"),mn("2",fontweight = "bold"),linethickness = "1"))`*LinearAlgebra:-Transpose(L[1/2])+a*`&mu;__e`*sin(theta)*LinearAlgebra:-Transpose(L[1/2])/(lambda+a*`&mu;__e`*cos(theta))+2*(lambda^2+a^2*`&mu;__e`*cos(theta)^2)]*S(theta) = 0

Definitions:

L__n = a*omega*`sin&theta;`+m*`cosec&theta;`+n*`cot&theta;`+`&PartialD;__&theta;`

LinearAlgebra:-Transpose(L[n]) = -a*omega*`sin&theta;`-m*`cosec&theta;`+n*`cot&theta;`+`&PartialD;__&theta;`

Also I divided the differential equation in 3 parts (A, B, C).

F := (omega*a*sin(theta)+m/sin(theta)+1/(2*tan(theta)))(diff(S(theta), theta)-omega*a*sin(theta)*S(theta)-m*S(theta)/sin(theta)+S(theta)/(2*tan(theta)))+diff(S(theta), theta, theta)-omega*a*sin(theta)*(diff(S(theta), theta))-m*(diff(S(theta), theta))/sin(theta)+(diff(S(theta), theta))/(2*tan(theta))

G := a*`&mu;__e`*sin(theta)*(diff(S(theta), theta)-omega*a*sin(theta)*S(theta)-m*S(theta)/sin(theta)+S(theta)/(2*tan(theta)))/(lambda+a*`&mu;__e`*cos(theta))

H := (2*(lambda^2+a^2*`&mu;__e`*cos(theta)^2))*S(theta)

DF := F+G+H

dsolve(DF)

dsolve({DF, DF(0) = 1}, numeric, range = 0 .. 20)

Error, (in dsolve/numeric/type_check) insufficient initial/boundary value information for procedure defined problem

 

NULL

Download DIFFERENTIAL_EQUATION.mw

I have a very difficult time finding things in the documentation.

I would like to see documentation for `<,>` and `<|>`, as in `<,>`(1,2,3,4).

If I write ?<|> in Maple I get

If I do ?<,> I get

However, there is no mention of the usage of the backticks.

Finally, if I do ?`<,>` I get

No mention of the use of `<,>`. I only know about this from a previous reply to one of my questions and now I have forgotten the details so I am trying to look it up but it's hard to find.

Dear Power Users, I tried to make a model in Maple but after several hours of struggling I need some help. Can someone tell me what I did wrong or explain to me what is the meaning of the maplemint comment on procedure (optimalisatie)? Thank you in advance.

DiLeoCasper.mw 

First 160 161 162 163 164 165 166 Last Page 162 of 2218