MaplePrimes Questions

i want to plot form 15 to 0 and 5 to 1.....but my plot shows only increasing order..how can i change it?

Download plot.mw

Is there a way to suppress dependencies of a all functions appearing in a worksheet without having to do suppress([each function])?

out_parameter_for_Sample(_)_and_rtables.mw
When I generate sample data using a random variable, a normal distribution and the Sample command, I want to see the sample and the values used to randomly evaluate the distribution which generates said sample. The Maple help sheets say you can so this with the out parameter in the Sample command. I try several different calls and hoped the rtable is filled on execution. As seen below, the execution fails. Am a specifying the rtable incorrectly? See attached Maple sheet.

restart

NULL

NULL

with(Statistics)

with(LinearAlgebra)

NULL

When I generate sample data using a random variable, a normal distribution and the Sample commands, I want to see the sample and the values used to randomly evaluate the distribution. The Maple help sheets say you can so this with the out parameter in the Sample command.

NULL

NULL

NULL

From the help, we are required to specify a float rtable, that will be filled when the Sample command is succesfully executed.

NULL

Define a random variable.

X := RandomVariable(Normal(0, 1))

_R

(1)

Sample that random variable. We sample it 10 times and look at the output.

A := Sample(X, 10)

Vector[row](%id = 36893491057099280004)

(2)

The float rtable is defined with the sample number of elements as will be sampled.

tab1 := rtable(1 .. 10, datatype = float)

Array(%id = 36893491057099270012)

(3)

NULL

tab2 := convert(rtable(1 .. 10, datatype = float), list)

[HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0), HFloat(0.0)]

(4)

Now I try several different calls and hope the rtable is filled on execution.

A := Sample(X, 10, out)

Error, invalid input: no implementation of Statistics:-Sample:-ModuleApply matches the arguments in call, 'Statistics:-Sample:-ModuleApply(X,10,out)'

 

A := Sample(X, 10, out = tab1)

Error, invalid input: no implementation of Statistics:-Sample:-ModuleApply matches the arguments in call, 'Statistics:-Sample:-ModuleApply(X,10,out = tab1)'

 

A := Sample(X, 10, out = tab2)

Error, invalid input: no implementation of Statistics:-Sample:-ModuleApply matches the arguments in call, 'Statistics:-Sample:-ModuleApply(X,10,out = tab2)'

 

A := Sample(X, 10, out = true)

Error, invalid input: no implementation of Statistics:-Sample:-ModuleApply matches the arguments in call, 'Statistics:-Sample:-ModuleApply(X,10,out = true)'

 

NULL

As seen above, the execution fails. Am a specifying the rtable incorrectly?

NULL

NULL

NULL


 

Download out_parameter_for_Sample(_)_and_rtables.mw

 

I followed the code on the website https://de.maplesoft.com/support/help/maple/view.aspx?path=updates/Maple18/GraphTheory to convert a graph to LaTeX code. However, after compiling with pdflatex, I found that some edges of the graph are jagged.

restart:    
with(GraphTheory):
with(SpecialGraphs):
S:=SoccerBallGraph():
Latex(S,FileTools:-JoinPath([currentdir(), "soccer.tex"]),300,300,true)

soccer.pdf

I suspect it's because of the converted LaTeX code.

PS: The PDF conversion issue from last time still remains unsolved in Maple 2023; see

https://www.mapleprimes.com/questions/236142-How-To-Remove-The-Mosaic-Of-Vertices.

How to rectify this Error,in RK Method.Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 3, got 2.

I cound't plot p5,p6,p7.

If RK Method is suitable for this or not please tell the suitable numerical method code for this.Help me

IP-TEMP.mw

The example worksheet uses table to determine the covering relations in a POSET. This is old (Maple 11 and earlier code) that I have been trying to update into a package to explore calculations in an algebraic structure.  The example procedure seems to work well, but because of problems with similar procedures elsewhere, I have some concerns about the validity of using tables in this fashion. Specifically, this procedure initializes a table, then proceeds to modify the table entries, and then reformats the sequences into sets.  I have not been able to find documentation for modifying entries in tables after they have been defined. The documentation for tables only covers adding entries, removing entries, but not modifying entries.

