Hi,

I am wondering how we can edit a diagram plotted in maple. for example in editdiagram.mw I want name tips of diagram as "a", "b",.... and write them besides the tips on diagram.

Another question is how we can make color different parts of a single diagram by various colors? for example in above diagram we want make blue vertical line, make green horizontal line and make black the curve line?

Thanks a lot

In Maple 2017.2:

solve(x^3+x*a+2 > 0, [a, x]);
           [[a = -3, x < 1, -2 < x], [a = -3, 1 < x]]

solve(x^4+x*a+1 > 0, [a, x]);
Error, (in rootbound) 1st argument must be a polynomial with numeric coefficients

The first one clearly doesn't give the complete set of solutions. It can be made to work by adding conditions on a such as a>0.

I'm not sure if symbolic operations take index=real[i] into account or just silently ignore it:

limit(RootOf(_Z^3+a*_Z+2, index = real[1]), a = -infinity);
                               0

In fact the root grows as sqrt at -infinity.

Hello,

When starting Maple, i get an error message like:

All Drivers are up to date, Im running a system Win10-64Bit on a notebook.

Has anybody had the same problem?

Thanks in Advance

hi every one,i want to create some kind of numbering for index, how can i generate them using "seq" command or "for loop" ?
 

Download problem.mw

solve(z^(1+I) = 1, z, allsolutions = true);
                      exp((1 + I) Pi _Z1)
getassumptions(_Z1);
                         {_Z1::integer}

Since this is a single-valued power function, there is only a finite number of solutions. evalc correctly gives exp(-2*Pi) for _Z1=-1.

evalf doesn't help here regardless of the level of precision, I think because it always generates a non-zero imaginary part for exp(-(1+I)*Pi):

seq(print(evalf(evalf[d]((exp((-1-I)*Pi))^(1+I)))), d = 10 .. 3010, 300);
                                             -12  
              0.001867442732 - 1.361179007 10    I
                                            -312  
             0.001867442732 - 1.674479874 10     I
                                           -610  
               1.000000000 + 2.386571217 10     I
                                           -910  
               1.000000000 + 8.502509375 10     I
                                            -1212  
             0.001867442732 - 1.646483173 10      I
                                            -1514  
             0.001867442732 - 4.556560265 10      I
                                            -1812  
             0.001867442732 - 1.287611101 10      I
                                            -2112  
             0.001867442732 - 1.072784224 10      I
                                          -2410  
              1.000000000 + 8.162729354 10      I
                                            -2713  
             0.001867442732 - 7.375390371 10      I
                                          -3010  
              1.000000000 + 1.988371005 10      I

Unrelated, but it would be nice to have a simple way to display lists/matrices with specified width and alignment.

In my code I'm trying to catch the case when a certain limit doesn't exist:

limit(exp(-I*a*z)/a, a = infinity) assuming Im(z) > 0;
                               0

Doesn't work. Strangely, this works:

limit(exp(-I*a*z)/a, a = infinity) assuming Im(z) > 0, z::complex;
                           undefined

But this again doesn't:

limit(cos(a*z)/a, a = infinity) assuming Im(z) > 0, z::complex;
                               0

(because the limit of cos(a*z) returns -1..1).

So then I decide to check the exponents directly:

is(a*Im(z) > 0) assuming a > 0, Im(z) > 0;
                             false

Doesn't work, even though it works for arbitrary f:

is(a*f(z) > 0) assuming a > 0, f(z) > 0;
                              true

I found a construct that works:

signum(Im(a*z)) assuming a > 0, Im(z) > 0;
                               1

Would it be possible to improve the handling of the assumptions on Re/Im though?

In Maple 17, the following expression needs to be integrated with respect to q3, p3 and q. Here, mu is a real, positive scalar. 

a := 1/(sqrt(mu^2+(px-p3x-q3x)^2)*sqrt(mu^2+(-p3x+qx-q3x)^2)*sqrt(mu^2+q3x^2)*(sqrt(mu^2+(-p3x+qx-q3x)^2)+sqrt(mu^2+q3x^2)))

However, the integration will not work with the "int" command (e.g. wrt q3). The indefinite integration will work if the integral is evaluated using the steps: highlight expression -> right click -> Integrate -> wrt q3 command.

The output of the integral (using the above method) is very long, it's impossible to manipulate the answer (on my i5, 8GB machine running Maple 17) because it is very tough to copy such a long output. Also, there is no way to specify that mu is a positive scalar. 

Is there a better way to perform the integration, e.g. between 0 and lambda, -1 through 1, or -infinity to +infinity?  

