what does _Z and Z_^2 mean?...

I'm solving this equation:

eq3 := 2*lambda*(q[1]+p[1]*exp(t))*(q[2]+p[2])^2*p[1]*exp(t) = s[1]

eq4 := 2*lambda*(q[1]+p[1])^2*(q[2]+p[2]*exp(u))*p[2]*exp(u) = s[2]

solve({eq3, eq4}, {t, u})

This is what I get:

{t = ln(RootOf(-s[1]+(2*lambda*p[2]*q[1]+2*lambda*q[1]*q[2])*_Z+2*lambda*_Z^2)/((q[2]+p[2])*p[1])), u = ln(RootOf(-s[2]+(2*lambda*p[1]*q[2]+2*lambda*q[1]*q[2])*_Z+2*lambda*_Z^2)/((q[1]+p[1])*p[2]))}

I'm completely new at maple.

I can't find what _Z  and _Z^2 mean. Can anybody help? and how I can solve this _z and _Z^2

Hello,

i would like to continue to work on a code which was done With CodeBuilder some years ago.

i often like to print my code to be able to read without any computer.

Problem : the print function on a worksheet with several coderegion doesn t work well.

1) Do you have some ideas to circumvet this issue ,

2) i m thinking about using a text éditor which can give a nice printing with maple. Emacs is far to much complex for me. Consequently, i havé seen that some of you have testéd Sublime text. I think that it may be a good alternative. Do you have ideas about simple text editor that i can use to edit maple code ?

thank you for your help

Empty output from pdsolve...

hi, I wonder why this output is empty because I tried to get solution below partial equation:

restart;
with(PDEtools);
pe := diff(u(t, x, y), t) = u(t, x, y) - diff(u(t, x, y), x \$ 2) - diff(u(t, x, y), y \$ 2) + (1 + u(t, x, y)*I)*abs(u(t, x, y))*u(t,x,y)^2;
inc := u(0, x, y) = cos(Pi*x/50) - cos(Pi*x/150)*I;
sys := [pe, inc];

pdsolve(sys);

end after this pdsolve(sys); Maple doesn't show anything

Matlab code generation ... ???...

When I want to generate matlab code of the following symbolic expression I got very strange result:

```with(CodeGeneration); Matlab(Sum(a(i), i = 1 .. N)); Warning, the function names {Sum, a} are not recognized in the target language Warning, precedence for Range unspecified Warning, cannot translate range cg = Sum(a(i), i == (1..N));```

Problem:`Sum(a(i), i == (1..N));`

Right matlab code should look like:

`cg = sum(a);`

Where can I found compete list of supported matlab built-in function? Is there any possibility how to get proper matlab code?

Is there any way how to convert more complex symbilic expressions to the Matlab???

Is "Delete As Spam" broken?...

The last several times that I've tried to use the Delete As Spam feature I've gotten an error message Error Generating Page, and the spam was not deleted. Anyone know what's happening?

How do I solve the heat equation in spherical coor...

I'm trying to understand what I need to do to get Maple to solve the heat equation in spherical coordinates. I'm considering a test problem: a spherical shell with inner radius 1 and outer radius 2. Boundary conditions require the temperature of the inner and outer surfaces are to always equal zero. The initial condition requires the temperature inside the shell to be -sin(pi*r).

The attached worksheet provides my attempt to solve this problem. Executing this worksheet will quickly give a general solution but this does not include the initial condition. If this is included then Maple fails to give a solution even after a 20 min wait.

Could someone offer me guidance on how to proceed?

shell.mw

about table and MultiSet for the construction of d...

Hello,

I have code which contains a data structure. This data structure was defined with Types and a Table.

I would like to know which is the best manner to do to build a data structure.

Is it better to use a Table or a MultiSet ?

May you precise me the pros and the cons of the 2 solutions?

Or indicate me if you have other ways to do.

Thanks a lot for your help.

Is it possible to Add/sum by mapping multiple valu...

I want to sum a series by getting multiple values from a list.  Based on the list L I was trying

add(A, [m2, m3] in L)   This seem to have dissapeared from the document but the error message remains.

 (1)

 (2)

I found this code by Kitonum 16740  https://www.mapleprimes.com/questions/138955-How-Do-I-Create-A-Variable-Number-Of-Nested-Loops-

 (3)

 (4)

 (5)

Instead of using I would like to map m2 and m3 to the values in list L and sum them. This should give a method for having a variable number of loops

I tried variants off without success

for the next higher order would need to map m2,m3, m4 to the generated list.

 (6)

 (7)

 (8)

how to create the function that finds the the sum ...

sum2N := proc (N::integer) local summation, i, k; summation := 0; k := 0;

for i from 0 to N do if i = 1 then k := k+1

end if; summation := (i+k)^2 end do

end proc;
sum2N(10);

the result shows 121 but its not correct i think the right answer show be 2985

Transformer les couples (x,y) de résultat d'une éq...

Bonjour,

Comment transformer cette liste de résultat:

a := {x = -12, y = -23}, {x = -34, y = 87}, {x = 18, y = 98}, {x = 27, y = 93}, {x = 45, y = 68}