Is modifying tables as my  procedure does an undocumented feature?

In addition, the documentation does not explain how to clear a table. It only describes how to clear a table entry. Older code sometimes purported to clear a table by assigning its name (with uneval quotes) to itself, but this does not seem to work.

exampletableoperations.mw

Hi to All! I have a big expression including 10 constant parameters (k_1, k_2, k_3, omega_1, omega_2, omega_3, a_3, b_2, c_0, c_2) to be solved for a_3, b_2, c_0, c_2. I wonder if there is any relation between the constants a_3, b_2, c_0, c_2 not depending on (k_1, k_2, k_3, omega_1, omega_2, omega_3) satisfying this big expression e.g. a_3=2, b_2=c_0/5, c_2 free. 

I also tried fsolve but I could not get an answer.

Thank you in advance.

Regards,

question1.mw

In accordance with this statement obtained by Чебышёв (1853), each of 

simplify(int(x^(1/2)*(x^2 + 1)^(-3/4), x), symbolic);
simplify(int((x^(1)*(1 - x^2))^(1/3), x), symnolic);
simplify(int(x^(-1)*(x^6 + 1)^(-1/6), x), symnolic);
simplify(int(x^(17/2)*(x^2 + 1)^(1/4), x), symnolic);

can be reduced to an integral of rational functions, which can be expressed in terms of elementary functions. But it appears that Maple 2023.0 is still unable to completely calculate them. For instance: 
 

restart;

interface(version)

`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`

(1)

timelimit(0.1e4, `assuming`([simplify(int(x^(1/2)/(x^2+1)^(3/4), x))], [x > 0]))

(2/3)*x^(3/2)*hypergeom([3/4, 3/4], [7/4], -x^2)

(2)

timelimit(0.1e4, `assuming`([simplify(int((x*(-x^2+1))^(1/3), x))], [`or`(`and`(x <= 1, x >= 0), x <= -1)]))

(3/4)*x^(4/3)*hypergeom([-1/3, 2/3], [5/3], x^2)

(3)

timelimit(0.1e4, `assuming`([simplify(int(1/(x*(x^6+1)^(1/6)), x))], [x <> 0]))

-(1/36)*x^6*hypergeom([1, 1, 7/6], [2, 2], -x^6)-(1/12)*Pi*3^(1/2)-(1/3)*ln(2)-(1/4)*ln(3)+ln(x)

(4)

timelimit(0.1e4, `assuming`([simplify(int(x^(17/2)*(x^2+1)^(1/4), x))], [x > 0]))

(1/81920)*(8192*x^10*(x^2+1)^(1/4)+512*x^8*(x^2+1)^(1/4)-640*x^6*(x^2+1)^(1/4)+880*x^4*(x^2+1)^(1/4)-1155*ln(2*RootOf(_Z^2+1)*x^(1/2)*(x^2+1)^(3/4)-2*RootOf(_Z^2+1)*x^(3/2)*(x^2+1)^(1/4)+2*(x^2+1)^(1/2)*x-2*x^2-1)*x^(1/2)*RootOf(_Z^2+1)-1540*x^2*(x^2+1)^(1/4)-1155*x^(1/2)*ln(-2*x^(1/2)*(x^2+1)^(3/4)-2*x^(3/2)*(x^2+1)^(1/4)+2*(x^2+1)^(1/2)*x+2*x^2+1))/x^(1/2)

(5)

gc()


 

Download Chebyshev_theorem_on_the_integration_of_binomial_differentials.mw

However, closed-form (and readable) solutions in elementary forms exist (cf. Regression reports for Computer Algebra Independent Integration Tests. Summer 2022 version (12000.org)); in fact, Mathematica returns: 

So, why can't Maple find these compact antiderivatives (expressed by elementary functions) directly here? In other words, is there a way to resolve them in Maple without applying some change of the variable to these indefinite integrals manually?

Trying to solve in Maple:

