in this link   https://maplesoft.com/applications/Preview.aspx?id=123874  there is paper/application 

which is supposed to have Maple package called with(PSsolver);

But when downloading the application (by clicking on the download link) and opening it, it only has one file 

              an_extension_of_prelle-singer_method_and_maple_implementatio.mws

Where will the package (source code) itself PSolver be?  I can't see any other link in the Maple page above to download additional files which might include the package itself.

The worksheet above says this

Before actualy introducing the ideas concerning our program, let us first load it to the Maple section. (Please, put the library file into directory `c:/DirectoryWhereIputtheLibrary`).

I googled and not able to find it also. This was written almost 24 years ago, so trying to contact the author now might not be good idea.

Should not the Maple application be all self contained with all the files needed to run it?

Any one has copy of this package by any chance?

the most paper use another function to get the result and then do substitute i try to get by the way of them but i fail so how about if we can get the results in direct function there is any way for finding thus parameter after substitution of our function in ode?

>  >  >  >  >  >  (1) >  >  Warning, A new binding for the name `gamma` has been created. The global instance of this name is still accessible using the :- prefix, :-`gamma`.  See ?protect for details.   >  (2) >  (3) >  (4) >  (5) >  Error, (in unknown) incorrect use of identity(<expr>,<name>)

Download Find_params.mw

Is there any function that returns a boolean value that tells me if a point is within a polygon or not?

I tried to contact an author using the menu item more->contact autor.

I got the following back from postmaster@maplesoft.com

The IP-adress 199.71.183.16 (domain maplesoft.com) does not match the IP-adress of my mailprovider.

The authentification of the mailprovider works when I send directly from this mailprovider to gmail accounts.

Is this a known issue when sending via maplesoft.com or only temporary?

Is gmail too restrictive?

Problem Q15 in the book Parabolic Problems by David Angell and Thomas Britz describes a large circle (LC) and several smaller circles (SCs) which are each tangent to its neighbour SC(s), and externally to LC. All circles are tangent to the x axis and above it.

Section one of this worksheet displays the LC and six of the SCs based on the book's formula for the diameter of the latter in terms of the diameter of the LC and the largest SC, which is determined by the user.

Section two finds and displays that all of the displayed SCs' centers lie on the diameter of a circle closely related to the LC and larger than it.

Can this be proved to be the case for any sizes of the LC and SCs in the same formation as that displayed?

Parabola_Problems_Q15.mw

I need assistance with building the homotopy analysis method to solve the system of odes. here is the attempt to do it. I'm still new to maple restart; # Declare functions for the system PDEtools[declare]([f(x), g(x)], prime = x): # Order of expansion N := 4; # Define series for each function f1(x) := sum(p^i*f[i](x), i = 0..N): f2(x) := sum(p^i*g[i](x), i = 0..N): # Define the system of ODEs using Homotopy HPMEq1 := (1-p)*(diff(f1(x), x$3)) + p*(diff(f1(x), x) + f2(x)): HPMEq2 := (1-p)*(diff(f2(x), x$3)) + p*(diff(f2(x), x) - f1(x)): # Extract coefficients for the system for i from 0 to N do equl[1][i] := coeff(HPMEq1, p, i) = 0: equl[2][i] := coeff(HPMEq2, p, i) = 0: end do: # Define boundary conditions for the system cond= 0, D(f1[0])(5) = 1]: cond = 0, D(f2[0])(0) 1]: for j from 1 to N do cond = 0, D(f1[j])(0) = 0, D(f1[j])(5) , D(f2[j])(0) = 0, D(f2[j])(5) = 0]: end do: # Solve the system iteratively for each order for k from 0 to N do dsolve([equl[1][k], cond[1][k]], f1[k](x)): dsolve([equl[2][k], cond[2][k]], f2[k](x)): end do:

i try to get same result by substituation but i don't know what is mistake after i take second derivative is wronge i don't know how get same result as in paper did can anyone help  to calculate this part is not hard but is complicated ,How calculated second derivative and put in our ode to get the parameters?

>  >  >  >  >  >  (1) >  >  >  (2) >  >  (3) >  (4) >  (5) >  (6) >  >  (7) >  (8)

Download complex-issue.mw

How do I access the Maple Fundamentals guide video? Every time I click on the video, it tells me I must log in to view the video. But there is no where on that page to log in. When I log in to Maple Primes, I can't find the video.

I think the worksheet below is enough to define the problem I'm facing.

(question updated by adding a simple 2D case ; the expected value is about 0.67684)

