C_R

3412 Reputation

21 Badges

5 years, 316 days

MaplePrimes Activity


These are replies submitted by C_R

@sursumCorda 

I am sorry, I can't relate https://www.mapleprimes.com/questions/235749-Better-Way-To-Remove-Entries-In-Piecewise#answer292495 to this quesion and the transformed minus sign. I also cannot spot any substantial changes in the expressions in both pdf files.

@acer Thank you for answering all my questions and the procedure that reveals the name of the procedure that issued the warning. The procedure worked well on a few examples with the exception of e^2 and du/dy warnings. My interpretation is that those warnings are not using the WARNING function.

@sursumCorda 

printing (via printpreview) or exporting your file on my machine is still fine.

I guess the print error has something to do with the print driver (assuming that everyone is using Maple 2022. There have been substantial improvements with the introduction of the print layout mode). On my machine I used the Adobe print driver

confirmed_2.pdf

I can't reproduce this with Maple 2022.2. Which version do you use?

Untitled_(16).pdf

@WhiteNoise 

the complicated formating works in almost all instances. If for some reasons the user has defined a name kOmega (which is unlikely I have to admit) the output will not be as intended.

Can you provide Maple code that crashes?

I have an open question about volume integrals that I would like to at least ask without the need of getting an answer right away:
Is it possible to compute the moment on inertia about an arbitrary oriented axis? A stopped at the point where I could not manage to define the domain of a cylinder.

Maybe such an example fits better in the “Physics Courseware Support: Mechanics”.

@Thomas Richard 

Thank you!
With AllSolutions it becomes a question of equaltiy of integrals where I am stuck:
 

expr0 := int(1/(sqrt(x__0-x)*sqrt(-x^2+1)), x = 0 .. x__0, allsolutions); `assuming`([expr0], [0 < x__0 and x__0 < 1])

 

2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), I*(x__0+1)^(1/2)/(1-x__0)^(1/2))/(1-x__0)^(1/2)

(1)

Remove the imaginery unit by hand

expr1_new := 2*EllipticF(sqrt(x__0)/sqrt(x__0+1), -sqrt(x__0+1)/sqrt(x__0-1))/sqrt(1-x__0)

2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), (x__0+1)^(1/2)/(x__0-1)^(1/2))/(1-x__0)^(1/2)

(2)

expr2 := sqrt(2)*EllipticF(sqrt(1-1/(x__0+1))*sqrt(2), (1/2)*sqrt(x__0+1)*sqrt(2))

2^(1/2)*EllipticF((1-1/(x__0+1))^(1/2)*2^(1/2), (1/2)*(x__0+1)^(1/2)*2^(1/2))

(3)

plot(expr1_new-expr2, x__0 = 0 .. 1)

 

queryequal(expr1_new, expr2)

queryequal(2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), (x__0+1)^(1/2)/(x__0-1)^(1/2))/(1-x__0)^(1/2), 2^(1/2)*EllipticF((1-1/(x__0+1))^(1/2)*2^(1/2), (1/2)*(x__0+1)^(1/2)*2^(1/2)))

(4)

difference := op(1, queryequal(2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), (x__0+1)^(1/2)/(x__0-1)^(1/2))/(1-x__0)^(1/2), 2^(1/2)*EllipticF((1-1/(x__0+1))^(1/2)*2^(1/2), (1/2)*(x__0+1)^(1/2)*2^(1/2))))-op(2, queryequal(2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), (x__0+1)^(1/2)/(x__0-1)^(1/2))/(1-x__0)^(1/2), 2^(1/2)*EllipticF((1-1/(x__0+1))^(1/2)*2^(1/2), (1/2)*(x__0+1)^(1/2)*2^(1/2))))

2*EllipticF(x__0^(1/2)/(x__0+1)^(1/2), (x__0+1)^(1/2)/(x__0-1)^(1/2))/(1-x__0)^(1/2)-2^(1/2)*EllipticF((1-1/(x__0+1))^(1/2)*2^(1/2), (1/2)*(x__0+1)^(1/2)*2^(1/2))