restart;
f := 15;
fk := 7;
zm := 350;
ym := 200;
eps := 1 - fk^2*exp(-((z - zm)/ym)^2)/f^2;
dz := diff(z(x, bn), bn);
db := diff(b(x, bn), bn);
eq1 := diff(z(x, bn), x) = cot(b(x, bn));
eq2 := diff(b(x, bn), x) = subs(z = z(x, bn), -1/(2*eps)*diff(eps, z));
eq3 := diff(dz, x) = -bd/sin(b(x, bn))^2;
eq4 := diff(bd, x) = subs(z = z(x, bn), dz/(2*eps)*(diff(eps, z)^2/eps - diff(eps, z $ 2)));
sys := eq1, eq2, eq3, eq4;
cond := z(0, bn) = 0, b(0, bn) = bn, zd(0, bn) = 0, bd(0, bn) = 1;
dsolve({cond, sys}, [z(x, bn), b(x, bn)], numeric);
Error, (in dsolve/numeric/process_input) dependent variables must be functions of a single unknown, the independent variable. Got [z(x, bn), b(x, bn)]

What does "Error, (in dsolve/numeric/process_input) dependent variables must be functions of a single unknown, the independent variable. Got [z(x, bn), b(x, bn)] " mean?

How can this system be solved?

So that computed results do not take up many lines in Euler transforms and in various tensor math, how can I present the sine and cosine functions in Maple with this compact form:

  • sin(x(t)) = s(x) or sx
  • cos(x(t)) = c(x) or cx

The input and computed output of symbolic calculations equations needs to be output in this compact notation. Derivatives with the compact notation are understood to be function of time.

Here is an example implimentation in Mathematica:

Matricies Rx and Ry are similarly defined and the dot product can be computed:

Here's example derivative:

Any help or pointers on this is appreciated. I'm new on Maple Primes and am not sure how to find out if this is already posted somewhere.

Thanks,

David

The issue arises from solving the following ODEs in Maple (where a is a suitable real parameter): 

ode__1 := a*(diff(y(x), x) + 1)^2 + (y(x) - x)^2*diff(y(x), x) = 0: # dsolve(ode__1);
ode__4 := a*(x*diff(y(x), x) + y(x))^2 - (y(x) + x)^2*diff(y(x), x) = 0: # dsolve(ode__4);

However, dsolve cannot give fully simplified solutions, so I have to compute these unevaluated integrals (i.e., expr1) manually: (For the sake of completeness, I list some related ODEs below.) 
 

restart;

ode__1 := a*(diff(y(x), x)+1)^2+(y(x)-x)^2*(diff(y(x), x)) = 0
ode__4 := a*(x*(diff(y(x), x))+y(x))^2-(y(x)+x)^2*(diff(y(x), x)) = 0

dsolve(ode__1, y(x), explicit)

expr__1 := convert(useInt(eval(selectfun([%], Intat), _Z = z)), list)

[Int(1/(z^2+(z^4+4*a*z^2)^(1/2)+4*a), z), Int(-1/(z^2-(z^4+4*a*z^2)^(1/2)+4*a), z)]

(1)

value(expr__1)

[(z^4+4*a*z^2)^(1/2)*((1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)+2*(-a)^(1/2)*ln(z+((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)-2*(-a)^(1/2)*ln(z+((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)-2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)^2*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2))*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z*(z^2+4*a)^(1/2))+(1/2)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/2)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(1/2))-4*a*(-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*a*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/8)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(3/2))), -(1/2)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))+(1/2)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(1/2))+(z^4+4*a*z^2)^(1/2)*((1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)+2*(-a)^(1/2)*ln(z+((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)-2*(-a)^(1/2)*ln(z+((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)-2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)^2*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2))*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z*(z^2+4*a)^(1/2))+4*a*(-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*a*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/8)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(3/2)))]

(2)

verify(diff([-z/(z^2+sqrt(z^2*(z^2+4*a))), z/(z^2-sqrt(z^2*(z^2+4*a)))], z), `~`[op](1, expr__1), simplify)

true

(3)

dsolve(ode__4, y(x), explicit)