>  restart >  # How can I define CharFunc so that int(CharFunc, x=...) gives me the expected value? # # Motivation: # Let Dom some 2D domain and Env another 2D domain which contains Dom. # The interior of Dom is defined by a sequence L(x, y) of inequalities. # The area of Dom can be expressed by the integral over Env of the characteristic function # CharFunc(x, y) of Dom which returns 1 if (x, y) belongs to Dom and 0 otherwise. # # Here is a simple 1D example CharFunc := proc(x)   if x::numeric then     piecewise(is(And(x >= 0, x <= 1)), 1, 0)   else     'procname( _passed )'   end if: end proc: Env := -1..2: plot(CharFunc(x), x=Env, thickness=5); >  # Expected value = 1 int(CharFunc(x), x=Env); (1) >  # Expected value = 1/2 int(CharFunc(x), x=1/4..3/4); (2) Works with floats >  # Expected value = 0.5 int(CharFunc(x), x=0.25..0.75); (3) >  # Expected value = 1 int(CharFunc(x), x=-1.0..2.0); (4) 2D example >  Env := x=0.8..3, y=0..1.3; (5) >  Dom := proc (x, y) options operator, arrow; And(y <= 1/(1+sinh(2*x)*ln(x)^2), .8 <= x, x <= 3, 0 <= y, y <= 1.1) end proc (6) >  CharFunc := proc(x, y)   piecewise(is(Dom(x, y)), 1, 0) end proc: >  plot3d(   [0, 'CharFunc'(x, y)], Env   , grid=[40, 40]   , style=surface, color=[gray, blue], transparency=[0, 0]   , title=typeset('CharFunc'(x, y)) ); >  int(CharFunc(x, y), Env) (7)

Download Integration_issue.mw

Thanks in advance

Does Maple provide any command/examples that can help verfying whether an analytic function or a time series (a sequence of numbers) is periodic?

I had a look at signal processing. In principle auto-correlation could be used for that purpose but I am interested in a command that outputs a period for a given allowable error or a graphical display that can help in an assessment.

Analytic convolution and analysis of the convolution terms for periodic conditions would be great for analytic functions.

I have no time (and limited skills) to program such things myself.

Hello everyone, first time poster here.

I am trying to calculate xDLVO theory for a surface, which is simply 3 equations with 3 unknown. I have values for everything except x, y and z and Maple gives intermediate results for all the equations when defining them. The cos(x) is of a contact angle, and the problem is that Maple will not calculate x, y and z when all three cos(x) have decimal numbers. when changing just one of them to a whole number, such as changing cos(66.969) to cos(67), then it solve the unknowns. This is the same for the other two numbers, there just always have to be a atleast one whole number in one of the cos(x). I have tried to just use the results of the cos(x), so I enter decimal numbers instead of the entire (1+cos(x) bracket, and it still will not work with decimal numbers.

I am not the brightest at math or Maple, so maybe it is a dumb mistake from my part :)

 

eqw := Ytotw*(1 + cos(66.696)) = 2*(sqrt(Ylww*x) + sqrt(y*Ypw) + sqrt(Ymw*z));

           eqw := 17.99142992 = 2*sqrt(26)*sqrt(x) + 10.09950494*sqrt(y) + 10.09950494*sqrt(z)

eqf := Ytotf*(1 + cos(42.497)) = 2*(sqrt(Ylwf*x) + sqrt(y*Ypf) + sqrt(Ymf*z));

           eqf := 62.95291273 = 2*sqrt(39)*sqrt(x) + 3.019933774*sqrt(y) + 12.58570618*sqrt(z)

eqd := Ytotd*(1 + cos(30.405)) = 2*(sqrt(Ylwd*x) + sqrt(y*Ypd) + sqrt(Ymd*z));

          eqd := 77.77865435 = 14.25482374*sqrt(x)

solve([eqw, eqf, eqd], [x, y, z]);
                               []

 

Thank you for your time!

Best regards

Jakob

Hello, I need help to find the real and Imaginary components of following complex model using maple, here q=u+iv,

The maple worksheet shows an incorrect evaluation of the integral in (1) which is a standard integral representation of a Bessel function.  Equations (2)-(5) along with the graph show the incorrectness of the evaluation.  What is going on?

Bessel.mw

This worksheet defines two physics problems and fails in the attempt to solve the first one.

How can these problems be solved?

Rolling_circle.mw

Hello, can somebody give some suggestions to speed up this summations? I am summing up to 10 and it take ages (about 10min).

Digits := 10;
n1 := 122;
x := 54;
n2 := 65;
y := 14;
alpha1 := 1.1;
beta1 := 1.1;
alpha2 := 1.1;
beta2 := 1.1;

Sigma2 := (j, l, psi0) -> Sum(Sum(Sum(psi0^(-y - alpha2 - l - alpha2*h - v)*pochhammer(beta2 - 1, h)*pochhammer(beta1 - 1, u)*pochhammer(y + alpha2 + l + alpha2*h, v)*2^(-y - alpha2 - l - alpha2*h - x - alpha1 - j - alpha1*u - v)*(psi0 - 1)^v/(h!*(y + alpha2 + l + alpha2*h)*u!*v!*(y + alpha2 + l + alpha2*h + x + alpha1 + j + alpha1*u + v)), u = 0 .. upto), h = 0 .. upto), v = 0 .. upto);

F2 := psi0 -> Sum(binomial(n1 - x, j)*Sum(binomial(n2 - y, l)*(-1)^(j + l)*Sigma2(j, l, psi0), l = 0 .. n2 - y), j = 0 .. n1 - x);

upto:=10;

F2(5.620);

Thanks in advance.