> restart;
> with(plots); with(StringTools); with(plottools);
> INF := 999999999999999999999;
                     999999999999999999999
> NULL;
> MinoxAngle := 200; MikromaAngle := 350; MinicordAngle := 290; GamiAngle := 280; GamiFocal := 0.25e-1; SummitarDial := [1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 2, 2.2, 3, 4, 5, 6, 7, 10, 20, INF]; Minox35MLDial := [3, 4, 6, 10, 20, INF]; Minox35Angle := 100; MinicordDial := [.35, .4, .5, .6, .7, .8, .9, 1, 1.2, 1.5, 2, 3, 4, 8, INF]; Mini := nops(MinicordDial); MikromaDial := [.5, .6, .7, .8, .9, 1, 1.2, 1.5, 2, 2.5, 3.5, 6, INF]; MinoxLXDial := [.2, .24, .3, .4, .6, 1, 2, 30]; MinoxLXAngle := 270; GamiDial := [.5, .6, .7, .8, 1, 1.2, 1.5, 2, 3, 5, 99990000000000]; MinoxBDial := [8*(1/12), 10*(1/12), 1, 1.5, 2, 3, 6, INF]; MinoxBAngle := 270;
                              200
                              350
                              290
                              280
                             0.025
[1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 2, 2.2, 3, 4, 5, 6, 7, 10, 

  20, 999999999999999999999]
            [3, 4, 6, 10, 20, 999999999999999999999]
                              100
 [0.35, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.2, 1.5, 2, 3, 4, 8, 

   999999999999999999999]
                               15
    [0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.2, 1.5, 2, 2.5, 3.5, 6, 

      999999999999999999999]
              [0.2, 0.24, 0.3, 0.4, 0.6, 1, 2, 30]
                              270
   [0.5, 0.6, 0.7, 0.8, 1, 1.2, 1.5, 2, 3, 5, 99990000000000]
         [2  5                                        ]
         [-, -, 1, 1.5, 2, 3, 6, 999999999999999999999]
         [3  6                                        ]
                              270
> 
> NULL;
> dd := GamiDial; N := nops(dd); dstx := [seq(convert(dd[i], string), i = 1 .. N)];
   [0.5, 0.6, 0.7, 0.8, 1, 1.2, 1.5, 2, 3, 5, 99990000000000]
                               11
  [".5", ".6", ".7", ".8", "1", "1.2", "1.5", "2", "3", "5", 

    "99990000000000"]
> NULL;
> MinicordAngle := 290;
                              290
> NULL;
> 
> LensDial := proc (LensName, focal, Angle, scale, dr) local p1, p2, p3, p8, pk, pt, rk, R, R2, R3, Rc, c1, ds2; R := 1600; R2 := 1400; R3 := 1200; Rc := 1500; CaptionTail1 := "EXTENSION ANGLE"; CaptionTail2 := "LENS FOCUSING DIAL"; Caption1 := Join([LensName, CaptionTail1]); Caption2 := Join([LensName, CaptionTail2]); q := 180/Pi; rotation := 90; dir := dr; ds := scale; N := nops(ds); dstx := [seq(convert(ds[i], string), i = 1 .. N)]; ds2 := subs(dstx[N] = infinity, dstx); MaxAngle := Angle; f := focal; degr := -(-ds[1]+f)*Angle/(D-f)+rotation; g1 := degr/q; for j to N do deg[j] := subs(D = ds[j], degr) end do; for i to N do rdn[i] := evalf(deg[i]/q); xv[i] := R2*cos(rdn[i]); yv[i] := R2*sin(rdn[i]); wx[i] := R3*cos(rdn[i]); wy[i] := R3*sin(rdn[i]) end do; pk := {seq([ds[i], deg[i]], i = 1 .. N)}; rk := {seq([dir*xv[i], yv[i]], i = 1 .. N)}; txt := {seq([dir*wx[i], wy[i], ds2[i]], i = 1 .. N)}; p3 := pointplot(rk, caption = Caption2, captionfont = ["ROMAN", bold, 22], symbol = solidcircle, symbolsize = 15, color = red, axes = none); c1 := circle([0, 0], Rc, thickness = 8); p8 := textplot(txt, 'font' = ["times", "bold", 14]); display(p3, c1, p8, scaling = constrained) end proc;
Warning, `CaptionTail1` is implicitly declared local to procedure `LensDial`
Warning, `CaptionTail2` is implicitly declared local to procedure `LensDial`
Warning, `Caption1` is implicitly declared local to procedure `LensDial`
Warning, `Caption2` is implicitly declared local to procedure `LensDial`
Warning, `q` is implicitly declared local to procedure `LensDial`
Warning, `rotation` is implicitly declared local to procedure `LensDial`
Warning, `dir` is implicitly declared local to procedure `LensDial`
Warning, `ds` is implicitly declared local to procedure `LensDial`
Warning, `N` is implicitly declared local to procedure `LensDial`
Warning, `dstx` is implicitly declared local to procedure `LensDial`
Warning, `MaxAngle` is implicitly declared local to procedure `LensDial`
Warning, `f` is implicitly declared local to procedure `LensDial`
Warning, `degr` is implicitly declared local to procedure `LensDial`
Warning, `g1` is implicitly declared local to procedure `LensDial`
Warning, `j` is implicitly declared local to procedure `LensDial`
Warning, `deg` is implicitly declared local to procedure `LensDial`
Warning, `i` is implicitly declared local to procedure `LensDial`
Warning, `rdn` is implicitly declared local to procedure `LensDial`
Warning, `xv` is implicitly declared local to procedure `LensDial`
Warning, `yv` is implicitly declared local to procedure `LensDial`
Warning, `wx` is implicitly declared local to procedure `LensDial`
Warning, `wy` is implicitly declared local to procedure `LensDial`
Warning, `txt` is implicitly declared local to procedure `LensDial`
> ;
> NULL;
> LensDial("MEOPTA MICROMA  HELGOR 25mm", 0.25e-1, 350, MikromaDial, 1);