We have the following sequence of natural numbers
1, 2, 4, 7, 11, 16, 67, 83, 46, 73, 47, 85, 70, 20, 16, 76, 83, 55, 73, 56, 85, 79, 119, 934, 463, 389, 1009, 9028, 8237, 7357, 7567, 7688, 8899, 10021, 12035, 53056, 65071, 17093, 39109, 90232, 23249, 94273, 37291, 19316, 61435, 53461, 16481, 18508, 80629, 92657, 75679, 97708, 80831, 13861, 16885, 58916, 62041, 14083, 38099, 99142, 24259, 95303, 30421, 12466, 66485, 58531, 13651, 15698, 89719, 91867, 76889, 98938, 84061, 16121, 12235, 53296, 69311, 11473, 37489, 98552, 25669, 96733, 33851, 15916, 62035, 53111, 11221, 12298, 89309, 90487, 78499, 99578, 87691, 19771, 17885, 58966, 67081, 18173, 37279, 97372, 27479, 97573, 37681, 18776, 67885, 58981, 19091, 19198, 89299, 99407, 70609, 90718, 81821, 12931, 14035, 53156, 65251, 15373, 37469, 96592, 29689, 98813, 32011, 11146, 64235, 53371, 17461, 16598, 89689, 98827, 73019, 91168, 86251, 15401, 10585, 58636, 63821, 12973, 38059, 95222, 22399, 99463, 36641, 14806, 60985, 59051, 15241, 14398, 89489, 98647, 74839, 93998, 90091, 19162, 26345, 54517, 71701, 10874, 47959, 96133, 33329, 92494, 49591, 19757, 75955, 56122, 22331, 13489, 98599, 99758, 85969, 97129, 92351, 15502, 20725, 52877, 78001, 10264, 46379, 97543, 34759, 95924, 43141, 14317, 71525, 52702, 20911, 12089, 98209, 90478, 87599, 99769, 96991, 20162, 26296, 69457, 75692, 29854, 46090, 9263, 3829, 9484, 5051, 1708, 8275, 5933, 3601, 1270, 929, 1138, 8521, 1469, 9853, 3802, 2297, 8137, 7534, 4574, 4972, 3013, 3323, 3454, 4765, 5897, 8209, 9253, 3755, 5800, 313, 542, 475, 805, 740, 280, 316, 848, 1084, 5038, 8543, 3697, 8203, 3269, 9865, 5932, 2639, 9607, 7315, 5384, 5083, 4054, 4754, 4825, 5536, 6608, 8320, 493, 650, 313, 571, 434, 694, 757, 1019, 9364, 4903, 3359, 9799, 10246, 64469, 96715, 52039, 93296, 69511, 11869, 97085, 58354, 45661, 16931, 14239, 93520, 2819, 9463, 3931, 1676, 7045, 5692, 3251, 1810, 469, 1253, 3811, 1474, 5033, 3598, 9247, 7724, 4573, 4051, 1802, 2380, 1132, ... .
 What is the next term? Is it possible to find that with Maple? The Lagrange polynomial is not taken into account. I'd like to recall a great answer by Carl Love to a similar question.  However, the current situation seems to be different. The Predict command of Mma fails here.

Hi!

I have solved an ODE and defined some physical values with Units.

Now i want to calculate the speed at a time of 120s for example, but the evalustion does not show a result.

I want to have a result like ve=63m/s for example and  not ve=v(120s)

Further more if calculating the complete worksheet, the diagramm settings are set to default and i can not see the curve like before. Is there a way to keep user defined settings?

Has somebody an idea what's going wrong?

Calculating completely without units works.

Thank you in advance

Volker
 

Download Berechnungen.mw

Hi

I use some commands like the evalf[10]

But they don't work for the definite double integral in attached file.

It leads to "Error, (in type/algfun) too many levels of recursion"

Please help me to find the answer

Thank You Very Much

Double-Int.mw

Hi,

is there a way to control the way axes tickmark numbers are formated, especially number of decimal places? I have 3 consecutive plot comands with the same options for similar functions in the same x range (and similar y range) and the thickmarks on each plot have differently formated numbers on axes:

1)

2)

3)

in each case I'm calling:

plot(func..., opt2d, labels = [x, U...])

the opt2d are:

fntsz2 := 14;
opt2d := (x = 90 .. 160, size = [500, 200], font = ["Arial", fntsz2], labelfont = ["Arial", fntsz2], labeldirections = [horizontal, vertical]);

I'd like to get all the plots formatted like in the 3rd case: x axis: ticks starting at min (90), numbers rounded to integers; y axis: "0" before decimal point, 2 decimal places.

regards

Wojtek

singular(csc(z), z);
                          {z = _Z1 Pi}
getassumptions(_Z1);
       {_Z1::(AndProp(integer, RealRange(1, infinity)))}

discont(csc(z), z);
                            {Pi _Z2}
getassumptions(_Z2);
                         {_Z2::integer}

I'm not sure why singular adds the condition _Z1>=1.

Assume I had a 2D line

how to put and draw this line into a new geometric world defined by patch?

 Can anyone help me with this problem. I am trying to solve it and its given me an error message "Error, (in dsolve/numeric/bvp) cannot determine a suitable initial profile, please specify an approximate initial solution".

Download Numerical.mw

 Looking forward to your favorable response.

Kind regard

In Maple 2017.2:

if true then print(1);print(2);print(3);print(4);print(5) end if;
                               1
                               3
                               5