Thomas Richard

Mr. Thomas Richard

3556 Reputation

13 Badges

14 years, 160 days
Technical professional in industry or government
Aachen, North Rhine-Westphalia, Germany

MaplePrimes Activity


These are answers submitted by Thomas Richard

This ODE is beyond the scope of the Student package. You get a better error message when expanding the product first:

Student:-ODEs:-ODESteps(expand(ode));
Error, (in Student:-ODEs:-OdeSolveOrder2) ODE is not supported

 

Check Tools -> Options -> Interface -> Large toolbar icons (first checkbox in the "Show" section).

I guess you want a space curve as follows:

plots:-spacecurve([t*sin(t), t*cos(t), 2/3*sqrt(2)*t^(3/2)], t=0..Pi);

To reproduce the graph from your screenshot, you will have to rotate the output a bit, or add options accordingly.

I'm too lazy for that. ;-)

By the way, since you asked for the length:

AL := VectorCalculus:-ArcLength(<t*sin(t), t*cos(t), 2/3*sqrt(2)*t^(3/2)>, t=0..Pi);

 

This is a good example for Maple's solve/identity feature:

restart:

poly := x^5 - 5*x^4 - 35*x^3 + a*x^2 + b*x + c; # partially given polynomial

x^5-5*x^4+a*x^2-35*x^3+b*x+c

(1)

fpoly := (x-r)*(x-r-d)*(x-r-2*d)*(x-r-3*d)*(x-r-4*d); # roots in arithmetic progression

(x-r)*(x-r-d)*(x-r-2*d)*(x-r-3*d)*(x-r-4*d)

(2)

efpoly := expand(fpoly);

-24*d^4*r+24*d^4*x-50*d^3*r^2+100*d^3*r*x-50*d^3*x^2-35*d^2*r^3+105*d^2*r^2*x-105*d^2*r*x^2+35*d^2*x^3-10*d*r^4+40*d*r^3*x-60*d*r^2*x^2+40*d*r*x^3-10*d*x^4-r^5+5*r^4*x-10*r^3*x^2+10*r^2*x^3-5*r*x^4+x^5

(3)

# see ?solve,identity

sol := solve(identity(poly = efpoly, x));

{a = 125, b = 194, c = -280, d = -3, r = 7}, {a = 125, b = 194, c = -280, d = 3, r = -5}

(4)

 


Download: roots-wanted.mw

There is no such thing as a Purchase Code for the Maple Calculator. If you want to upgrade from the free version to Premium, that's easily done from within the app.

Perhaps you mean something different?

If your question is about the purchase as such, please contact our Customer Service.

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):

`&Delta;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.

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