In a previous post, I promised to write about testing the quality of pseudo-random number sequences.  I'll post later about some of the statistical tests often used, but I first wanted to mention a sort of practical test one can do. One of the many things you might want to do with pseudorandomly generated numbers is Monte Carlo integration/simulatation/etc.  As mentioned by acer in this comment, Monte Carlo integration can be shown to work better with some of the pseudorandom number generators (PRNGs) which are considered inferior in a statistical sense.  In this post, we will play with a simple Monte Carlo approximation of π.

I am performing reliability analysis.  I have two types of data:

1.  Suspended data (data that has not failed)

2.  Failed data

I require to calcaulte how many ways the data can fail for example:

1. F1   Time = 10

2. S1  Time = 100

3. F2  Time = 200

4. S2  Time = 250

The number of combinations for this data is 3:

F1, S1, F2, S2

F1, F2, S1, S2

F1, F2, S2, S1

There are two pieces of extended functionality that I quite often want from the Maple Compiler. The first (task A) is to be able to link in and use an arbitrary function from some other external ("3rd party") shared library, within my Compile'd Maple procedure. The second (task B) is to directly call the compiled Maple procedure from within some computational routine in a 3rd party shared library (which I would then access using define_external). This post is about the first of those, task A.

Hello all,

I'm working through an optimization problem, and I'm running into a problem adding constraints on the variables I am trying to optimize. I'll attempt to ask the question in a general sense, but if it's not clear, I'll be more specific.

I have a procedure called EIG that calculates the sum of the squared residuals of a model compared to experimental data by inputting 6 variables. I'm trying to minimize this function by using the Optimization...

If you work with Maple on a daily basis sometimes you come across
things that seem impossible to do. Yesterday I tried to have two colors
in a histogram ie -3..-1 = red and -1..3 =blue but I could not figure
out how to do it...hummm.....any suggestions?

restart:
with(Statistics):
with(plots):

n:=5000:
A:=Sample(RandomVariable(Normal(0,1)),n):

Histogram(A,bincount=100, color=blue);

help_for_maple_quest.docplease give a explain how to work for operations.

 I would like to select the color of the arrow for T,N and B.

Thank you for any suggestions.

Ralph Wallace

restart;
with(VectorCalculus):with(plots):
Display the tangent, principal normal and binormal vectors on a curve.
c:=sqrt(a^2+b^2);
                               ...

Hello all,

I have a little problem with an expression I have to develop in taylor/series.

Here it is:

I want to develop it in "nu", giving that "nu" is a small parameter, and it says that it doesn't exist a series expansion.

But when I can obtain the series expansion of the "sqrt" part only of the expression, which yields:

E.g. I want to modify the Binomial distribution in a way, that k is not the number of successes, but the outcome.

Let's assume each success wins $100, while a failure results in $0.

My thought was to substitute k for k/100, but that produces an error message.

Can anyone tell me, how to implement that in Maple?

Thanks a lot.

Hi,

i use fsolve to find appraxmaite solution. there are some parameter in my equastion wher i do loop in this pramater form 0.. 1000,000,00

i have other problem with "fsolve and digits". when Digits =12;  i cannot find  approx solution using "fsolve"

when i decrease or increase  the digit "Digits:=11 , Digits:=13",  i find the  approx solution using "fsolve".

i donot know what is wrong!!!

Hi all,

I have Maple 11  ( Student Version - I am a UC Berkely Distance Learning Student ). I need the Dynamic Systems Package, can I down load it and install it in my Maple 11 install ?? If so how ??

Thanks,

Rajat Mitra

There was an interesting question from Christopher2222:

Maple keeps changing the variable xi to the greek letter

Now I need a solution for a specific application, Here is the code:

with(Statistics):

L:=[Normal, Beta, Gamma]:

for i to 3 do
f[i]:=RandomVariable(L[i](4,2):
end do:

This will work for L[1], but not for L[2] and L[3]...

hello all, I'm trying to solve ODE's but my computer takes a very long time evaluationg. and still not giving an answer for more than 4 hours.

so please is there anything error in my solution

deqs := {NULL
         , diff(u2(y), y, y)-((h^2/K)+m*h^2*M^2)*u2(y)+m*b*rho*h^2*GR*(c3*y+c4) = 0
         , diff(u1(y), y, y)-((1/(K*(1+k1)))+(M^2/(1+k1)))*u1(y)+k1*(diff(N(y), y))/(1+k1)+GR*(c1*y+c2)/(1+k1) = 0

I am not able to format the number along the axis of a 2D plot to use period as the decimal separator - and not to have different number of positions after the decimal on the x- vs. the y-axis either.  I have set the Windows language/region to the US.  Anyone who might help?  Regards,