@Anthrazit 

You should have attached this file when you asked your question, this would have saved me from taking off in the wrong direction and wasting my time.

@Anthrazit 

That doesn't change anything: your code still runs without errors with Maple 2015:

Download 2015_a.mw

Maybe you should add this line after the last local declaration

print(detail(segm));

to see what segm really is.
You should get this`

[name of the object          segm
form of the object           segment2d
the two ends of the segment  [[0,0],[1,1]]

that nothing is wrong BEFORE these lines?
For instance that you didn't write

segment(segm[A,B]):  # no comma after sem

?

It works well in Maple 2015

Download 2015.mw

@Carl Love 

Just for fun: I suppose you knew Ramanujan's anecdote about the number 1729?
I not here it is https://vedicmathschool.org/story-of-hardy-ramanujan-number-1729/

@vv 

You're right, I made a confusion between the elements of the sigma-algebra, which are sets, and the set of elements from which this sigma-algebra is build.

@Zeineb

Taking into account @vv's remark, here is an edited version of the code I delivered previously.
I believe there is a simpler way to build the sigma algebra generated by a subset S of a set of sets.
Some examples are given in the attached file. 

Download SigmaAlg_b_mmcdara.mw

@vv 

I don't understand.
A sigma-algebra T(S) of a set S is a subset of P(S) (the power set of S), whose cardinal is of course a power of 2.
But there is no reason for the cardinal of T(S) to be a power of 2 too.
Am I right?

In any event I didn't verifiy if the set produces by ALG was indeed closed under complementation and countable union, nor if it computes correctly the complement of a set of the form {p <= z < q} (which I doubt for, as you quoted, ALG is essentially written for finite sets).
This is why I adviced the OP to check if its procedure ALG works correctly for non countable sets.

@Saha 

the error message
Warning, expecting only range variable t in expression sin(x)-sin(x)*t+1/2*sin(x)*t^2-1/6*sin(x)*t^3 to be plotted but found name x

means that you want to plot wrt t something which depends on t and x.
So either give x a value or use plot3d (see the help pages instead).

Have you observed that the function you want to plot is of the form

f(t, x) = sin(x)*g(t)

?
So why not trying to plot simply g(t) ... which is the beginning of the Taylor expansion of exp(-t) at t=0 ? 
Look to these plots

plot3d( sin(x)*mtaylor(exp(-t), t=0, 4) , t=-5..5, x=-2*Pi..2*Pi);
plot3d( sin(x)*exp(-t) , t=-5..5, x=-2*Pi..2*Pi);

@NanaP 

There are always questions to which several players contribute and others that remain more "confidential".

You'll quickly understand all this with a little practice.

The most upsetting thing (even for me) is that you'll never know why your question remains unanswered (is it that no one has an answer to offer, or is it that the problem is so uninteresting or so poorly presented that no one cares to answer it?)

So a good practice is probably to present your problem in the most consise, but nevertheless detailed, form and to deliver what you have been able to achieve with Maple (generally people don't like an image or a reference, or a piece of some article to figure out what the OP [the acronym used to designate the person who asked the question] wants).

Concerning successive questions, I agree that when no one have given you a reply, it's very tempting to bring up a question again. But you have to account for jet lag and half a day may separate a question and its answer.

Concerning the difficulty to find Maple solutions or algorithms in the web, I agree that it is extremely difficult.
Theitems you get refer generally to Mapleprimes. An advice: even if the Search/Advance Search feature could be improved, do not hesitate to browse amonq the Questions and Post to seek if similar problems have not been alerady answeerd.
Another interesting link is https://www.maplesoft.com/applications/ 

Hope to see you soon.

To conclude, I believe, this thread, here is a more complete version of BDS_2.mw where I added several estimators of the left and right endpoints of the underlying distribution, given a sample from it and the lnowledge that it is left and right bounded.
Tou will see that the estimators are quite good is the distribution is "smooth" and quite rough if it is sharp (here left peaked).
BSD_2b.mw 

If you'd like to dig up further on the subject, here are some companions articles (their understanding resquires some knowledge in Extreme Values Theory).​​
https://boris.unibe.ch/36886/1/10687_2007_Article_52.pdf
https://arxiv.org/pdf/1412.3972.pdf
https://projecteuclid.org/journals/annals-of-statistics/volume-12/issue-4/Estimating-an-Endpoint-of-a-Distribution-with-Resampling-Methods/10.1214/aos/1176346811.full
https://www.d.umn.edu/~yqi/downloads/19sinica.pdf
https://mendel-journal.org/index.php/mendel/article/view/28/30

@NanaP 

UPDATED
(simpler definition of the SubjectiveBeta Distribution)

COMMENT 1
You should stop posting replies while the previous one doesn't receive an answer.

COMMENT 2:
Note from @NanaP: The Use_formulas.mw file by @mmcdara stopped at the 25 percentile "theoretical" calculation and did not go beyond this point until a histogram of simulated values is created

You're partially wrong: the first file I sent you (BDS.mw) "goes as far as the histogram", maybe you didn't pay attention?

COMMENT 3: (in response to your next reply)
My file Use_formulas.mw contains indeed some errors: I used the moments to compute the Skewness and Kurtosis instead of the central moments.
Apologies.

Here is a file which "[achieve] the ultimate goal of being able to simulate from a custom defined distribution of SB(minimum, mode, mean, maximum) .",  although I didn't particularly like the way you worded it.

The procedure SubjectiveBeta builds a "Subjective Beta Distribution" given its mean, mode, min value and max value

SubjectiveBeta := proc(me, mo, P, Q)
   # me = mean
   # mo = mode
   # P  = xmin
   # Q  = xmax

Download BSD_2b.mw

I hope that's all right with you, because I'm done with this subject.

@NanaP 

I'm not sure I understood correctly so correct me if I'm wrong:

• You have a random variable X with a SB (short for "subjective Beta") distribution of parameters xmin, xmax, alpha, beta (I'm not sure of this point, maybe they are xmin, xmax, mode, mean ?).
• You have a sample S of X, let's say of size N, from which you compute some empirical statistics, for instance
  the mean 

  m := (1/N)*add(s[n], n=1..N)
  

• Some facts are:
  • the above empirical statistics is the best unbiased asymptotic estimator of the Mean(X) (of course you already knew that);
  • the classical estimation of the variance 

    v := (1/(N-1))*add((s[n]-m)^2, n=1..N)

    also has this same property.

  • The case of Skewness and Kurtosis is a little bit more complex. If I'm not mistaken there are no formulas for unbiased estimator (nearly unbiased estimators do exist).

  • The case of xmin and xmax is even more complex.
    I have some material concerning what is called "the endpoint distribution", that is the estimation of xmin or xmax from sample S.
    Here is a first reference, https://arxiv.org/pdf/1412.3972.pdf, do not hesitate to ask for some more.

  • Concerning the estimation of the mode there is no unbiased estimator, but a lot of estimators do exist.

@NanaP 

Thank you for introducing me to something I didn't know about, even though I've been working in a parallel domain for a long time.

Hope my contribution will be useful

@acer 
I had observed on my own that isolate didn't work for this example.
So you are right saying that the imaginary unit is not the cause of my difficulty in the sense that I could faced them without having I.

It's just that for the particular case b=1 isolate works if c <> I, and not instead, so my quick jump to an erroneous conclusion.

@acer 

Thank you acer for all these informations.

@NanaP 

The case 

Min := 1000; Mlikely := 1400; Avg := 1500; Max := 2100;

can be solved by using the definitions of Mean, Variance, Skewness, Kurtosis, Mode and Quantile:

Download Use_formulas.mw

Note that the mode is not necesarily  unique (Beta(a, a) has two modes at x=0 and x=1 when a < 1 and a continuum of modes whan a=1).
So, I don't know what you want to do, but this re-parameterization you want to use is extremely dangerous.

Download Distribution-Beta-Subjective_Failure.mw

About the Beta distribution there are only, up to my knowledge, these parameterizations:

• shape (a) and scale (b), the most common ant the one Maple adopts;
• mean (M) and variance (V);
• mean and precision (inverse of variance, noted tau below);

As you can see there is a 1-to-1 map between each of pair of parameterization

m := Mean(Beta(a, b)):
v := Variance(Beta(a, b)):
solve({m=M, v=V}, {a, b});

solve({m=M, 1/v=tau}, {a, b});
        /        / 2        \      / 2        \        \ 
        |      M \M  - M + V/      \M  - M + V/ (M - 1)| 
       < a = - --------------, b = -------------------- >
        |            V                      V          | 
        \                                              / 
  /      3        2              / 2                \        \ 
 { a = -M  tau + M  tau - M, b = \M  tau - M tau + 1/ (M - 1) }
  \                                                          / 

This VOSE article you refer about seems to be oriented to risk assesment, for which this Subjective Beta Distribition might be, according to the authors, and likely under some restrictions an interesting tool to achieve this goal.
I'm working on the domain of Reliability Analysis (assessment of probability of failures due to extreme events in engineering, something probably not that far from the "risk assessment" VOSE talk about?]) for more than 25 years and I had never heard of the approach developped by VOSE (but each domain develops it's own method with little inter-operability).

Some intersting links:

Probability ontology
http://www.math.wm.edu/~leemis/2008amstat.pdf
https://en.wikipedia.org/wiki/Relationships_among_probability_distributions

Different distributions
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5013898/bin/supp_btw170_Swat_et_al_Supplementary1_ProbOnto_KBprintout.pdf

Different parameterizations
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5013898/bin/supp_btw170_Swat_et_al_Supplementary2_ProbOnto_useInPharmML.pdf

@NanaP 

Does it mean you want to re-parameterize a Generalized Beta distribution by min, max, mode, min instead of alpha, beta, min and max ?

If it so you will have a problem because the relation (alpha, beta) --> (mode, mean) is not a 1-to-1 mapping.
Take the simplest case alpha > 1, beta=alpha, xmin=0, xmax=1 (thus the Generalized Beta is just Beta(alpha, alpha)).
As Mean(Beta(alpha, alpha)) = Mode(Beta(alpha, alpha)) = 1/2 whatever alpha > 1, a re-parametrization of the Beta distribution in term of mode and mean doesn't say anything about the pdf, cdf, variance, ... of this distribution.