Maple 2019 Questions and Posts

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

I'm using Maple2019.

Here is my code:

with(geometry);
point(o, 0, 0);
point(A, 0, 1);
point(d, 0, 2);
point(F, 0.8944271920, 1.4472135960);
line(lOD, [o, d]);
line(lAF, [A, F]);
alpha := FindAngle(lOD, lAF);
 =
                      alpha := 2.034443936


FindAngle is supposed to return the smaller angle between two lines.

Here it is greater than Pi/2

How can that be ?

Thank you for your advice

Hello, I was wondering about the following.

cos(Pi/2^n) can be converted to a radical for each integer n, because cos(Pi/2)=0 is known and successively

cos(x)=sqrt( (1+cos(2x))/2 )

so in principle radical expressions should exist for all integers n. But maple only converts cos(Pi/8) to a radical, but not cos(Pi/16) and higher powers. Is there a simple way without writing an own program to obtain these radical expressions using convert(*,radical)?

restart;
with(LinearAlgebra):
with(VariationalCalculus):
# Energie cinetique/potentiel
T := 0.5*m*(diff(r(t),t)^2+r(t)^2*diff(theta(t),t)^2)+(1/6)*m[0]*l^2*diff(theta(t),t)^2;
Ve := 0.5*k*(r(t)-r[o])^2;
Vg := -m*g*r(t)*cos(theta(t)) - 0.5*m[o]*g*l*cos(theta(t));
# Lagrangien
L := T - (Ve + Vg);
E := EulerLagrange(L,t,[theta(t),r(t)]);
EQ1:=simplify(E[1]);
EQ2:=simplify(E[2]);
dsolve({EQ1,EQ2,D(theta)(0)=0,D(r)(0)=0,theta(0)=0,r(0)=10},{r(t),theta(t)});
dsolve({EQ1,EQ2},{r(t),theta(t)});

Hello ! 

I'm currently trying to learn some advanced mechanics with the Lagrangian method. I'm ok with the physics part. I can solve small problems by hand but when I try it on maple, it won't work.

The equations I'm using have no mistakes, I'm following a book step by step. 

What I'm stuck with is to get a formula of theta and r that are only function of time (t) to then be able to plot it. And I need dsolve to work to get that.

 

 

 

 

As of last week using Maple 2019, and now with Maple 2020, I have begun to have issues trying to export rendered plots using plot3d. Even trying a simple example such as:

plot3d(x^2+y^2, x=0..1, y=0..1); 

will return the correct 3D plot, but right-clicking the plot closes the program immediately. Has anyone else experienced a similar issue? 

I am trying to graph x^(1/2) and (x^2/8) revolving around the y-axis using the ring method but I have not been able to do it. I only get the shell method. I have adjusted the functions to read x=y^2 and x=sqrt(8y). The bounds for this method is 0 to 2. No matter the commands I have tried I have not been able to get the washer method to show. I do not want to use the tutor program.

Can someone help me with the following question?

 

I would like to plot the function f(x)=a*x*(1-x)+x*ln(x)+(1-x)ln(1-x) in the interval [x1,x2], where x1 and x2 are points that satisfy the following: df/dx(x=x1)=df/dx(x=x2)=0 and f(x1)=f(x2), where a is some parameter that is greater than 0.

I guess I need to use here fsolve, but I am not sure how.

x:=Proc
Local a,
b,

C;
use XML Tools in 
end use;

end proc;

Hello,

I have a problem, please help me to solve it.

I am trying to extract coefficients of some polynomial expressions but I receive Error.

For example:

>with(Physics):

>Setup(mathematicalnotation = true)

>S := a*(x^2)+b*x+c*(1/(x-h)):

>Coefficients(S, x, 1)

>Error, (in Physics:-Coefficients) unable to compute coeff.

**** Note that if instead of the term (1/(x-h)) we have (1/x), then we get the true result. The error occurs only when we have some extra things in the denominator namely "-h".

Regards

How to learn Maple Programming effectively, whether Python will help?

My objective is to write a procedure that will read data from xml at different mentioned condition, read ICBO value at different condition of VGS i.e ICBO value at VGS=5V

BJT_ICBO := proc(parsedXML, VGS at condition )
 

 

end proc;

can some one help me to write procedure 

 

In maple, how to change strings style? Thanks!  

In Mathematica , we can choose Style to change strings style.

Style["A", Blue, Italic, 24]

 

source of the  problem : I make an User Interface.

the color of   "mass "  may consider be changed to  red . 

 

Hi, can someone explain me why CodeGeneration in C throws the error : Error, (in Print) improper op or subscript selector ? Interestingly, CodeGeneration in Python or Matlab works. Source code is attached. Thanks in advance!

membrane_energy.mw
 

 

Membrane Energy

 

