Thomas Richard

Mr. Thomas Richard

3496 Reputation

13 Badges

14 years, 134 days
Maplesoft Europe GmbH
Technical professional in industry or government
Aachen, North Rhine-Westphalia, Germany

MaplePrimes Activity


These are answers submitted by Thomas Richard

The two linear equations you feed to solve are contradictory for f2. Check this in a separate step:

lineq1 := -2*A*xc - B*yc = D; lineq2 := -B*xc - 2*C*yc = E;
4*lineq1 + 3*lineq2;

For f1, they have a solution.

The use ... end use construct is newer than those control structures affected by interface(longdelim). One could think about introducing end as a short form for end use, but that would be counterproductive, as the recommendation is to use the long form for better readability.

I do not fully understand your question. You really should upload your worksheet to make it easier for others to help.

Anyway, here's my attempt. The algsubs command won't work because of the presence of sqrt. Fortunately, we can use subs here (simpler syntactic substitution), and by providing an assumption on c, we can simplify the result.

Please elaborate on the remaining steps of your draft.

# https://www.mapleprimes.com/questions/238930-How-To-Do-The-Algebraic-Extension-In

restart:

U := -2*a[0]*((x^3*a[1] + 42)*diff(y(x), x, x) - 6*y(x)*x*a[1])*sqrt(-a[0]*(a[1]^2 + 3)) - 3*x^2*((1/3*a[1]^3 + a[0] + a[1])*x + (8*a[1]^2 + 24)*a[0])*diff(y(x), x, x);

-2*a[0]*((x^3*a[1]+42)*(diff(diff(y(x), x), x))-6*y(x)*x*a[1])*(-a[0]*(a[1]^2+3))^(1/2)-3*x^2*(((1/3)*a[1]^3+a[0]+a[1])*x+(8*a[1]^2+24)*a[0])*(diff(diff(y(x), x), x))

(1)

V := diff(U,x);

-2*a[0]*(3*x^2*a[1]*(diff(diff(y(x), x), x))+(x^3*a[1]+42)*(diff(diff(diff(y(x), x), x), x))-6*(diff(y(x), x))*x*a[1]-6*y(x)*a[1])*(-a[0]*(a[1]^2+3))^(1/2)-6*x*(((1/3)*a[1]^3+a[0]+a[1])*x+(8*a[1]^2+24)*a[0])*(diff(diff(y(x), x), x))-3*x^2*((1/3)*a[1]^3+a[0]+a[1])*(diff(diff(y(x), x), x))-3*x^2*(((1/3)*a[1]^3+a[0]+a[1])*x+(8*a[1]^2+24)*a[0])*(diff(diff(diff(y(x), x), x), x))

(2)

# resultants:

R0 := resultant(U,V,a[0]);

Error, (in resultant) invalid arguments

 

R1 := resultant(U,V,a[1]);

Error, (in resultant) invalid arguments

 

#

S := simplify(subs((-a[0]*(a[1]^2+3))=c^2,U)) assuming c>=0;

((-2*c*a[0]*a[1]-a[1]^3-3*a[0]-3*a[1])*x^3+(-24*a[1]^2-72)*a[0]*x^2-84*a[0]*c)*(diff(diff(y(x), x), x))+12*y(x)*c*x*a[0]*a[1]

(3)

#

 

 

Download ode-parameter-elimination.mw

Just supply the output=DataFrame option to the ExcelTools:-Import command. The default type is Matrix.

Please see the ?_ (or ?underscore) help page.

The command for grouping multiple terms into a common fraction is normal:

normal(e);

 

The defaults for scaletorange haven been improved in Maple 2023.

Independently of that, I suggest assigning A:=abs(j) (much shorter) and inserting another plot option: scaling=constrained (otherwise the circle looks too elliptic).

restart:

with(plots):

a := 1/2:

j := BesselJ(0,sqrt(2*I)*sqrt(x^2+y^2)*a)/BesselJ(0,sqrt(2*I)*a);

BesselJ(0, (1/2+(1/2)*I)*(x^2+y^2)^(1/2))/BesselJ(0, 1/2+(1/2)*I)

(1)

A := abs(j); # shorter than A := sqrt(Re(j)^2 + Im(j)^2);

abs(BesselJ(0, (1/2+(1/2)*I)*(x^2+y^2)^(1/2))/BesselJ(0, 1/2+(1/2)*I))

(2)

densityplot(A, x = -1 .. 1, y = -sqrt(-x^2 + 1) .. sqrt(-x^2 + 1), grid = [68, 68], style = patchnogrid, scaling=constrained, scaletorange = 0 .. 1);

 

 


 

Download scaletorange-M23M24.mw

 

Color functions are available for 3D plots only.

This is fairly simple if we use the Matrix constructor along with your mapping (shown here in 1D Math aka Maple notation):

`ΔPV1` := Matrix(12, 8, (i, j) -> 100*(100*H0[i] - 100*H1[i,j])/((100 - H0[i])*H1[i,j]))

See the modified worksheet. I deleted some stuff at the end - perhaps you still need that, sorry. Also, I had to remove the absolute path from the Import call to make it work on my machine.

QuestionPrimes_TR.mw

Since your worksheet contains no dsolve call, it is unclear where your solution is coming from, and what assumptions on parameters are needed. However, the odetest result can be simplified as follows:

ot2 := odetest(G1, pde2);
simplify(ot2,'symbolic');

Note that the 'symbolic' option must be handled with care, as it may not be valid for all parameter values in the complex domain.

If we assume(x0>1) then it works in Maple 2024.1 as well. Perhaps the solution process gets stuck without this information.

According to the documentation, the Crout algorithm is for tridiagonal matrices only. Perhaps you have some other definition (or naming convention)? Anyway, the result can be obtained by the MatrixDecomposition command of the Student:-NumericalAnalysis package, using method=LDU option. Worksheet attached:

Download MatrixDecompositionLDU.mw

You can verify the IVP solution as follows (using Maple 2024.1 here, currently with no Physics updates):

simplify(combine(residual)) assuming positive;

It will return [0,0].

I don't think proper fading is supported, but if I interpret your question correctly, it is sufficient to let both endpoints depend on the animation parameter:

plots[animate]( plot, [sin(t), t=x..x+Pi], x=0..4*Pi );

That should apply to space curves as well, but I did not actually try (today is very busy for me, sorry).

This is the expected behaviour of expand and combine for trig expressions. Since you started with the latter command, take a look at the ?combine help page, notably the 2nd bullet item of the Description section. The subpage ?combine,trig has more details.

In your example, expand returns a polynomial in cos(x). [You may omit the trig option, by the way.] In this particular sense, it is the inverse operation of combine which yields a cos of a multiple of x as your input expression. For ordinary polynomials, the inverse of expand is factor.

A different question is which representation is simpler - that's debatable, of course.

I cannot say much about the Context Panel, as I'm rarely ever using it.

1 2 3 4 5 6 7 Last Page 1 of 45