Maple 2017 Questions and Posts

These are Posts and Questions associated with the product, Maple 2017

Let us consider

sol1 := dsolve({diff(y(x), x) = solve((1/2)*(diff(y(x), x))^2 = (1-ln(y(x)^2))*y(x)^2, diff(y(x), x))[1], 
y(0) = 1}, numeric);
sol1 := proc(x_rkf45) ... end proc

The problem under consideration has the symbolic solution:

sol2 := dsolve({diff(y(x), x) = solve((1/2)*(diff(y(x), x))^2 = (1-ln(y(x)^2))*y(x)^2, diff(y(x), x))[1], 
y(0) = 1});

sol2 := y(x) = exp(x*sqrt(2)-x^2)

Let us compare the plots of sol1 and sol2 (which should coincide):

A := plots:-odeplot(sol1, x = 0 .. 1, color = navy, style = point):
B := plot(rhs(sol2), x = 0 .. 1, color = red):
plots:-display([A, B]);

The plots differ after approximately 0.707. Bug_in_dsolve_numeric.mw

Edit. The title and one of the tags.

evalf(Int(x*(1-2*x^(3/10))^(10./3),x=0..1));  # Crashes Maple

Note that:

int(x*(1-2*x^(3/10))^(10./3),x=0..1);
int(x*(1-2*x^(3/10))^(10/3),x=0..1);

are OK.

(Windows 7, Maple 2017.3, 64 bit)

Application to generate vector exercises using very easy to use patterns. These exercises are for sum of vectors in calculations of magnitude, direction, graph and projections. Basically to evaluate our students on the board and thus not lose the results. For university professors and engineering students.

Generator_of_exercises_with_vectors.mw

Lenin Araujo Castillo

Ambassador of Maple

This is my first time in this forum, so I hope I use the correct conventions. If not please notice me. 

  1. When a light body orbits a heavy body under the influence of gravity (e.g. a planet around the Sun), Newton’s laws show that the orbit is restricted to a two-dimensional plane and is given by the differential equation

    d2/d(φ)2(1/r(φ)) + 1/r(φ) = GM/h2

    Here, (r, φ) is the path of the light body in polar coordinates, M is the mass of the heavy body, G is the gravitational constant, and h is a constant related to the angular velocity of the light body (h = r2φ ̇). The heavy body can be considered to be approximately stationary and located at the origin.

    Use Maple to solve this differential equation numerically, taking M = 1, G = 1, h = 1 with initial conditions

    r(0)=2/3, r′(0)=dr/dφ (0) = 0

    Using polar coordinates, create a plot of the orbit (r(φ), φ) for
    0 ≤ φ < 2π. You should observe a perfect ellipse.

  2. Since I am not a frequent maple user, I hope somebody can help me here

Hi.... I'd want to numerically solve a system of  n  polynomial equations of degree 2 with respect to n unknowns. It seems to me that there is a software in Maple (which deals only with the equations of degree 2) that solves this problem. I cannot find it anymore. Do you know it ? Thanks in advance.

For a numerical integration in Maple, one only adds evalf(...) to the integration commend. Now my question is that is it also the case for the multi-integration? Also consider that the multi-integration of the Maple doen't give any answer for some examples. So I think putting evalf around them, doesn't help. How can someone compute a numerical approximation of a multi-integration by Maple? Is there any prepared commend or package for it?

I have a code in Maple that solves a system of polynomial equations (18 equations, 18 variables) using the command 'RootFinding[Isolate]' (it takes to Maple about 40 minutes to find all the solutions to the equations). Is it possible to see the implementation of the Isolate command? if not, is there any other open source code that implements the equivalent to Isolate command and also able to solve such equations? (I'm talking about the  numeric solutions. not symbolic).

I already tried to use the command 'solve_poly_system' in python (sympy package) which is open source. but the function didn't return a solution after a very long time (and I used a system of 9 equations and 9 variables which were much simpler than the equations Maple was able to solve using Isolate command).

I'm attaching the equations Maple solved using Isolate command. and also the equivalent in python with the simpler equations (that it wasn't able to solve).

solving_18_polynomial_equations.mw

trying_to_solve_polynomial_equations_python.txt

I'd appreciate any idea of how to use open source code that produces the same results as the code in Maple (if it's possible).

Thanks

 

Hey everyone
 

I am trying to create a procedure, which includes a part, where a function is plotted. I want to be able to declare the range of the plot within the procedure parameters as x=range - example when executing: procedureName(function,x=range). I have isolated a small part of the procedure and the different approaches I have tried:

For simplicity's sake let's call the procedure test.

test:=proc(function,{x::range:=0..1})

#I would then like to be able to plot the following:

plot(function,'x'=x);

end proc:
 

However, the single quotes do not prevent x from being evaluated within the procedure, and there interprets the second argument of 'plot' as 0..1=0..1. Another approach I tried was the following:

 

test:=proc(function,{x::range:=0..1})

local xRange:=x;

unassign('x');

plot(function,x=xRange);

end proc:

 

However, I run into the same problem as above - the elements within the single quotes are still evaluted within the procedure. A workaround would be to simply rename the range parameter, for example:

 

test:=proc(function,{xInr::range:=0..1})

plot(function,x=xInr);

end proc:

 

I would really like to be able use x as the parameter, though.

 

Thanks in advance!

Here is my code:

with(DEtools): with(plots): with(linearAlgebra):

DE1 := diff(y(x), x) = d^3*y(x)/dx^3+3*d^2*y(x)/dx^2+4*d*y(x)/dx+12*y(x);

DEplot(DE1, [y(x), x = 0 .. 5][[y(0) = -3, (D(y))(0) = 0, ((D^(2)(y))(0) = 0, ((D^(3)(y))(0)]]);

I get the Error message: Error, (in DEtools/DEplot) called with too few arguments

Any help would be greatly appreciated, thanks.

My question is within the worksheet.

assignation.mw

In the context of homework, there are some calculations that I do through Maple, and others that I do mostly by hand, and just show my work with Maple's engine that is great at rendering special mathematical formats (like matrices, for example).

I would like to show something like that:

 

It's pretty easy to make the matrices, right now my line looks like this:

 

But I would like to refrain for executing the "+" operator and the multiplication by x6 and x7, to have the output like in the first picture. If I disable all execution and take the statement as plain text, then obviously I will not get the correctly rendered matrices.

Is there a way, like an escape character, to prevent specific operations from being executed in one line in particular?

how I can  select 2D function (sigma) that these boundary conditions are satisfied?

thanks

sigma(1, z) = 0, sigma(r, 0) = 1, diff(sigma(0, z), r) = 0, diff(sigma(r, 1), z) = 0

Hi,

I have an inverse I solved for with RootOf, which I transfrom and then try to numerically integrate the transformation. I get the following error: Error, (in convert/radical) too many levels of recursion

Integrating the inverse itself works, it is after applying the transformation 1/(1+x2) that I get the issue. Eventually I want to integrate something less simple, but here the problem already starts. 

Thanks for the help! Its for science! 
 

g1_inveq := proc (t, x, alp) options operator, arrow; tan((x-t)/(alp*(1+t^2)))-t end proc

proc (t, x, alp) options operator, arrow; tan((x-t)/(alp*(1+t^2)))-t end proc

(1)

g1_inv := proc (v, alp) options operator, arrow; evalf(RootOf(g1_inveq(t, v, alp) = 0, t, 1)) end proc

proc (v, alp) options operator, arrow; evalf(RootOf(g1_inveq(t, v, alp) = 0, t, 1)) end proc

(2)

evalf(Int(1/(1+g1_inv(v1, .5)^2), v1 = 0 .. 1))

Error, (in convert/radical) too many levels of recursion

 

``


 

Download recursionerror.mw

 

Hello people in mapleprimes,

I have a question.
I want to modify

a := rho^(epsilon-1)*a__01^(-k)*(N__E2/Omega+N__E1*mu)*(K+1)/mu;
to
b:=rho^(epsilon-1)*a__01^(-k)*(N__E2/(Omega*mu)+N__E1)*(K+1);

But, I can't do that modification.

The key term is
c:=N__E2/Omega+N__E1*mu;

I know how to use applyop, but even if I use it, I couldn't obtain b.

Please teach me about this.

Nov_22_2017_1_mapleprimes.mw

 

 

Hello people in mapleprimes,

I want to know some way, if there is, to do the modification of the following, that is

The original expression: (beta*L__1/P)^beta*((1-beta)*L__1)^(1-beta)
The expression after modification: beta^beta*L__1*(1/P)^beta*(1-beta)^(1-beta)

I did this with freezing (1-beta) in (1-beta)*L__1.
But, this is a little long code. I will be glad if you tell me some more simple way to do this.

Thanks in advance.

nov_21_2017_mapleprimes.mw

 

 

First 23 24 25 26 27 28 29 Last Page 25 of 40