restart; with(VectorCalculus); with(LinearAlgebra); with(CodeGeneration)

Define variables

 

Vectors for vertices of current position

v1 := Vector(3, symbol = v1_i) = Vector[column]([[v1_i[1]], [v1_i[2]], [v1_i[3]]], ["x", "y", "z"]) 

v2 := Vector(3, symbol = v2_i) = Vector[column]([[v2_i[1]], [v2_i[2]], [v2_i[3]]], ["x", "y", "z"])NULL

v3 := Vector(3, symbol = v3_i) = Vector[column]([[v3_i[1]], [v3_i[2]], [v3_i[3]]], ["x", "y", "z"])NULL

n := `&x`(v2-v1, v3-v1)

v4 := v1+n/norm(n, 2)^.5

Vector for vertices of next position

v1n := Vector(3, symbol = v1n_i) = Vector[column]([[v1n_i[1]], [v1n_i[2]], [v1n_i[3]]], ["x", "y", "z"])NULL

v2n := Vector(3, symbol = v2n_i) = Vector[column]([[v2n_i[1]], [v2n_i[2]], [v2n_i[3]]], ["x", "y", "z"])NULL

v3n := Vector(3, symbol = v3n_i) = Vector[column]([[v3n_i[1]], [v3n_i[2]], [v3n_i[3]]], ["x", "y", "z"])NULL

nn := `&x`(v2n-v1n, v3n-v1n)

v4n := v1n+nn/norm(nn, 2)^.5``

``

 

Define Transformation

 

V := LinearAlgebra:-Transpose(Matrix([v2-v1, v3-v1, v4-v1]))

Dimension(V) = 3, 3 

Vn := LinearAlgebra:-Transpose(Matrix([v2n-v1n, v3n-v1n, v4n-v1n]))

Dimension(Vn) = 3, 3NULL

Note we have Vn = T*V and if the current triangle is not degenerate, then T = Vn/V. As we can pre-compute 1/V we define a new matrix for it:

Vinv := Matrix(3, 3, symbol = Vinv_ij) =

Matrix(%id = 18446746713267339006)

(1.2.1)

T := MatrixMatrixMultiply(Vn, Vinv)

Dimension(T) = 3, 3``

``

``

Define Energy

 

E := Trace(MatrixMatrixMultiply(T, LinearAlgebra:-Transpose(T)))``

``

``

Gradient and Hessian

 

gradE := Gradient(E, [v1n[1], v1n[2], v1n[3], v2n[1], v2n[2], v2n[3], v3n[1], v3n[2], v3n[3]])

Dimension(gradE) = 9NULL

