I need to change the following continued fractions into rational numbers using the cfrac function:

1)  [1 , 2]               2) [1, 1, 2]                3) [1, 1, ,1, 2]

Thanks,

Sara

I know there is some code to compute the symmetry of pde, but i want to compute the symmetry of some differential expression, as f(x,t,ux,ut).

The difference is that when solve the infinitesimal, we don't have the identity f(x,t,ux,ut)=0, we just take each variable in the jet space as independent.

Is there any package to do that?

help and thanks!

Greetings,

I have been attempting to plot a couple of functions and I keep getting the same error.  I will copying and paste what I have entered and the error messages

with(plots);

with(VectorCalculus);

sphereplot*(rho = theta, theta = 0 .. 4*Pi, phi = 0 .. (1/2)*Pi);

Error, (in simpl/reloprod) invalid terms in product: (rho = theta)*(theta = 0 .. 4*Pi)

with(plots);
with(VectorCalculus);
fieldplot*(<-x, y, cos z>, x = 0 .. 5, y = 0 .. 5, z = 0 .. 2*Pi);
Error, (in simpl/reloprod) invalid terms in product: (x = 0 .. 5)*(y = 0 .. 5)

Any assistance would be most appreciated.  Thanks in advance.

Hello,

how can I use procedures saved in ".m" files in my own repository?

I've created

LibraryTools:-ShowContents("C:\\Users\\Stefan\\maple\\toolbox\\myfactor\\lib");
[["Auxiliary.m", [2016, 10, 10, 15, 53, 20], 1024, 3333]]

and registered it in the maple.ini
libname;
     "C:\Program Files\Maple 18\Maple 18\lib",
       "C:\Users\Stefan\maple\toolbox\Iterator\lib", ".",
       "C:\Users\Stefan\maple\toolbox\myfactor\lib"

but Maple don't know the procedure squarefree() from the repository. What is wrong?

Thanks

Hi all,

I am trying to plot in semilog scale a function involving products of exponential integrals and complex exponentials. For small and moderate values of the argument, the plot is well shown. However, for larger values the plot shows strong fluctuations. I was wondering how one can deal with such a problem. Any help is highly appreciated.

Please refer to the attached script for the functio of interest.

Thanks
F

Question.mw

I understand that if I want to use a scientific constant in a Maple worksheet or document, I first have to declare it. Here is an example for the speed of light:

However if this worksheet is re-executed (e.g. using !!!) this happens:

Obviously there is a kind of recursion for c which is listed in the Variable Palette as a name for a variable but which is also the scientific constant. The problem can be solved e.g. by associating another name to the constant and unassigning the old definition in the Variable Palette:

This works but it makes the document harder to read since everybody is used to lowercase c for the speed of light. The same is of course true for other scientific constants.

Is there a way to use scientific constants in their "usual" notation in Maple, including using the unit of the constant?

Dear all!.       

I have an expression with ramdom variables. Can i use NLP to optimize it?. In this case, does Maple take the histogram of theses ramdom variables?. Thank!!!

What is the difference between the Expression Palette and the Layout Palette? First I thought that the Layout Palette is supporting the entering of names with e.g. sub- and superscripts, but I learned (the hard way) that the templates of the Layout Palette are "mathematically active".

What is the meaning of the colours green and violett in both palettes?

 Dear All ! 

I really need to solve this problem as soon as possible, As you know the downside equation is not exact, but I can not find its integration factor, blease help me !

                                                                 ∫{ ( ωx + σy ) d x + (ωy −σx)dy}=0

 Regards ,

Hi, so I am a newbie with simulation on maplesim. I've made a simulation of the lower part of a humanoid robot and I attempt to input angles to the joints from an excel file.

However, after adding my file .xls to the data set and linking it to a time lookup table, I get this error:

Havent installed any add-ins, I am working with excel 2007, Maplesim 6.4 and Maple 18.

Any insight on the problem would be of great help.

Thank you

Hi everybody, 

I generate random variables of the form X__E__n, n=1..N where N is assigned to some integer value (the double index is not a crucial point for the issue already appears with a single level of index).
Next I want to express a likelihood based upon these random variables in two different ways : a very symbolic one which uses the function "Product" an a second one that uses "product".

My problem is that I can't obtain the desired expressions.
As "Product" returns me a completely satisfactory expression, "product" (as I use it) refuses to return me the desired expression.

Without doubt this is due to a poor knowledge of the detailed way Maple proceeds to evaluate indices in seq, sum or product.
All explanation will be largely appreciated.

Thanks in advance

Download WhatHappensHere.mw

I'm trying to make 3D-animation and the procedure takes an enormous amount of time to make frames for it. It took over 10 minutes to create animation consisting of 50 frames (over t=0..1000). If I increase the number of frames or range of time, it just wouldn't reach the end of calculation (keeps evaluating till I lose patience or my faith in Maple)