expr__4 := convert(useInt(eval(selectfun([%], Intat), _Z = z)), list)

[Int((z^2-4*a*z+(-4*a*z^3+z^4-8*a*z^2+4*z^3-4*a*z+6*z^2+4*z+1)^(1/2)+2*z+1)/(z*(-4*a*z+z^2+2*z+1)), z), Int(-(z^2-4*a*z+2*z+1-((-4*a*z+z^2+2*z+1)*(z+1)^2)^(1/2))/(z*(-4*a*z+z^2+2*z+1)), z)]

(4)

value(expr__4)

[(-4*a*z^3+z^4-8*a*z^2+4*z^3-4*a*z+6*z^2+4*z+1)^(1/2)*(ln(z-2*a+1+(-4*a*z+z^2+2*z+1)^(1/2))+arctanh((2*a*z-z-1)/(-4*a*z+z^2+2*z+1)^(1/2)))/((z+1)*(-4*a*z+z^2+2*z+1)^(1/2))+ln(z), ((-4*a*z+z^2+2*z+1)*(z+1)^2)^(1/2)*(ln(z-2*a+1+(-4*a*z+z^2+2*z+1)^(1/2))+arctanh((2*a*z-z-1)/(-4*a*z+z^2+2*z+1)^(1/2)))/((z+1)*(-4*a*z+z^2+2*z+1)^(1/2))-ln(z)]

(5)

verify(diff([2*arctanh(sqrt((z+1)^2*(z*(z-2*(2*a-1))+1))/(z^2-1))+ln(z), 2*arctanh(sqrt((z+1)^2*(z*(z-2*(2*a-1))+1))/(z^2-1))-ln(z)], z), `~`[op](1, expr__4), simplify)

true

(6)

NULL


 

Download senseless_results_of_int.mw
 

restart;

ode__1 := a*(diff(y(x), x)+1)^2+(y(x)-x)^2*(diff(y(x), x)) = 0
ode__4 := a*(x*(diff(y(x), x))+y(x))^2-(y(x)+x)^2*(diff(y(x), x)) = 0

dsolve(ode__1, y(x), explicit)

expr__1 := convert(useInt(eval(selectfun([%], Intat), _Z = z)), list)

[Int(1/(z^2+(z^4+4*a*z^2)^(1/2)+4*a), z), Int(-1/(z^2-(z^4+4*a*z^2)^(1/2)+4*a), z)]

(1)

value(expr__1)

[(z^4+4*a*z^2)^(1/2)*((1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)+2*(-a)^(1/2)*ln(z+((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)-2*(-a)^(1/2)*ln(z+((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)-2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)^2*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2))*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z*(z^2+4*a)^(1/2))+(1/2)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/2)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(1/2))-4*a*(-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*a*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/8)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(3/2))), -(1/2)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))+(1/2)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(1/2))+(z^4+4*a*z^2)^(1/2)*((1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)+2*(-a)^(1/2)*ln(z+((z-2*(-a)^(1/2))^2+4*(-a)^(1/2)*(z-2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)*(((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)-2*(-a)^(1/2)*ln(z+((z+2*(-a)^(1/2))^2-4*(-a)^(1/2)*(z+2*(-a)^(1/2)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))+(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)-2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^2-4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2)))-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)^2*(((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*ln(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)-2*(z^4+4*a*z^2)*a*ln((8*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)*((z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))^2+4*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)*(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1))/((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)+4*(z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2))/(z+2*(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)/(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*((z^4+4*a*z^2)*a/(z^2*(z^2+4*a)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)))^(1/2)))/((-(-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)+(-a)^(1/2))*((-a)^(1/2)*(z^4+4*a*z^2)/(z^2*(z^2+4*a))+(a*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))-1)*((z^4+4*a*z^2)^(1/2)/(z*(z^2+4*a)^(1/2))+1))^(1/2)-(-a)^(1/2))*(-(z^4+4*a*z^2)/(z^2*(z^2+4*a))+1)))/(z*(z^2+4*a)^(1/2))+4*a*(-(1/8)*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)*z^2*(z^2+4*a)*arctanh((1/2)*z*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1)/(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))/((z^4+4*a*z^2)*a*(a*((z^4+4*a*z^2)/(z^2*(z^2+4*a))-1))^(1/2))-(1/8)*z^2*(z^2+4*a)*arctan((1/2)*z/a^(1/2))/((z^4+4*a*z^2)*a^(3/2)))]

