MaplePrimes Questions

The netlist in the attached Maple file will not solve. 
It will solve if I remove

.CONS Vx=v[D]-v[S]
.VARS Vx

from the netlist.

Syrup returns a 2 element list of [eqs, vars] that can be easily solved using the solve command as shown in the attached file

ModelDebug.mw

It will be great if  DEtools:-convertAlg  which is designed to find coefficients of a linear ode, not to expand the coefficients it finds.

This caused my patterns to fail, as after putting the coefficients in a certain form, and then calling DEtools:-convertAlg later on somewhere else in the program, the patterns no longer worked due to this unexpected expansion.

It will be better if DEtools:-convertAlg returns the coefficients it find as is and not expand them. I could not find an option to use for this. Here is an example.

restart;
ode:=b*(x+x^2)*diff(y(x),x$2)+x*(2*a - 1)*diff(y(x),x)+c*y(x)=0;

After the call the coefficients become

L:=DEtools:-convertAlg(ode,y(x))

You might ask that one can easily "fix" this afterwords. Sure. I can use collect on x. But the point is that in general, it might not be clear what to collect on or what the changes that were made are.

In principle, it is better to keep the coefficients in same form as they were in the ode so that the ode do not change form.

Maple 2022.1 on windows.

For example my friend emailed me.
His email address is oneman@gmail.com, the email content is HELLO.
How does Maple print HELLO to the screen so I can see it?
Thanks for your help!

restart

eq1 := (2*(r^2+a^2*cos(theta)^2))*(M*r-(1/2)*a^2-(1/2)*r^2)*(diff(f(r, theta), r, theta))+(2*(a^2*(M-r)*cos(theta)^2-M*r^2+a^2*r))*(diff(f(r, theta), theta))

2*(r^2+a^2*cos(theta)^2)*(M*r-(1/2)*a^2-(1/2)*r^2)*(diff(diff(f(r, theta), r), theta))+2*(a^2*(M-r)*cos(theta)^2-M*r^2+a^2*r)*(diff(f(r, theta), theta))

(1)

eq2 := sin(theta)*(r^2+a^2*cos(theta)^2)*(diff(f(r, theta), theta, theta))-cos(theta)*(diff(f(r, theta), theta))*(a^2*cos(theta)^2-2*a^2-r^2)

sin(theta)*(r^2+a^2*cos(theta)^2)*(diff(diff(f(r, theta), theta), theta))-cos(theta)*(diff(f(r, theta), theta))*(a^2*cos(theta)^2-2*a^2-r^2)

(2)

eq3 := -2*(r^2+a^2*cos(theta)^2)^2*(M*r-(1/2)*a^2-(1/2)*r^2)*sin(theta)*(diff(g(r, theta), r, r))+sin(theta)*(r^2+a^2*cos(theta)^2)^2*(diff(g(r, theta), theta, theta))+(4*(-(1/4)*cos(theta)^4*a^4+a^2*r*(M-(1/2)*r)*cos(theta)^2-M*a^2*r-(1/4)*r^4))*cos(theta)*(diff(g(r, theta), theta))-2*M*sin(theta)*(diff(g(r, theta), r))*(a^2+r^2)*(cos(theta)*a-r)*(cos(theta)*a+r)

-2*(r^2+a^2*cos(theta)^2)^2*(M*r-(1/2)*a^2-(1/2)*r^2)*sin(theta)*(diff(diff(g(r, theta), r), r))+sin(theta)*(r^2+a^2*cos(theta)^2)^2*(diff(diff(g(r, theta), theta), theta))+4*(-(1/4)*cos(theta)^4*a^4+a^2*r*(M-(1/2)*r)*cos(theta)^2-M*a^2*r-(1/4)*r^4)*cos(theta)*(diff(g(r, theta), theta))-2*M*sin(theta)*(diff(g(r, theta), r))*(a^2+r^2)*(cos(theta)*a-r)*(cos(theta)*a+r)

(3)

pdsolve([eq1, eq2, eq3])

[{f(r, theta) = _F1(r)+(Int((r^2+a^2*cos(theta)^2)/((cos(theta)+1)^(1/2)*(cos(theta)-1)^(1/2)), theta))*_C1/(2*M*r-a^2-r^2)}, [diff(diff(g(r, theta), r), r) = (sin(theta)*(sin(theta)^2*a^2-a^2-r^2)^2*(diff(diff(g(r, theta), theta), theta))-4*cos(theta)*((1/4)*a^4*sin(theta)^4+(-(1/2)*a^2+r*(M-(1/2)*r))*a^2*sin(theta)^2+(1/4)*(a^2+r^2)^2)*(diff(g(r, theta), theta))+2*M*sin(theta)*(diff(g(r, theta), r))*(a^2+r^2)*(sin(theta)^2*a^2-a^2+r^2))/(sin(theta)*(sin(theta)^2*a^2-a^2-r^2)^2*(2*M*r-a^2-r^2))]]

(4)

``

Download pde1.mw

 

restart

A := Omega; assume(Omega > 0); B = Omega