CodeGeneration[C](gradE, defaulttype = numeric, optimize = tryhard, functionprecision = double, precision = double, deducetypes = false, resultname = 'gradE')t1 = v2n_i[0] - v1n_i[0];
t2 = v2n_i[1] - v1n_i[1];
t3 = v2n_i[2] - v1n_i[2];
t4 = t1 * Vinv_ij[0][0] + t2 * Vinv_ij[1][0] + t3 * Vinv_ij[2][0];
t5 = t1 * Vinv_ij[0][1] + t2 * Vinv_ij[1][1] + t3 * Vinv_ij[2][1];
t6 = t1 * Vinv_ij[0][2] + t2 * Vinv_ij[1][2] + t3 * Vinv_ij[2][2];
t7 = v3n_i[0] - v1n_i[0];
t8 = v3n_i[1] - v1n_i[1];
t9 = v3n_i[2] - v1n_i[2];
t10 = t7 * Vinv_ij[0][0] + t8 * Vinv_ij[1][0] + t9 * Vinv_ij[2][0];
t11 = t7 * Vinv_ij[0][1] + t8 * Vinv_ij[1][1] + t9 * Vinv_ij[2][1];
t12 = t7 * Vinv_ij[0][2] + t8 * Vinv_ij[1][2] + t9 * Vinv_ij[2][2];
t13 = t2 * t9 - t3 * t8;
t14 = fabs(t13);
t15 = t1 * t9 - t3 * t7;
t16 = fabs(t15);
t17 = t1 * t8 - t2 * t7;
t18 = fabs(t17);
t19 = pow(t14, 0.2e1) + pow(t16, 0.2e1) + pow(t18, 0.2e1);
t20 = pow(t19, -0.5e1 / 0.4e1);
t19 = t19 * t20;
t21 = Vinv_ij[0][0] * t13;
t22 = Vinv_ij[1][0] * t15;
t23 = Vinv_ij[2][0] * t17;
t24 = (t23 + t21 - t22) * t19;
t25 = -v2n_i[2] + v3n_i[2];
t26 = fabs(t15) / t15;
t27 = -v2n_i[1] + v3n_i[1];
t28 = fabs(t17) / t17;
t21 = t23 + t21 - t22;
t22 = (t16 * t25 * t26 + t18 * t27 * t28) * t20;
t23 = 0.1e1 / 0.2e1;
t29 = Vinv_ij[0][1] * t13;
t30 = Vinv_ij[1][1] * t15;
t31 = Vinv_ij[2][1] * t17;
t32 = (t31 + t29 - t30) * t19;
t29 = t31 + t29 - t30;
t30 = Vinv_ij[0][2] * t13;
t15 = Vinv_ij[1][2] * t15;
t17 = Vinv_ij[2][2] * t17;
t31 = (t17 + t30 - t15) * t19;
t15 = t17 + t30 - t15;
t17 = t6 + t12;
t30 = t5 + t11;
t33 = t4 + t10;
t13 = fabs(t13) / t13;
t34 = -v2n_i[0] + v3n_i[0];
t35 = (t14 * t25 * t13 - t18 * t34 * t28) * t20;
t35 = t17 * Vinv_ij[1][2] - t24 * (t23 * t35 * t21 - t19 * (t25 * Vinv_ij[0][0] - t34 * Vinv_ij[2][0])) + t30 * Vinv_ij[1][1] - t31 * (t23 * t35 * t15 - t19 * (t25 * Vinv_ij[0][2] - t34 * Vinv_ij[2][2])) - t32 * (t23 * t35 * t29 - t19 * (t25 * Vinv_ij[0][1] - t34 * Vinv_ij[2][1])) + t33 * Vinv_ij[1][0];
t36 = (t14 * t27 * t13 + t16 * t34 * t26) * t20;
t34 = Vinv_ij[2][2] * t17 - t24 * (-t23 * t36 * t21 + t19 * (t27 * Vinv_ij[0][0] - t34 * Vinv_ij[1][0])) + t30 * Vinv_ij[2][1] - t31 * (-t23 * t36 * t15 + t19 * (t27 * Vinv_ij[0][2] - t34 * Vinv_ij[1][2])) - t32 * (-t23 * t36 * t29 + t19 * (t27 * Vinv_ij[0][1] - t34 * Vinv_ij[1][1])) + t33 * Vinv_ij[2][0];
t36 = (t16 * t9 * t26 + t18 * t8 * t28) * t20;
t36 = t24 * (-t23 * t36 * t21 - t19 * (-t8 * Vinv_ij[2][0] + t9 * Vinv_ij[1][0])) + t31 * (-t23 * t36 * t15 - t19 * (-t8 * Vinv_ij[2][2] + t9 * Vinv_ij[1][2])) + t32 * (-t23 * t36 * t29 - t19 * (-t8 * Vinv_ij[2][1] + t9 * Vinv_ij[1][1])) + t4 * Vinv_ij[0][0] + t5 * Vinv_ij[0][1] + t6 * Vinv_ij[0][2];
t37 = (t14 * t9 * t13 - t18 * t7 * t28) * t20;
t9 = t24 * (-t23 * t37 * t21 + t19 * (-t7 * Vinv_ij[2][0] + t9 * Vinv_ij[0][0])) + t31 * (-t15 * t23 * t37 + t19 * (-t7 * Vinv_ij[2][2] + t9 * Vinv_ij[0][2])) + t32 * (-t23 * t29 * t37 + t19 * (-t7 * Vinv_ij[2][1] + t9 * Vinv_ij[0][1])) + t4 * Vinv_ij[1][0] + t5 * Vinv_ij[1][1] + t6 * Vinv_ij[1][2];
t37 = (t14 * t8 * t13 + t16 * t7 * t26) * t20;
t4 = t24 * (t23 * t37 * t21 - t19 * (-t7 * Vinv_ij[1][0] + t8 * Vinv_ij[0][0])) + t31 * (t15 * t23 * t37 - t19 * (-t7 * Vinv_ij[1][2] + t8 * Vinv_ij[0][2])) + t32 * (t23 * t29 * t37 - t19 * (-t7 * Vinv_ij[1][1] + t8 * Vinv_ij[0][1])) + t4 * Vinv_ij[2][0] + t5 * Vinv_ij[2][1] + t6 * Vinv_ij[2][2];
t5 = (t16 * t3 * t26 + t18 * t2 * t28) * t20;
t5 = t10 * Vinv_ij[0][0] + t11 * Vinv_ij[0][1] + t12 * Vinv_ij[0][2] + t24 * (t21 * t23 * t5 + t19 * (-t2 * Vinv_ij[2][0] + t3 * Vinv_ij[1][0])) + t31 * (t15 * t23 * t5 + t19 * (-t2 * Vinv_ij[2][2] + t3 * Vinv_ij[1][2])) + t32 * (t23 * t29 * t5 + t19 * (-t2 * Vinv_ij[2][1] + t3 * Vinv_ij[1][1]));
t6 = (-t18 * t1 * t28 + t14 * t3 * t13) * t20;
t3 = t10 * Vinv_ij[1][0] + t11 * Vinv_ij[1][1] + t12 * Vinv_ij[1][2] + t24 * (t21 * t23 * t6 - t19 * (-t1 * Vinv_ij[2][0] + t3 * Vinv_ij[0][0])) + t31 * (t15 * t23 * t6 - t19 * (-t1 * Vinv_ij[2][2] + t3 * Vinv_ij[0][2])) + t32 * (t23 * t29 * t6 - t19 * (-t1 * Vinv_ij[2][1] + t3 * Vinv_ij[0][1]));
t6 = (t16 * t1 * t26 + t14 * t2 * t13) * t20;
t1 = t10 * Vinv_ij[2][0] + t11 * Vinv_ij[2][1] + t12 * Vinv_ij[2][2] + t24 * (-t21 * t23 * t6 + t19 * (-t1 * Vinv_ij[1][0] + t2 * Vinv_ij[0][0])) + t31 * (-t15 * t23 * t6 + t19 * (-t1 * Vinv_ij[1][2] + t2 * Vinv_ij[0][2])) + t32 * (-t23 * t29 * t6 + t19 * (-t1 * Vinv_ij[1][1] + t2 * Vinv_ij[0][1]));
t2 = 0.2e1;
gradE[0] = -t2 * (-t24 * (t23 * t22 * t21 + t19 * (t25 * Vinv_ij[1][0] - t27 * Vinv_ij[2][0])) + t33 * Vinv_ij[0][0] + t30 * Vinv_ij[0][1] + t17 * Vinv_ij[0][2] - t32 * (t23 * t22 * t29 + t19 * (t25 * Vinv_ij[1][1] - t27 * Vinv_ij[2][1])) - t31 * (t23 * t22 * t15 + t19 * (t25 * Vinv_ij[1][2] - t27 * Vinv_ij[2][2])));
gradE[1] = -t2 * t35;
gradE[2] = -t2 * t34;
gradE[3] = t2 * t36;
gradE[4] = t2 * t9;
gradE[5] = t2 * t4;
gradE[6] = t2 * t5;
gradE[7] = t2 * t3;
gradE[8] = t2 * t1;