(2)

verify(diff([-z/(z^2+sqrt(z^2*(z^2+4*a))), z/(z^2-sqrt(z^2*(z^2+4*a)))], z), `~`[op](1, expr__1), simplify)

true

(3)

dsolve(ode__4, y(x), explicit)

expr__4 := convert(useInt(eval(selectfun([%], Intat), _Z = z)), list)

[Int((z^2-4*a*z+(-4*a*z^3+z^4-8*a*z^2+4*z^3-4*a*z+6*z^2+4*z+1)^(1/2)+2*z+1)/(z*(-4*a*z+z^2+2*z+1)), z), Int(-(z^2-4*a*z+2*z+1-((-4*a*z+z^2+2*z+1)*(z+1)^2)^(1/2))/(z*(-4*a*z+z^2+2*z+1)), z)]

(4)

value(expr__4)

[(-4*a*z^3+z^4-8*a*z^2+4*z^3-4*a*z+6*z^2+4*z+1)^(1/2)*(ln(z-2*a+1+(-4*a*z+z^2+2*z+1)^(1/2))+arctanh((2*a*z-z-1)/(-4*a*z+z^2+2*z+1)^(1/2)))/((z+1)*(-4*a*z+z^2+2*z+1)^(1/2))+ln(z), ((-4*a*z+z^2+2*z+1)*(z+1)^2)^(1/2)*(ln(z-2*a+1+(-4*a*z+z^2+2*z+1)^(1/2))+arctanh((2*a*z-z-1)/(-4*a*z+z^2+2*z+1)^(1/2)))/((z+1)*(-4*a*z+z^2+2*z+1)^(1/2))-ln(z)]

(5)

verify(diff([2*arctanh(sqrt((z+1)^2*(z*(z-2*(2*a-1))+1))/(z^2-1))+ln(z), 2*arctanh(sqrt((z+1)^2*(z*(z-2*(2*a-1))+1))/(z^2-1))-ln(z)], z), `~`[op](1, expr__4), simplify)

true

(6)

NULL


 

Download senseless_results_of_int.mw

 

As you can see, the lengthy output of is nearly meaningless! (And if you want to simplify it, Maple will simply return: Error, (in simplify/recurse) indeterminate expression of the form 0/0.) So, how do I get the simplified results in Maple?
The integrals are: 

expr__1 := [Int(1/(z^2 + sqrt(z^4 + 4*a*z^2) + 4*a), z), Int(-1/(z^2 - sqrt(z^4 + 4*a*z^2) + 4*a), z)]: # (value(expr__1));
expr__4 := [Int((z^2 - 4*a*z + sqrt(-4*a*z^3 + z^4 - 8*a*z^2 + 4*z^3 - 4*a*z + 6*z^2 + 4*z + 1) + 2*z + 1)/(z*(-4*a*z + z^2 + 2*z + 1)), z), Int(-(z^2 - 4*a*z + 2*z + 1 - sqrt((-4*a*z + z^2 + 2*z + 1)*(z + 1)^2))/(z*(-4*a*z + z^2 + 2*z + 1)), z)]: # (value(expr__4)):

Note. By the way, Mma can solve the original ODEs directly and explicitly: 