> LensDial("GOERZ MINICORD  25mm", 0.25e-1, 335, MinicordDial, 1);

> 
> ;
> LensDial("MINOX LX MINOX 15mm", 0.15e-1, 270, MinoxLXDial, 1);

> LensDial("GAMI ESAMITAR 25mm", 0.25e-1, 290, GamiDial, 1);

infinity symbol looks like an  "8" lying down.

How to use Maple to plot  the  infinity( lay flat  "8") symbol at 12 o'clock position of a circle

 Hi, I got an error (the code below) when using subsop in procedure. But when I replace "NULL" with an integer, it would work.                             

> phephoi := proc (A, B)
local check, i, j, C, E;
C := A;
E := B;
for i from 1 to nops(C) do
   check := 0; 
   for j from 1 to nops(E) do
      if E[j] = C[i] then check := 1; break; end if;
   end do; 
   if check = 0 then C := subsop(i = NULL, eval(C)); end if;
end do;
return C;
end proc;

> A := [1, 2, 3, 5, 6, 7];
B := [2, 4, 7, 11, 8];
phephoi(A, B);

                             [1, 2, 3, 5, 6, 7]
                              [2, 4, 7, 11, 8]
Error, (in phephoi) invalid subscript selector
 

How to put number 1,2,3...12  around a circle, like a clock face using Maple code ?

Right now I can only do this manually by plotting a circle, then save the file and use Windos paint to

mannual put in the numbers

Hi all

I would like to plot a quadrilaterial (  spherical quadrilaterial)  on the sphere ( the vertices are known)

How can I do this

Many thanks

Let us consider 

with(Statistics);
U := RandomVariable(DiscreteUniform(-10, 10)):
V := RandomVariable(DiscreteUniform(-10, 10)):
Probability(U^2-V^2 <= 1/9, numeric);
  0.

, whereas a positive number greater than 1/21 is expected. 

I have not been using Maple for the last couple of years because I have Maple installed on my big and heavy laptop. Such laptop forces me to sit at a desk when I work which I dont like. This has resulted in that I have only (for the most part) used my ipad for the last couple of years hence dropping my laptop and Maple. The problem with such approach has been that is more difficult to write scientific papers because you cant run a full version of word on the ipad hence you cant use a citation manager. There does not exist an app on the ipad that offer one. This leads to a lot of extra work!

I have therefore now bought a Microsoft Surface Pro 4 which is a tablet (like the ipad) but its a not a "toy" tablet like the ipad or android tablets because it has Windows 10 Pro (64-bit) which means that you can run full versions of pc software like word 2016 and Maple. The version of Surface Pro 4 I bought has an Intel Core m3-6Y30 processor with 2 cores, 4GB RAM, 128 GB of storage, no fan and a processor base frequency of 0.9 GHz. The m3-6Y30 processor has something called dynamic frequency scaling which means (as I understand it) that the processor can go up to 2.2 GHz when it is under heavy load. Since Maple recommend that you have a system with a cpu with at least 1 GHz to run Maple 2016 I am not anticipating any problems here.