(5)

`assuming`([convert(difference, Int)], [0 < x__0 and x__0 < 1])

2*(Int(1/((-_alpha1^2+1)^(1/2)*(1-(x__0+1)*_alpha1^2/(x__0-1))^(1/2)), _alpha1 = 0 .. x__0^(1/2)/(x__0+1)^(1/2)))/(1-x__0)^(1/2)-2^(1/2)*(Int(2/((-_alpha1^2+1)^(1/2)*(-2*_alpha1^2*x__0-2*_alpha1^2+4)^(1/2)), _alpha1 = 0 .. (1-1/(x__0+1))^(1/2)*2^(1/2)))

(6)

`assuming`([simplify(%)], [0 < x__0 and x__0 < 1])

-2*2^(1/2)*(Int(1/((-_alpha1^2+1)^(1/2)*(4+(-2-2*x__0)*_alpha1^2)^(1/2)), _alpha1 = 0 .. x__0^(1/2)*2^(1/2)/(x__0+1)^(1/2)))+2*(Int(1/((-_alpha1^2+1)^(1/2)*((x__0+1)*_alpha1^2+1-x__0)^(1/2)), _alpha1 = 0 .. x__0^(1/2)/(x__0+1)^(1/2)))

(7)

``

 

NULL

 

Can someone help?
 

Download Equality_of_Elliptic_expressions_02.mw

What you indend to do sounds intersting.

Could you give or refer to examples of what the final result you want might look like?

I didn't find anything useful that matched the keywords in your question.

@tomleslie 

Regrettably I forgot to mention that my simplified example is an excerpt from a larger worksheet where I took care of unit inconsistency risk mitigation (Testdimension in this case).

Your answer illustrates the step that I do for the time being: substituting s by a dimensionless s (s=s*Unit(1/s)) and doing the “simplifications”.

The result is an expression which is partially dimension less. To resolve this mixture of dimensionless expressions and expression with dimension, the dimensionless s has to be replaced by s=s*Unit(s). In summary, two substitutions I am looking to make expendable.

The reason why I have to restore a transfer function with units of s in Hertz is that I want to take advantage of recent plot functionality with units (the DynanmicSystems package does not provide this).

Off course an expand option is only one way. That’s why I added some other thoughts maybe worth considering the day an enhancement of the DynamicSystems package is planned.

Thanks for taking the time to highlight potential risks.

 

@dharr 

To complete the method you have given, it might be required to delete unnecessary characters (a colon in your example).

On a 4k display the box placeholders come with addtional information. I don't know if that is the case on 2k. fo could mean "function of". The multiplication box before the integration variable is interesting because mathematically there is no multiplication. However, when you drag an drop the integral (in the ouput of q) to a new input line (2 times), the new input shows that Maple correctly interprets the expression.

int(GAMMA(s + a)*GAMMA(b - s)*z^(-s), s = c - i*infinity .. c + i*infinity)

 

Overall, with a bit of caution errors and time can be saved with your method.
Thank you!

@AmirHosein Sadeghimanesh I have no idea. I think it should work without the style option.

It would be helpful if you upload your Maple worksheet. For that use the green arrow in the menu of the editor.

@acer Thank you for the example. It perfectly shows the differences.

@EngM 

Have you considered using a C code block? Technically, you could export Maple expressions to C code once the results of your matrix calculation are derived in an explicit form so that they can become the outputs of a C code block.

In the image above I have recently done this for large expressions. Since the sequence of essential steps and commands is hard to remember I have summarized them here for myself and others.

Consider asking how your matrix calculations can be done in a separate question (preferably with attached example files). It's and interesting one. It’s not excluded that there are other solutions.

Also your question about connecting components inside a Modelica component (and @WhiteNoise's answer) would have been worth branching off to make the content easier to find. 

 

First 53 54 55 56 57 58 59 Last Page 55 of 67