One can use display command to combine several plots. For example,

plots:-display(
    plot(PDF(X,t),t=-3..3,thickness=3),
    plot(PDF(X,t),t=0.5..3,filled=true),
    plottools:-line([-1,0],[-1,PDF(X,-1)]));

Alec

binomial(7,6)*binomial(5,3);
                                  70
use combinat in 
numbcomb(7,6)*numbcomb(5,3) 
end;
                                  70

Alec

30/366.;
                            0.08196721311
evalf(30/366);
                            0.08196721311

Alec

combinat:-numbperm(4,2);
                                  12

Alec

Also, 3d plots in Maple can be exported in the AutoCAD DXF format, which can be imported by Solidworks.

Theoretically speaking. It may not work in practise - last time I checked that (few years ago, with Maple 10) , I couldn't open those Maple produced DXF files even in AutoCAD.

Alec

The speed is not the main problem here. The entire procedure seems to be wrong.

Alec

The classical example - Lorenz attractor, can be seen on that page.

Also, the MapleSim version is available from the Maple Application center.

Alec

See ?matrixplot and ?plot3d,option

Alec

Opening a worksheet in the same folder is not enough. The currentdir should be set to it. See ?currentdir

Or you can use the path to that file, "C:/something..." in Windows, or "/.../" in Linux.

Alec

Save it in a library (repository), and you'll be able to use it from other worksheets, just by adding your library (repository) location to the libname.

Alec

For example,

expr:=2*a^32*b^3 + 3*c^2*a^3 + 4*a^5*c^2:

applyrule([a^n::posint=a,b^n::posint=b],expr);

                                       2
                            2 a b + 7 c  a

Alec

I found a way around the infinite loop bug in the applyrule:

`convert/ET`:=proc(expr) 
    local H;
    eval(applyrule(hypergeom(
            [a::algebraic, b::algebraic],
            [c::algebraic], 
            z::algebraic) =
            (1-z)^(-a)*H([a,c-b],[c],z/(z-1)),expr), 
        H=hypergeom)
end;

Alec

The conversion can be done using the following procedure,

`convert/ET`:=expr->subsindets(expr,
    'specfunc(anything,hypergeom)', 
    f->`if`(nops(op(1,f))=2 and nops(op(2,f))=1,
        (1-op(3,f))^(-op([1,1],f))*
            hypergeom([op([1,1],f),
                op([2,1],f)-op([1,2],f)],
                op(2,f),
                op(3,f)/(op(3,f)-1)),
        f));

For example,

convert(hypergeom([a,b],[c],z),ET);

                   (-a)                              z
            (1 - z)     hypergeom([a, c - b], [c], -----)
                                                   z - 1

simplify(convert(%,ET),symbolic);

                      hypergeom([a, b], [c], z)

convert(hypergeom([1,2],[3],z),ET);
                                             z
                    hypergeom([1, 1], [3], -----)
                                           z - 1
                    -----------------------------
                                1 - z

Alec

It should be possible to implement that Euler transformation by defining something like `convert/ET` using applyrule. However, when I tried to use applyrule for that, it went into a seemingly infinite loop. It looks like a bug.

Here is a more simple example of that bug. First - when applyrule works as it should,

applyrule(sin(x::algebraic)=cos(x/2),sin(3*t));

                                   3 t
                               cos(---)
                                    2

Now, replace cos with sin,

applyrule(sin(x::algebraic)=sin(x/2),sin(3*t));

An infinite loop. It looks as if applyrule is applying recursively to the result - producing first sin(x/2), then sin(x/4), sin(x/8) etc. Normally, such things shouldn't happen.

Alec

Use other than e index. For example,

for i in E1 do i od;

Alec