B = Omega

(1)

``

Download assumetest.mw

``

In a physics problem, I came across the following triple integral:

exp := -sin(alpha)*i*r*(-sin(alpha)*cos(phi)*cos(theta)+sin(theta)*cos(alpha))/(4*sqrt(-2*sin(theta)*sin(alpha)*cos(phi)*p*r-2*cos(alpha)*cos(theta)*p*r+p^2+r^2)*Pi(-2*sin(theta)*sin(alpha)*cos(phi)*p*r-2*cos(alpha)*cos(theta)*p*r+p^2+r^2)*(-2+sqrt(2))*Pi)

`assuming`([int(int(int(exp*p^2*sin(alpha), p = 0 .. 1), alpha = 0 .. (1/4)*Pi), phi = 0 .. 2*Pi)], [alpha > 0, alpha < (1/4)*Pi, r > 0, r < 1, phi > 0, phi < 2*Pi, theta > 0, theta < (1/4)*Pi])

``

I tried to perform each integration separately, but got no result. Is there any transformation or procedure, that I am not aware of, to accomplish this task?
Grateful,
Oliveira

Download Triple_integral.mw

Hello,

   A question on programming Procedures (proc).

   What is the best, or most common method, to invoke packages within proc or at the worksheet level?

   I have always been confused about Proc scope for packages.   I use Scientific Constants a lot.   It and several other "with(package)" are part my main loaded pakages right after restart command.   LinearAlgebra is another common used package.

   It has been a question for me of whether to implement a "with()" statement inside of the proc to bound the scope or ensure it is seen/available by the procedure for the available commands inside the packages.  

  All thoughts appreciated.

Thanks,
Bill

After substitution of (10) into (4), how to collect the terms of like powers of eta (i.e., eta^-3, eta^-2,eta^-1, eta^0, eta^1,eta^2 ), and equate the coefficients to
zero, get a system of algebraic equations for A[m]?

 

PA.mw 

g := proc() use a=a+b in use b= a-b in a*b end use end use end proc:

Maple shows g := proc() (a+b)*(a+b) end proc. 

My guess is (a+b)*a. What is wrong with me?

I faced a problem with the following syntax which provided me we partial circles inequality plotting, and how able to improve the syntax:
restart: with(plots):
inequal({y > -1, y >= x^2+1, (x-1)^2+(y-1)^2 <= 16}, x = -5 .. 8, y = -6 .. 6, optionsfeasible = (color = grey), optionsexcluded = (color = white));

restart: with(plots): inequal({y > -1, y >= x^2+1, (x-1)^2+(y-1)^2 <= 16}, x = -5 .. 8, y = -6 .. 6, optionsfeasible = (color = grey), optionsexcluded = (color = white));

Hello, 

How can I force Maple to perform division on a fractional polynomial?  Here's an exmaple of what Im trying to do - I want rr and nsr to be divided and simplified and return a polynomial.

I can see from this that its doing what i want for integer exponents - just not fractions!

I've tried everything I can think of... Let me know! Thanks.

Hello! 

Ive stubled into something odd. Here it is: 

Now I really wonder how they figured how i was suppose to figure out how this is suppose to go. They did not explain all that much on how to find this out. This specifically is a "Computerized Question." 

It says "find the smallest number "n", so that A^n=I"

This was the result when i tried to solve it:

"Error, (in Engine:-Dispatch)" I have no idea what is wrong. 

Any way, 

Greetings,

The Function 

#Opdracht 9

A := `<,>`(`<|>`(0, 1, 0, 0), `<|>`(0, 0, 1, 0), `<|>`(1, 0, 0, 0), `<|>`(0, 0, 0, -1))

Matrix(%id = 18446746584512210574)

(1)

solve(A^n = I, n)

Error, (in Engine:-Dispatch) badly formed input to solve: not fully algebraic

 

``

Download Mapleprimes_Question_Book_2_Paragraph_4.1_Question_9.mw

I make new cone puzzle. However, I cat't make function l(θ). Can maple solve this puzzle?

θ=90 degree is a YouTube problem I found. 

 

Have to draw the graph

 

What's going on in the following? Why can't I restore the default behavior of diff after using Physics:-diff or even just using with(Physics)? Is it because of Physics:-ModuleLoad()?

restart:

diff(conjugate(f(x)), x);

(diff(f(x), x))*(-conjugate(f(x))/f(x)+2*abs(1, f(x))/signum(f(x)))

Physics:-diff(conjugate(f(x)), x);

diff(conjugate(f(x)), x)

forget(.., conjugate(f(x)));

diff(conjugate(f(x)), x);

diff(conjugate(f(x)), x)

restart:

with(Physics):

:-diff(:-conjugate(f(x)), x);

diff(conjugate(f(x)), x)

 

Download PhysicsDiff.mw

I don't know why the prettyprinted output of my worksheet is shown the way that it is above. I didn't do anything differently than I usually do to upload a worksheet. Anyway, the output is simple enough that I think that my Question is still clear.

First 273 274 275 276 277 278 279 Last Page 275 of 2425