In[1]:= DSolve[a*(y'[x]+1)^2+(y[x]-x)^2*y'[x]==0,y[x],x,IncludeSingularSolutions->Automatic]

                                   2                3                    2
                  a - x C[1] - C[1]             16 a  - 4 a x C[1] - C[1]
Out[1]= {{y[x] -> ------------------}, {y[x] -> --------------------------}}
                       x + C[1]                     4 a (4 a x + C[1])

In[2]:= DSolve[a*(x*y'[x]+y[x])^2-(y[x]+x)^2*y'[x]==0,y[x],x,IncludeSingularSolutions->Automatic]

                     2 a C[1]       2 a C[1]     2  2 a C[1]
                  a E         (-(a E        ) + a  E         + x)
Out[2]= {{y[x] -> -----------------------------------------------}, 
                                     2 a C[1]
                                  a E         - x
 
                2 a C[1]    2 a C[1]
               E         (-E         + 2 a x)
>    {y[x] -> --------------------------------}}
                    2 a C[1]              2
              2 a (E         - 2 a x + 2 a  x)

Unfortunately, Maple fails to do so.

(x^(3))^(1/3) doesn't simplify to x.  I am missing something.

The ODE is: 

eqn := y(x)*(2*x*diff(y(x), x) + y(x)*(diff(y(x), x)^2 - 1)) = -1: # How about another ODE 'lhs(eqn) = +1' ?

Maple can solve it, but I find that (to get all four solutions) I have to execute the dsolve command twice
 

restart;

eqn := y(x)*(2*x*(diff(y(x), x))+y(x)*((diff(y(x), x))^2-1)) = -1

dsolve(eqn, {y(x)}, 'parametric', 'singsol' = all)

y(x) = (c__1^2+2*c__1*x+1)^(1/2), y(x) = -(c__1^2+2*c__1*x+1)^(1/2)

(1)

dsolve(eqn, {y(x)}, 'singsol' = all)

y(x) = (-x^2+1)^(1/2), y(x) = -(-x^2+1)^(1/2), Int(-((_a^2+y(x)^2-1)^(1/2)*_a*y(x)^2-_a^2*y(x)^2-y(x)^4+2*_a^2+3*y(x)^2-2)/((y(x)^2+2*_a-2)*(-y(x)^2+2*_a+2)*(_a^2+y(x)^2-1)), _a = _b .. x)+Intat(-_f/(2*(_f^2+x^2-1)^(1/2)*x^2+(_f^2+x^2-1)^(1/2)*_f^2+2*x^3+2*x*_f^2-2*(_f^2+x^2-1)^(1/2)-2*x)-(Int(-(_a*_f^3/(_a^2+_f^2-1)^(1/2)+2*(_a^2+_f^2-1)^(1/2)*_a*_f-2*_a^2*_f-4*_f^3+6*_f)/((_f^2+2*_a-2)*(-_f^2+2*_a+2)*(_a^2+_f^2-1))+2*((_a^2+_f^2-1)^(1/2)*_a*_f^2-_a^2*_f^2-_f^4+2*_a^2+3*_f^2-2)*_f/((_f^2+2*_a-2)^2*(-_f^2+2*_a+2)*(_a^2+_f^2-1))-2*((_a^2+_f^2-1)^(1/2)*_a*_f^2-_a^2*_f^2-_f^4+2*_a^2+3*_f^2-2)*_f/((_f^2+2*_a-2)*(-_f^2+2*_a+2)^2*(_a^2+_f^2-1))+2*((_a^2+_f^2-1)^(1/2)*_a*_f^2-_a^2*_f^2-_f^4+2*_a^2+3*_f^2-2)*_f/((_f^2+2*_a-2)*(-_f^2+2*_a+2)*(_a^2+_f^2-1)^2), _a = _b .. x)), _f = y(x))+c__1 = 0, Int((_a^2*y(x)^2-2*_a^2+y(x)^4-3*y(x)^2+(_a^2+y(x)^2-1)^(1/2)*_a*y(x)^2+2)/((y(x)^2+2*_a-2)*(-y(x)^2+2*_a+2)*(_a^2+y(x)^2-1)), _a = _b .. x)+Intat(_f/(2*(_f^2+x^2-1)^(1/2)*x^2+(_f^2+x^2-1)^(1/2)*_f^2-2*x^3-2*x*_f^2-2*(_f^2+x^2-1)^(1/2)+2*x)-(Int((2*_a^2*_f+4*_f^3-6*_f+_a*_f^3/(_a^2+_f^2-1)^(1/2)+2*(_a^2+_f^2-1)^(1/2)*_a*_f)/((_f^2+2*_a-2)*(-_f^2+2*_a+2)*(_a^2+_f^2-1))-2*(_a^2*_f^2-2*_a^2+_f^4-3*_f^2+(_a^2+_f^2-1)^(1/2)*_a*_f^2+2)*_f/((_f^2+2*_a-2)^2*(-_f^2+2*_a+2)*(_a^2+_f^2-1))+2*(_a^2*_f^2-2*_a^2+_f^4-3*_f^2+(_a^2+_f^2-1)^(1/2)*_a*_f^2+2)*_f/((_f^2+2*_a-2)*(-_f^2+2*_a+2)^2*(_a^2+_f^2-1))-2*(_a^2*_f^2-2*_a^2+_f^4-3*_f^2+(_a^2+_f^2-1)^(1/2)*_a*_f^2+2)*_f/((_f^2+2*_a-2)*(-_f^2+2*_a+2)*(_a^2+_f^2-1)^2), _a = _b .. x)), _f = y(x))+c__1 = 0

(2)

NULL


 

Download dsolve_twice.mw

However, in MATLAB®, the complete solutions can be found just in one go

>> dsolve('y*(2*x*Dy + y*(Dy^2 - 1)) = -1', 'x') % require the Symbolic Math Toolbox™
ans =
                         1
                        -1
 -(-(x - 1)*(x + 1))^(1/2)
  (-(x - 1)*(x + 1))^(1/2)
 (C1^2 + 2*x*C1 + 1)^(1/2)
-(C1^2 + 2*x*C1 + 1)^(1/2)

Does anyone know why?

Ideally, I would like to find all roots of this RootOf expression for a given interval.

I tried defining a function from the argument of the RootOf expression and using fsolve to find solutions, but could not get all of them.

What I managed to do skips the interval, is not really elegant and raises additional questions.
I would be grateful for any hints and improvements.

RootOf(_Z*cos(_Z)-sqrt(sin(_Z)^2))

RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))