``

``

``

``

``

Hessian

 

``

hessE := Hessian(E, [v1n[1], v1n[2], v1n[3], v2n[1], v2n[2], v2n[3], v3n[1], v3n[2], v3n[3]])

Dimension(hessE)

9, 9

(1.5.1)

 

CodeGeneration[C](hessE, optimize = tryhard, deducetypes = false, resultname = 'hessE')

Error, (in Print) improper op or subscript selector

 

``

````

``

``

``

``


 

Download membrane_energy.mw

 

 

 

Hello experts,

I am trying to use solve to find a solution to a system of two equations. The two equations involved are quite complicated, and so sometimes Maple has difficulty with this. In particular, the solve command tries evaluating but never comes up with a solution (I've waited at least an hour, before giving up). 

I recently learned about fsolve, which gives approximate numeric solutions (which would be fine for my purposes), but fsolve too struggles with a solution and simply returns my input to me. I tried plotting the system of equations using plots:-implictplot to see if a solution existed, and as expected it does. I was hoping to get some tips on trying to solve a difficult system like this, perhaps given the knowledge that a solution definitely exists. Unfortunately, I need a solution for many variations of the same system, so simply reading off the approximate solution isn't really an option.

In my attatched code, the system with W = 49 is the first one where Maple really begins to struggle, and I believe that solutions for W>49 are also difficult.

06042020_Predicting_w_AB_Ratio_Maple_Primes.mw

Thanks!

Hello again, it's my first time using maple, so I have more problems :(

I need solve three equations, but maple shows an error:

Error, (in solve) a constant is invalid as a variable, gamma

gamma.mw

I have used g instead of gamma and solve works, but I don't understand what happen.

 

The second problem is when I use g instead of gamma. In some tau values solve() doesn't show the solution. In each solve() always there are one unnecessary equation. Maybe could be that. But I don't know.

g.mw

 

Thank you in advance!!!

 

P.S. Sorry for my bad English

 Hello! I'm trying to solve the following: 

pde1 := (y+z)*(diff(u(x, y, z), x))+(z+x)*(diff(u(x, y, z), y))+(x+y)*(diff(u(x, y, z), z)) = 0;
{pdsolve(pde1,u(x,y,z))}

Unfortunately, after calling pdsolve , I get an empty result set . Can you help me figure out what's going on?  Does it  really have no solution? 

 

 

 

 

 

First 20 21 22 23 24 25 26 Last Page 22 of 47