En une liste simple telle que L: = [-12,-23,-34,87,18,98,27,93,45,68 ]

Merci

HELLO,

How to transform this result list :

a := {x = -12, y = -23}, {x = -34, y = 87}, {x = 18, y = 98}, {x = 27, y = 93}, {x = 45, y = 68}

In a simple list such that :  L: = [-12,-23,-34,87,18,98,27,93,45,68 ]

Thank You

Sums as hypergeometric. Why?...

I have two summations that Maple converts to hypergeometric outputs. Is there a way to stop that?

I dont know how to get Maple to expand the Hyper geometric output for a given walue of N

if the value of N is defined first the output is as I would expect.

The real problem, the sumations are nested and when I get to four nested sums Maple sits there for hours trying to evaluate the general case.

I have only included the first 2 summations.

Any insights on this would be appreciated.

 (1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

Plotting the first and second derivative of the nu...

Hello, I am solving an ODE numerically:

restart: deq1 := diff(y(t), t, t)+y(t)*abs(y(t)) = 0; ic1 := y(0) = 1.0, (D(y))(0) = 0.;
dsol1 := dsolve({deq1, ic1}, numeric, range = 0 .. 10);
odeplot(dsol1);

I need to get the plots of the first and second derivatives of the solution. Ok, I can get the first derivating by transforming the equation to system with two unknown functions and reducing the order, but I also need the second and possibly higher derivatives.

highly anisotropic problems with phase field model...

I am a student who is learning to use Maple to perform numerical analysis of phase field. There has been no solution to the problem of high anisotropy. The reference code is as follows

W :=piecewise(beta>=Pi/2+thetaM and beta<=Pi-thetaM, 1+eps*cos(4*beta),
beta>=Pi/2-thetaM and beta<=Pi/2+thetaM, (1+esp*cos(4*thetaM))/cos(thetaM)*cos(4*beta));

However, there is no problem at all when simply performing the following settings.

W := 1+ eps*cos(4*beta);

Please tell me how to write such a conditional maple statement. I have referenced many examples on the official website, but none of them succeeded. This includes writing the code in the following form.

W :=betta->piecewise(beta>=Pi/2+thetaM and beta<=Pi-thetaM, 1+eps*cos(4*beta),
beta>=Pi/2-thetaM and beta<=Pi/2+thetaM, (1+esp*cos(4*thetaM))/cos(thetaM)*cos(4*beta));

Remarks, beta is an angle, set as follows.

subs(beta=t_atan2(diff(pho(x,y),y), diff(pho(x,y),x))

Straight line with fixed direction...

We give a line (D) and a point A located at a distance AH=h from D. A constant angle of magnitude alpha pivots to its apex A and we call B and C the points where its sides cut the line D. Let O be the center of the circle circumscribed to the triangle ABC.
Demonstrate that the B and C tangents to the O circle keep a fixed direction.
Here is my code which don't work for slopes are not equal.

restart; with(plots): with(geometry):unprotect(D):
_EnvHorizontalName := 'x':_EnvVerticalName := 'y':
line(D, y = (1/2)*x-1); point(A, 5, 5); PerpendicularLine(lp, A, D); h := distance(A, D); intersection(H, D, lp);
alpha := (1/16)*Pi;
rotation(lp1, lp, (1/6)*Pi, 'clockwise', A); rotation(lp2, lp1, (1/6)*Pi-alpha, 'clockwise', A); FindAngle(lp1, lp2); evalf(%);
intersection(B, D, lp1); intersection(C, D, lp2);
triangle(T, [A, B, C]);
circumcircle(Elc, T, 'centername' = OO);
TangentLine(tgB, B, Elc); TangentLine(tgC, C, Elc);
evalf(slope(tgB)); evalf(slope(tgC));
dr := draw([D(color = blue), lp(color = red), Elc(color = green), A, B, C, T(color = black), H, tgB, tgC], printtext = true);

display([dr], axes = none, scaling = constrained);
Fig := proc (k) local dr, Elc, B, C, lp1, lp2; global D, A, lp, H, alpha; geometry:-rotation(lp1, lp, (1/6)*Pi+k, 'clockwise', A); geometry:-rotation(lp2, lp1, (1/6)*Pi-alpha+k, 'clockwise', A); geometry:-intersection(B, D, lp1); geometry:-intersection(C, D, lp2); geometry:-triangle(T, [A, B, C]); geometry:-circumcircle(Elc, T, 'centername' = OO); geometry:-TangentLine(tgB, B, Elc); geometry:-TangentLine(tgC, C, Elc); dr := geometry:-draw([D(color = blue), lp(color = red), Elc(color = green), A, B, C, T(color = black), H, tgB, tgC], printtext = true); plots:-display([dr], axes = none, scaling = constrained) end proc;
iframes := 10;

display([seq(Fig((1/12)*Pi+i/(10*iframes)), i = 1 .. iframes)], insequence, scaling = constrained);
How to improve this code ? Thank you.

 First 159 160 161 162 163 164 165 Last Page 161 of 2173
﻿