Hi,

Although everyone one might know, but still I want to clarify the concept of Sensitivity analysis. Suppose we have an objective function Z in which we have two variables X & Y and two parameters a and b. So in sensitivity analysis, we systematically change the values of parameters ( for eg. changing value of "a" by 20%, then 40%, then 60%, similiarly for parameter "b"). At a time, we only change the value of one parameter and the values of other parameters remain fixed to their original value. Then we would again solve the model ( ie. objective function Z) and calculate the new values of objective function and variables.  Then we try to see the impact of parameter change on the value of objective function Z and on variables X & Y.

The formula for sensitivity analysis for parameter a is  for objective function Z =

 (new values of objective function  - old value of objective function)*100 / (old value of objective function)

Here new value of objective function Z denotes when we use changed value of paramter a (keeping other parameter values fixed to their original value) and old value of objective function denotes value of Z with all parameters (including parameter a) at their original value

Although I can do it by a for loop but it is not an efficient procedure and also time consuming. Can anybody suggest whether any package suitable for sensitivity analysis is available in maple.

Thanks for your attentive reading of problem.

Thanks and Regards,

Nilesh

Is there any command or package in Maple that computes blow ups? I didn't find blow up in the help.

Please simplify the expression in image.Kind Regards,

Hi Maple friends,

I am trying to run a script I have written in Maple 2016 directly from my Terminal on my MacAir. I have OSSierra 10.12.1 and Maple 2016.

My script needs to read a file in and then run, whilst writing to some files on the way. I need this, since I want my script to be run on a large computer on my institute, where I have no interface/Maple open and has to access data on the way.

So, I go to Terminal, and stand in the directory of where the file/scrips is (lets call it terminal_test.mpl OR terminal_test.txt) and then write:

and get back

I also get this if I use .txt file AND if I write "Maple 2016" instead of ./maple. I have tried to write 

and then Maple opens in the right file, but does not run it. This is just to show, that my terminal can find the program. But this is not what I want. As I said, I would like it to:

start running a script from Terminal - read a file to be used in the script - on the way, write to 6-8 different files, that I can then access whilst the script is running. So my questions are: 

1. How do I run the script/file in the first place, and am I in the right place?
2. Where do I place the files to be read during the running of the script (Reading to script works fine already, so it should work also if I can just run the script/file)?
3. Where do the files I create in my script end, so I can access them?

Thank you so much for any answer - I have not been able to find anything on the net so far that could help me, but I know you guys are really good on this.

Tomas.

We have just released the 3rd edition of the Mathematics Survival Kit – Maple Edition.

The Math Survival Kit helps students get unstuck when they are stuck. Sometimes students are prevented from solving a problem, not because they haven’t understood the new concept, but because they forget how to do one of the steps, like completely the square, or dealing with log properties.  That’s where this interactive e- book comes in. It gives students the opportunity to review exactly the concept or technique they are stuck on, work through an example, practice as much (or as little) as they want using randomly generated, automatically graded questions on that exact topic, and then continue with their homework.

This book covers over 150 topics known to cause students grief, from dividing fractions to integration by parts. This 3^rd edition contains 31 additional topics, deepening the coverage of mathematical topics at every level, from pre-high school to university.

See the Mathematics Survival Kit for more information about this updated e-book, including the complete list of topics.

eithne

Suppose I want to revolve the curve given in the (X,Y)-plane by the set of parametric equations

x(t) = cos(t) + t sin(t)

y(t) = sin(t) - t cos(t)

for t in [0,Pi/2] around the X-axis. How can I plot the given surface of revolution? Similarly, the same question for

x(t) = exp(t)*cos(t)

y(t) = exp(t)*sin(t)

for t in [0,Pi/2]

How to plot an Elliptical Function of Third kind complete or incomplete, eg. EllipticPi(n,k) if n and k are not constants?

as The function I wish to plot and explore contains Elliptical Function of Third kind complete and incomplete with complicated form of n and k.

Please reply asap.

Regards

Hi everybody,

I have differential equation of the form eq := diff(x(t), t) = x(t) + 1  and I want to find the name of the dependent variable (more generally this issue concerns a differential system in N dependent variables).

How should I use “indets” in order to obtain {t, x(t)} ?

(indets(eq, name) returns {t}, while indets(eq, function) returns {t, x(t), d/dt x(t)}.

Thanks in advance

Perhaps I am being stupid, but I have a problem with functional operators, more specifically with looping through a collection of them. Consider the following code (just a study case):

nonIdMaps := []:
for f in [x -> x,x -> 2*x,x -> 3*x] do
   f,f(y);
   if f(y) <> y then nonIdMaps := [nonIdMaps[],f] end if
end do;
nonIdMaps,map(f -> f(y),nonIdMaps);

There are two things I do not understand about this output: 1.) Why the f's rather than the specific maps as given in the list? 2.) Why two times 3*y in the last line, rather than [2*y,3*y]? (this issue being of course quite likely a consequence of the first one). This behaviour can be compared with the following:

nonIdMaps := select(f -> f(y) <> y,[x -> x,x -> 2*x,x -> 3*x]):
nonIdMaps,map(f -> f(y),nonIdMaps);

the output of which I do understand, being as I expected it to be. I guess that there are some scoping rules I do not understand, or what?

how would you interpret the solutions to this:

>int(sin(y/2)^2/(x*(x-y)*y^2), x=omega__ir*t..omega__c*t) assuming t::positive, omega__ir::positive, omega__c::positive, omega__c>omega__ir;
 

which leads to the expression shown in the screen shot. In particular, I'm interested in the condition for the solution to be "undefined"

thanks.

This app shows the calculation of the final speed of a body after it made contact with a variable force; Taking as reference the initial velocity, mass and graph of the variation of F as a function of time. That is, given the variable forces represented by the lines in the time intervals, we will show the equation of momentum and momentum; With their respective values, followed by their response.

Momentum_with_two_variable_force.mw

https://www.youtube.com/watch?v=2xWkF7JhkpI

Lenin Araujo C.

Ambassador of Maple

If I make a single plot it will always be located in the centre of the worksheet. Is there an option to place it at the left border?

Dear friends,

I have to solve the BVP of the following type:

(x-1)*y''(x)+y'(x)=1+x, y(0)=1, y(1)=1

Expressing the highest derivative (y''(x)), we get a singularity at x=1.

So, as recommended, the method bvp[middefer] (or bvp[midrich]) is used:

dsolve({(x-1)*diff(y(x), x$2)+diff(y(x),x) = 1+x, y(0)=0, y(1)=1}, type = numeric, method = bvp[middefer], y(x), 'output' = Array([seq(k/10, k=0..10)]), 'abserr'=1.0e-3, 'maxmesh'=100)

Maple generates the following error:

Error, (in dsolve/numeric/bvp) unable to achieve requested accuracy of 0.10e-2 with maximum 100 point mesh (was able to get 0.82e3), consider increasing `maxmesh` or using larger `abserr`
But as I increase maxmesh, the achievable accuracy (reported in "(was able to get XXX)") decreases:

for

maxmesh = 200: 0.19e4

maxmesh = 500: 0.54e4

maxmesh = 1000: 0.12e5

maxmesh = 5000: 0.73e5

and so on.....

How to solve this equation with desirable accurancy and usage of Array as output?

i want to solve an ode , but maple return an integral in result, how can i have an answer?

Download dsolve.mw

 Dear Maple users,

In Maple 18, I want to  label y axis as $\hat{\sigma}_y$.

labels = [x, sigma[y]] works fine but I have no idea how to put a hat on sigma[y] .

regards,