(1)

allvalues(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2)))

Error, (in RootOf/sort1) cannot numerically evaluate the argument

 

RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), index = i)

RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2), index = i)

(2)

evalf(subs(i = 3, RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2), index = i)))

RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2), index = 3)

(3)

evalf(RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 3))

-4.493409458

(4)

rt := ''RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), i)''

'RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), i)'

(5)

subs(i = 3, rt)

RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 3)

(6)

evalf(RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 3))

-4.493409458

(7)

seq(subs(i = k, rt), k = 1 .. 5)

RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 1), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 2), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 3), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 4), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 5)

(8)

evalf(%)

0., 0., -4.493409458, 4.913180439, 4.913180439

(9)

{-4.493409458, 0., 4.913180439}[]

-4.493409458, 0., 4.913180439

(10)

seq(evalf(subs(i = k, rt)), k = 1 .. 5)

Error, (in evalf/RootOf) numeric exception: division by zero

 

evalf(seq(subs(i = k, rt), k = 1 .. 5))

Error, (in evalf/RootOf) numeric exception: division by zero

 

NULL

seq(subs(i = k, rt), k = -5 .. 5)

RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), -5), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), -4), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), -3), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), -2), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), -1), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 0), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 1), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 2), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 3), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 4), RootOf(op(RootOf(_Z*cos(_Z)-(sin(_Z)^2)^(1/2))), 5)

(11)

evalf(%)

-4.493409458, -4.493409458, -2.028757845, -2.028757838, 4.913180439, 0., 0., 0., -4.493409458, 4.913180439, 4.913180439

(12)

NULL

NULL

Download RootOf_a_periodic_function.mw

Is there any setting that controls the extent of a plot?

Left hand plot has defined extent of the plot, while the plot on the right hand side has not. When panning the graphics on the right side the plot is clipped.

Any idea how to make Maple to use the whole extent of the plot component as a boundary?

Download plotpoint2.mw

First 83 84 85 86 87 88 89 Last Page 85 of 2308