Maple 2019 Questions and Posts

Error, '_Inert_DELAYLESSTHAN' is not a valid inert...

Asked by: amanuit11 20

April 11 2020

1 3

x:=Proc
Local a,
b,

C;
use XML Tools in
end use;

end proc;

Error in coefficient...

Asked by: John2020 220

April 11 2020

2 6

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 especia...

Asked by: amanuit11 20

April 10 2020

0 0

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

Procedure: How to pass conditional value in proced...

Asked by: amanuit11 20

April 10 2020

1 0

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

Which function is good for displaying strings f...

Asked by: lcz 1029

April 09 2020

0 5

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 .

CodeGeneration in C throws Error: Error, (in Print...

Asked by: sirkris 5

April 08 2020

1 0

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

Define variables

Vectors for vertices of current position

=

Vector for vertices of next position

=

Define Transformation

=

Note we have and if the current triangle is not degenerate, then . As we can pre-compute we define a new matrix for it:

=

Matrix(%id = 18446746713267339006)

(1.2.1)

=

Define Energy

Gradient and Hessian

=

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

(1.5.1)

Error, (in Print) improper op or subscript selector

Download membrane_energy.mw

Looking for tips to solve system of two (complicat...

Asked by: ajfriedlan 55

April 07 2020

1 5

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!

Two problems with solve()...

Asked by: John Blacksad 10

April 06 2020

0 2

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

Empty result set on `pdsolve` in Maple2019...

Asked by: lcz 1029

April 06 2020

0 2

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?

Get specific previous result...

Asked by: John Blacksad 10

April 05 2020

2 5

Hello, I have this:

maple.mw

I need save the three solutions in three different variables.

Thanks a lot!!!

Multiple variable integration using ArrayInterpola...

Asked by: Ingvars 10

April 03 2020

2 2

I am using Maple to load well-behaved function data (512 points) and using ArrayInterpolation integrate it:

P:= r -> CurveFitting:-ArrayInterpolation(xydata, Array(1 .. 1, 1 .. 1, [r]), method = spline)[1]; .

While the simple integral int(P, 0 .. 100, numeric, digits = 4) outputs a number, if I add another variable like int(P*a, [r = 0 .. 1, a = 0 .. 1], numeric) it just outputs (obviously the integral is more complex but this just shows the principle).

int(P(r)*a, [r = 0 .. 1, a = 0 .. 1], numeric) gives and error 'invalid input: coordinates of xvalues must be of type numeric'. Does anyone know the syntax to make this integral evaluate numerically?

I really wouldn't want to use spline() or interpolation() because as far as I have tried they are really slow and incompatable respectively. Actual document below.

Much thanks in advance.

Download Au_Scattering.mw

Fsolve returns expression instead of numbers...

Asked by: zjervoj97 40

April 02 2020

1 8

Hi, when trying to solve 2 equations with 2 variables the fsolve function reutrns just expression instead of numerical solution. Can someone please help with this? Thanks a lot

NRTL.mw

>

restart

>

T:=325:

>	a12:=2305.28444347652 - 9.14490843016421T + 0.00680052257590234T^2:

>	a21:=-6665.24838284836 + 46.0897018087247T - 0.0694991633494123T^2:

>	alfa:=0.3:

>

x2:=1-x1:

>	tau12:=a12/T:

>	tau21:=a21/T:

>	G12:=exp(-alfa*tau12):

>	G21:=exp(-alfa*tau21):

>	lng1:=x2^2(tau21(G21/(x1+x2G21))^2+tau12(G12/((x2+x1*G12)^2))):

>	lng2:=x1^2(tau12(G12/(x2+x1G12))^2+tau21(G21/((x1+x2*G21)^2))):

>	lnga1:=subs(x1=xa1,lng1):

>	lngb1:=subs(x1=xb1,lng1):

>	lnga2:=subs(x1=xa1,lng2):

>	lngb2:=subs(x1=xb1,lng2):

>	r1:=lnga1+ln(xa1)=lngb1+ln(xb1)

(1-xa1)^2*(.4966874722/(.4073000470+.5926999530*xa1)^2+.1510891213/(-0.464212743e-1*xa1+1)^2)+ln(xa1) = (1-xb1)^2*(.4966874722/(.4073000470+.5926999530*xb1)^2+.1510891213/(-0.464212743e-1*xb1+1)^2)+ln(xb1)

(1)

>	r2:=lnga2+ln(1-xa1)=lngb2+ln(1-xb1)

xa1^2*(.1440753718/(-0.464212743e-1*xa1+1)^2+1.219463331/(.4073000470+.5926999530*xa1)^2)+ln(1-xa1) = xb1^2*(.1440753718/(-0.464212743e-1*xb1+1)^2+1.219463331/(.4073000470+.5926999530*xb1)^2)+ln(1-xb1)

(2)

>	fsolve({r1,r2})

$fsolve({xa1^2*(.1440753718/(-0.464212743e-1*xa1+1)^2+1.219463331/(.4073000470+.5926999530*xa1)^2)+ln(1-xa1) = xb1^2*(.1440753718/(-0.464212743e-1*xb1+1)^2+1.219463331/(.4073000470+.5926999530*xb1)^2)+ln(1-xb1), (1-xa1)^2*(.4966874722/(.4073000470+.5926999530*xa1)^2+.1510891213/(-0.464212743e-1*xa1+1)^2)+ln(xa1) = (1-xb1)^2*(.4966874722/(.4073000470+.5926999530*xb1)^2+.1510891213/(-0.464212743e-1*xb1+1)^2)+ln(xb1)}, {xa1, xb1})$

(3)

>

Download NRTL.mw

How I can normalized mode shape for the comparison...

Asked by: 9009134 110

April 01 2020

0 7

Hello,

I obtained a mode shape from a vibration problem.

I want to normalized mode shape for the comparison of responses corresponding to different modes.

How I can normalize the mode shape that provided in the maple file?

The figure corresponds to this mode shape is plotted that is attached.

Thanks

mode_shapes.mw

Download mode_shapes.mw

What happen? sqrt()= same result...

Asked by: John Blacksad 10

March 27 2020

0 2

Hi guys, Maple always shows the same result when I use sqrt() function

How can I fix it? Thanks

How to write the sequence...

Asked by: brian bovril 904

March 26 2020

1 4

This is a simple model of the covid infection as it applies here, down under, without the use of any de's like

https://www.mapleprimes.com/posts/212100-Exploring-The-CoVid19-Outbreak

I'm told a substantial number of people do recover within 45 days ....
So on the 45th day 1 person has recovered, on day 50 two people have recovered and so on.

Just want to know how to write this sequence...

Covid_dbling.mw

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post