http://www.maplesoft.com/products/system_requirements.aspx

However, I am not completly sure if maplesoft means that the base frequency needs to be 1 Ghz or not? Since my tablet has not been delivered yet I cant comment on how Maple runs on it but as I said I am not anticipating any problems.

Now to my question: When will Maple 2017 come out? Since it is already 2017 I dont want to buy an old version of Maple now that I want to update my Maple software.

I want to do matrix operation,vis-a-vis,addition,inverse of a 15 x 15 matrix. I already typed the matrix,but remains the operation.

Thank you

Raji

I have never seen DEtools:-odeadvisor hang before. I've seen dsolve itself hang many times, but not odeadvisor. Is this a bug? Maple 2016.1 on windows.

restart;
eq:=(2*y(x)*diff(diff(y(x),x),x)-diff(y(x),x)^2)^3+32*diff(diff(y(x),x),x)*(x*diff(diff(y(x),x),x)-diff(y(x),x))^3 = 0;
DEtools:-odeadvisor(eq);

Now it hangs. Waited for more than one hr. This is an ODE from a book. Do others see the hang as well?

I can't find in the Linear Algebra package the command for creating the duplication matrix.

The duplication matrix is the n * n rows and n * ( n + 1 ) / 2 columns that transforms vech(A)
to vec(A) where A is a symmetric n by n matrix (Magnus & Neudecker).

If it not yet in Maple itself, has someone already made a proc for it?

Harry Garst

I want  to find the parametrs that verify the below inequalities at the same time:

beta*(v+u)>phi*(alpha+lamda+u+delta)

phi>beta

beta> (alpha+lamda+u+delta)

phi* (alpha+lamda+u+delta)>beta*(u+v+alpha),

question : find the values of u,v,delta,lamda,phi,beta,alpha such as all these parameters are in the intervall ]0,1]

can anyone help me ?

Hi,

I have a second order, linear, non-homogeneous differential equation and for the solution Maple takes the particular solution under a indefinite integral form. After I substitute the values of the coefficients I want Maple to perform the integration. The integration is possible because I individually integrated one small part of the expression. The full expression has a lenghty sumation of different indefinite integrals so it would be cumbersome to perform each integration by hand.

Can somebody help me force Maple to perform these integrations?

I already tried eval, evalf, simplfy and it doesn't work.

Thanks a lot.

Hello

I have two vectors of angle Theta and Phi  given by

Theta=[ 0, Pi/3, 2*Pi/3, Pi];

Phi=[ 0, Pi/3 , 2*Pi/3 , Pi,4*Pi/3, 5*Pi/3, 2*Pi, 7*Pi/3];

I would like a general procedure how can I put on the unit sphere only the point 

x=cos(Phi).*sin(Theta);
 y=sin(Theta).*sin(Phi);
 z=cos(Theta);

where

(Theta, Phi) in the set {(0,Pi/3),  (0,Pi),  , (0,5*Pi/3),  

(Pi/3,0),   (Pi/3,2*Pi/3),   (Pi/3,4*Pi/3),   (Pi/3,2*Pi), 

(2*Pi/3,Pi/3),   (2*Pi/3,Pi),   (2*Pi/3,5*Pi/3),   (2*Pi/3,7*Pi/3), 

(PI,2*Pi/3),   (Pi,4*Pi/3),   (Pi,2*Pi) }

can we extract these set of point fro the definition of the two vectors Phi and Theta and then make the plot

of course we use 
Many thanks

Can somebody suggest an efficient way to specify / input a large list of binary variables in the LPSolve command, for example:

> LPSolve(objective function, constraints, binaryvariables={ x[0,0,1], .x[0,0,2], ...., x[i,j,k]}

Is it possible to assign a name to the set rather than input each element, x[i,j,k], manually?

Thanks in advance!

HI MaplePrimes.com and other watchers,

Please enjoy the attaced files about combinatorics.
You may already know what '4 choose 3' is.

an_excercise_in_combinatorics.mw

an_excercise_in_combinatorics.pdf

Hopefully this can be useful to the casual mathematical observer.

Regards,

Matt