hey everyone
i got a little problem here i really dont get it why my fsolve give me this error

Error, (in sqrfree) argument must be a polynomial or a rational function in {cf[-2], cf[-1], cf[0], cf[1]}
these are my variables :

cf[-2], cf[-1], cf[0], cf[1]}

thats part of my code :

res1;
(34.30563197 cf[-1] - 10.13498047 cf[0] + 5.197134649 cf[1]

   - 0.4714805434) (0.3095515346 cf[-2] + 0.3822521253) - (cf[-1](2.466022067

  ) + cf[0](-0.6605923783) + cf[1](1.076415647))^2

   - cf[-1](2.466022067) - cf[0](-0.6605923783)

   - cf[1](1.076415647) - 333.1166471 cf[-1] + 186.8672744 cf[0]

   - 128.6145289 cf[1] + 0.8737683108, (-39.73727883 cf[-1]

   + 26.25682759 cf[0] - 7.289811219 cf[1] - 0.3506934780)

  (0.4161980514 cf[-1] + 0.4402863507) - (cf[-1](0.3227083347)

   + cf[0](2.517826949) + cf[1](-0.7889944208))^2

   - cf[-1](0.3227083347) - cf[0](2.517826949)

   - cf[1](-0.7889944208) - 15.03436409 cf[-1]

   - 129.5665191 cf[0] + 68.25827819 cf[1] + 0.5888637790,

  (17.36111111 cf[-1] - 29.80788311 cf[0] + 17.36111111 cf[1]

   - 0.2500000000) (0.5000000000 cf[0] + 0.5000000000)

                                                            2
   - (cf[-1](-1.833333334) + cf[0](0.) + cf[1](2.333333334))

   - cf[-1](-1.833333334) - cf[0](0.) - cf[1](2.333333334)

   + 38.17276297 cf[-1] + 34.00261850 cf[0] - 77.23526300 cf[1]

   + 0.3750000000, (-3.548332053 cf[-1] + 12.78057007 cf[0]

   - 19.34221012 cf[1] - 0.1707008416) (0.5290914191 cf[1]

   + 0.5597136492) - (cf[-1](1.060931345) + cf[0](-1.540306477)

   + cf[1](-0.2297630675))^2 - cf[-1](1.060931345)

   - cf[0](-1.540306477) - cf[1](-0.2297630675)

   - 20.59209188 cf[-1] + 28.27859544 cf[0] + 40.45630289 cf[1]

         -10                      
   + 1 10    cf[-2] + 0.2254717519

fsolve({seq(res1[v] = 0, v = 1 .. 2*N+2)})     ( N here is 1 )

why i cant get the awenser and if you can plz solve it for me

i would be happy that anyone can give me an explanation with details

tnX a lot

Arian
 

I'm trying to figure out how to plot the following: if I have a region G in R^2, and two functions f <= g on G, I would like to plot the projection on the XY-plane in say one colour and the volume between f and g in another. For example, consider G the region between the parabolae y=x^2 and y=2-x^2, and f(x,y) = x+y+4 and g(x,y) = 25-x^2-y^2, just to name something. I'd also prefer the plot to be easy adaptable to other functions, e.g. with a different region G, different function descriptions for f and g but also, if possible, projection on one of the other two coordinate planes. Is there an easy way to do this?

I wonder - would it be possible to automate the following way of adding a legend to plots? If so - how?

plot([f(x), g(x), 7], legend = ['f(x)' = f(x), 'g(x)' = g(x), y = 7]);

In my opinion such a legend makes the plot much more readable - but most students (and others) will usually be too lazy to type this out, hence the wish for it be to automated.

Any help will be much appreciated.

Hello,

I would like to ask you for a help with solving a differencial equations for oscillator damped by a constant friction force (no damping force proportional to the velocity!)

General solultion:

Where 

The amplituce A is decreasing by C at each turning point and center of the ascillatory motion moves: x = -C/2 or x = C/2

Moving:

where

How to create a graphic output of movement according to equation x(t) similar to graph mentioned below?

Thank a lot for your help. I have been working on this task for several weeks without success.

Tom

Hi,

Although i can evaluate speed, displacement and accelleration, i can't plot these functions.

The error message is:

What am i doing wrong in my definitions?

Thank you in advance.

Download Berechnungen.mw

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