140 Reputation

7 Badges

9 years, 164 days

MaplePrimes Activity

These are replies submitted by maple2015

@Mac Dude 

Thank you for your hints

In the case of more than two variables, is it possible to use  "ChiSquareGoodnessOfFitTest" ?


Thank you

@Carl Love 

I hope you have a good sleep.

I attach a PDF to explain the proposed algorithm which is used in the "" worksheet. As I mentioned before, the convergence is decreased by moving from vicinity of the initial point.

I will be grateful, if I can know more about digging into the Maple code (method= classical[adambash]) ).


@Carl Love 

The mentioned formula is used for first order ODE. I modify it for 2nd order ODE.

I think Maple uses some controller steps or something like that to guarantee the convergence of the numerical solution.

I am very appreciate for your contribution.

Your hints make it possible to dig through the algorithm.

@Carl Love 

Thank you for your hints

The plot below shows difference between exact and approximated solutions based on my algorithm.

Plot (odeplot command) of the numerical solution for method= classical[adambash] shows a complete agreement between exact and approximated solutions.

If it is possible, can I find the correct algorithm which is used by Maple?


{eval(phi, y = -sqrt(-z^2+1)), eval(phi, y = sqrt(-z^2+1))}

Two last equations (above mentioned Eqs.) make balance between 39 equations and 41 unknowns.

The function phi is obtained in terms of two variables y and z, therefore the result of double integral is a number.

I just want to solve double integral in less possible time.

Is there a way to calculate it in less possible time?


Thank you


Your code rectifies the deficiencies of previous code


@Carl Love 


Thanks for your valuable hints


Thank you for your hints

I expected that P should be approached to zero. Now I can be sure that my assumptions are true!

If it is possible, please help me to find the amount of Omega in attached file

omega^2-.3846153846*(Int(Int((5/2)*(1182993435*omega^2-413891981)*(cos((1/650240945)*sqrt(1099314545041681865*omega^2-384615384486573943)*r*sin(t))-cos((1/65024094500000)*sqrt(1099314545041681865*omega^2-384615384486573943)*sqrt(-10000000000*r^2*cos(t)^2+11283791670)))*r/(cos((1/65024094500000)*sqrt(1099314545041681865*omega^2-384615384486573943)*sqrt(-10000000000*r^2*cos(t)^2+11283791670))*(8453132285*omega^2-2957483587)), t = 0 .. 6.283185308), r = 0 .. 1.062251932)) = 0





At y=0 the function 'psi' in attached file and its second derivative are vanished. Also by taking G as a constant amount, 'G,y' is zero. Therefore the problem has solution.

We haven't the integration constant. By assuming this constant equal to zero, I try to solve the problem (it is a guess!)

My questions are

1- The MAPLE can not solve the last equations in by fsolve. What is the reason? The problem hasn't answer or some codes like the Newton's iterative method must be used instead of fsolve?

2- Is it possible to minimize 'phi' by using two coupled nonlinear equations in If yes then which commands can be used?

3- Please suggest a method to solve the problem, if it is possible.

Thanks a lot



Please find the attachment

Yours sincerely





The ODE exp(B*x)*w''+A*w=0 is stability equation of a rod and we can not change it to

exp(B*x)*(diff(w(x),x,x)+A*w(x))=0 which is equivalent to diff(w(x),x,x)+A*w(x)=0 .

When B approaches to zero we have a rod with uniform mass distribution


Thank you anyway




The PDE is arisen from stress analysis of a solid


I have another question:

Consider the ODE  (w''+A*w=0)

Assuming (A>0) yields w=_C1*sin(sqrt(A*x))+_C2*cos(sqrt(A*x))

Now consider the ODE  (exp(B*x)*w''+A*w=0)

Why when B approaches to zero the answer of second ODE (Which is in terms of Bessel functions)  is not


and approaches to zero?

limit(rhs(dsolve(exp(B*x)*(diff(w(x), x, x))+A*w(x) = 0)), B = 0);

@Joe Riel 


There is no any backup unfortunately

Also I do the instruction which is described in Online Help:

But it dosent work again!


Thank you any way

1 2 3 4 5 6 Page 3 of 6