Is there a way to somehow precompute data for animation (maybe store it somewhere) so that it could gather it from there and the animation construction itself would be faster?

Here's the code: ClassicalTrajectoriesH2X_harm.mw (animation is in the end of the file)

P.S. I have Asus X555LJ; Intel Core i5-5200U, 2.2GHz; 12Gb RAM, for the last few months it presented itself as a rather fast piece of machinery. 

Hello everybody,

I need help to increase the speed of two procedures.

My first (recursive) procedure returns a list of all k - element subsets of a given set {1, ... , n}:

nkSubSetList := proc(n :: nonnegint, k :: nonnegint)
        option remember;          # boosts with factor 5
        local X, prev, nextOne;

        # Local function that augments a nonnegint to a given set, e.g.
        local augment := (SetA, intm) -> {op(SetA), intm};

        # Case k > n allowed & needed for recursive call
        if n < 0 or k < 0 or k > n then
            X := {};
            return X
        elif n = 0 or k = 0 then
            X :=  [{}];
            return X
        else
            prev := nkSubSetList(n - 1, k);
            nextOne := map(augment, nkSubSetList(n - 1, k - 1), n);
            return [op(prev), op(nextOne)]
        end if;
end proc:

# Starting the Profiler on my system the call nkSubSetList(20, 10) needs 0,5 sec.
The other algorithm returns a list of all k-element subsets of a given set T:

TkSubSetList := proc(T :: set(posint), k :: nonnegint) :: list(set(posint));
    local n := nops(T);                                 # Cardinality of index set T
    local callSet := A -> map(i -> T[i], A);     # Call the i-th entry of the given set T

    # Input ASSERTS
    ASSERT(k <= n, "k is bigger than the cardinality of index set T");

    return map(callSet, nkSubSetList(n, k));
end proc:

# Starting the Profiler on my system the call TkSubSetList({1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20}, 10) needs 1,8 seconds.

Who have an idea how to increase the speed of these algorithms?

Can somebody of Maple users execute the following command

restart; pdsolve({diff(w(x, y, z), x)+diff(w(x, y, z), y, y)+2*(diff(v(x, y, z), x)-(diff(u(x, y, z), y))-2*w(x, y, z)) = diff(w(x, y, z), z, z), 3*(diff(u(x, y, z), x, x))+2*(diff(u(x, y, z), y, y))+2*(diff(v(x, y, z), x, y))+2*(diff(w(x, y, z), y)) = diff(u(x, y, z), z, z), 3*(diff(v(x, y, z), y, y))+2*(diff(v(x, y, z), x, x))+2*(diff(u(x, y, z), x, y))-2*(diff(w(x, y, z), x)) = diff(v(x, y, z), z, z)}, {u(x, y, z), v(x, y, z), w(x, y, z)})

in Maple 2016.1.1 on a powerful comp and report the obtained result as an answer to the question?
 That would be very kind of her/him. Thanks in advance. 

Download pdsolve.mw

I have a system of pde as follow,

PDE := [(x*y+1)*(diff(f(x, y), y, y, y))+(x+(3/4)*f(x, y))*(diff(f(x, y), y, y))-(1/2)*(diff(f(x, y), y))^2+T(x, y) = (1/4)*x*(diff(f(x, y), y))*(diff(diff(f(x, y), y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(f(x, y), y, y)), (x*y+1)*(diff(T(x, y), y, y))/(.733)+(x/(.733)+(3/4)*f(x, y))*(diff(T(x, y), y)) = (1/4)*x*(diff(f(x, y), y))*(diff(T(x, y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(T(x, y), y))]

sys_ode := diff(g(y), y, y, y)+(3/4)*g(y)*(diff(g(y), y, y))-(1/2)*(diff(g(y), y))^2+h(y) = 0, (diff(h(y), y, y))/(.733)+(3/4)*g(y)*(diff(h(y), y)) = 0

ics := g(0) = 0, h(0) = 1, (D(g))(10) = 0, g(10) = 0, h(10) = 0

sol2 := dsolve([sys_ode, ics], numeric)

BC := {T(0, y) = h(y), T(x, 0) = 1, T(x, 10) = 0, f(0, y) = g(y), f(x, 0) = 0, f(x, 10) = 0, (D[2](f))(x, 0) = 0}

pds := pdsolve(PDE, BC, numeric)

module() ... end module

pds:-plot(T, y = 0 .. 10, x = 0);

Error, (in pdsolve/numeric/plot) unable to compute solution for x<HFloat(0.0):
solution becomes undefined, problem may be ill posed or method may be ill suited to solution

When I try to use the solution of the Ode as the boundary condition for PDE, by subbing g(y) and h(y) into BC. The plot returns me the error. Anyone knows the reason behind this and how to solve? Any help would be really greatly appreciated. Thanks