Suppose I have a excel file with say 20 columns of data with equal number of rows of data. All numeric data

How to read each column data to a seperate variable

Or 

The entire data into a matrix 

The graph we draw in the graph theory package appear to be vectorial. So the export pdf is also clear. 

with(GraphTheory):
DrawGraph(PathGraph(2), stylesheet=[vertexborder=false,
vertexpadding=5,edgecolor = blue,vertexcolor=black,edgethickness=2],size=[250,250])

draw1.pdf

We know that DrawGraph only offers straight forms for drawn edges.  Sometimes we need to draw graphs with curves.  But when I try to use maple, I find that the vertices (by pointplot) are blurred after zooming in (or even without zooming out).

with(plottools):    
with(plots):
Bezier:=proc(x0,y0,x1,y1,x2,y2,x3,y3)
local f,g,c,l1,l2,l3,p0,p1,p2,p3;
f:=t->x0*(1-t)^3+3*x1*t*(1-t)^2+3*x2*t^2*(1-t)+x3*t^3:
g:=t->y0*(1-t)^3+3*y1*t*(1-t)^2+3*y2*t^2*(1-t)+y3*t^3:
c:=plot([f(t),g(t),t=0..1],thickness=2,scaling=constrained, color=blue,axes=none):
p0:=pointplot([x0,y0],symbol= solidcircle,symbolsize = 50,color=black):
p1:=pointplot([x3,y3],symbol= solidcircle,symbolsize = 50,color=black):
display(c,p0,p1,size=[300,300]);
end:

x0:=0:   y0:=0:
x1:=-0.2:  y1:=0.2:
x2:=0.5:  y2:=0.2:
x3:=0.3:  y3:=0:
p:=Bezier(x0,y0,x1,y1,x2,y2,x3,y3)

 draw2.pdf 

How to make these vertices vectorially drawn (like DrawGraph)？

If I have a set S:= [a,b,c],

I want the loop to give me: [[a,b,c],[a+b,c],[a+c,b],[a,b+c],[a+b+c]]

Consider the following matrix

Q := <cos(t), sin(t)>; sin(t), -cos(t)>

If I multiply the first column by the second column (at first glance this would see to be a 2x1 column multiplied by a 2x1 column, which isn't possible), I still get a result. 

(Q[1 .. 2, 1]) . (Q[1 .. 2, 2])

Now this isn't really that big a deal. Apparently, something like Q[1 .. 2, 1] is actually a Vector, so we're actually computing the dot product between two vectors. However, the result is the following

What is the bar over the variable t?

Hello

I wonder if there is a way to parallelize the function listed below.   It deals with a set of sets (dynamical systems).  As far as I can see the problem is with subs and simplify.  In some cases a division by zero shows up (the reason I added try..catch and had to use for). 

The procedure should remove the systems with imaginary coefficients and return the valid indices.  

Many thanks

Ed

Obs.:  The number of dynamical systems can easily be greater than 10,000.  Any speed improvement is welcome.

Example of one such a dynamical system is

[(-x*rho*xi[6]*(beta - 1)*sqrt(((rho + beta - 1)*sigma - beta + 1)^2*xi[8]^2/(beta - 1)^2) - xi[8]*(-2*y*sigma*xi[7]*xi[6]^2*(beta - 1)^2 + (sigma*rho + (sigma + 1)*(beta - 1))*x*rho*xi[6] + 2*sigma*z*rho^2*xi[8]))/(2*xi[6]*xi[8]*rho*(beta - 1)), (-(beta - 1)*(-y*sigma*xi[7]*xi[6]^2*(beta - 1)^2 + rho^2*(sigma*z*xi[8] + x^2))*xi[6]*sqrt(((rho + beta - 1)*sigma - beta + 1)^2*xi[8]^2/(beta - 1)^2) - (2*sigma*(beta - 1)^2*(y*xi[7]*(rho + beta - 1)*sigma/2 + x^2*rho + y*(beta - 1)*xi[7]/2)*xi[6]^3 + 2*x*z*rho*sigma*xi[8]*(beta - 1)^2*xi[6]^2 + (-z*xi[8]*(rho + beta - 1)*sigma^2 + ((2*beta^2*z - 3*beta*z + z)*xi[8] + x^2*rho - (beta - 1)*(-2*beta*y*xi[7] + x^2))*sigma + x^2*(beta - 1))*rho^2*xi[6] + 2*x*z*rho^3*sigma*xi[8])*xi[8])/(2*xi[8]*sigma*(beta - 1)*(rho^2 + (beta - 1)^2*xi[6]^2)*xi[6]*xi[7]), (-(-x*sigma*(beta - 1)^2*xi[6]^3 + (beta - 1)^2*(sigma*y*xi[7] + x^2)*xi[6]^2 - rho^2*sigma*x*xi[6] - sigma*z*rho^2*xi[8])*(beta - 1)*rho*sqrt(((rho + beta - 1)*sigma - beta + 1)^2*xi[8]^2/(beta - 1)^2) + xi[8]*(-2*y*sigma^2*xi[7]*(beta - 1)^4*xi[6]^4 + 2*(rho*(rho + beta - 1)*sigma/2 + (beta - 1)*((-beta + 1/2)*rho + y*(beta - 1)*xi[7]))*sigma*(beta - 1)^2*x*xi[6]^3 + (beta - 1)^2*(((-y*xi[7] + 2*z*xi[8])*rho + y*(beta - 1)*xi[7])*sigma^2 + (x^2*rho + (beta - 1)*(-2*beta*y*xi[7] - 2*beta*z*xi[8] + x^2 + y*xi[7]))*sigma - x^2*(beta - 1))*rho*xi[6]^2 + 2*(rho*(rho + beta - 1)*sigma/2 + (beta - 1)*((-beta + 1/2)*rho + y*(beta - 1)*xi[7]))*sigma*x*rho^2*xi[6] + xi[8]*sigma*z*((rho - beta + 1)*sigma - beta + 1)*rho^3))/(2*xi[8]^2*sigma*(beta - 1)*rho*(rho^2 + (beta - 1)^2*xi[6]^2))]

coefmodel

{beta = 8/3, rho = 28, sigma = 10}

and phi and coe equal to 

[[], [], []]

cleanSystems := proc(sys::list,phi::list,coe::list,coefmodel::{list,set})
description "This function applies coefficient values to a system and remove it if the coef is somehow imaginary":
local i,aaa:=table(),bbb:=table(),ccc:=table(),ind:=table():
for i from 1 to nops(sys) do
  try
   aaa[i]:=simplify(subs(coefmodel,sys[i])):
   if has(aaa[i],I) then
      aaa[i]:=NULL:
      bbb[i]:=NULL:
      ccc[i]:=NULL:
      ind[i]:=NULL:
   else
      bbb[i]:=simplify(subs(coefmodel,phi[i])):
      ccc[i]:=simplify(subs(coefmodel,coe[i])):
      ind[i]:=i:
   end if:
  catch:
   aaa[i]:=NULL:
   bbb[i]:=NULL:
   ccc[i]:=NULL:
   ind[i]:=NULL:
  end try:
end do:
aaa:=convert(aaa,list):print(nops(aaa)):
bbb:=convert(bbb,list):print(nops(bbb)):
ccc:=convert(ccc,list):print(nops(ccc)):
ind:=convert(ind,list):print(nops(ind)):
return(aaa,bbb,ccc,ind):
end proc:

How to plot Nabla Discrete Mittag - Leffler Function in Maple?

How can we unveil all the possible transformations that reduce the number of independent variables of a given pde? I tried it by using "InvarientSolutions" which gives eight possible transformations. But under those transformations, I got only two different odes. Is there any other method/command from which we can get other possible odes of a given pde?

pde_ode.mw

The transformed boundary conditions are

Please and please, I have been having this challenge for some time now and I would be so so happy if I can get a solution to my challenge. I have been worried on:

how do I use these definitions (Volterra Integral Equation and Caputo definition) to solve

Helper, please. You can use any definition you are familiar with also. Only I need a fractional definition.

Thanks

My question is about the physics package: when transforming the metric tensor to new coordinates, sometimes the new coordinate can be expressed easily as a variable inside the new matrix entries,  but the new matrix is still displayed completely in terms of the old variables. How can I force it to display everything in terms of the new ones?

I have in internal procedure ppp that has an input prjpsn:=3. This set a Vector positions of x,y,z to either if prjpsn=1 

zpsn=1,   xpsn=2,  ypsn=3

or if prjpsn =3   

xpsn=1,   ypsn=2, zpsn=3

This then are used in another procedure. I can't get the values to transfere.  prjpsn would default to 3 in the module but can be set externally to 1 if needed.
The relevant parts are highlighted in green text. Cant get the worksheet to display.

Maple Worksheet - Errorprjpsn:=3


Failed to load the worksheet /maplenet/convert/Q_14-12-22_Module_internal_settings.mw .

Download Q_14-12-22_Module_internal_settings.mw

How do I calculate the conformal killing vectors of a metric tensor in GRTensorIII?

Why I get "Error, invalid subscript selector" error in my code?

y1 := Grid:-Seq(UP1(s, U, V, W, Phi, Xi, N, a, b, II, JJ, A, B, Dd, M, Ns), s = 1 .. 7);
UKt := add(y1[i], i = 1 .. 7);

Error, invalid subscript selector

As you can see y1 is defined without any problem and have 7 seqments, but the next line warns invalid subscript selector.

I have, for example,

> x(t) := sin(a)*cos(t) + cos(a)*sin(t)
            x := t -> sin(a) cos(t) + cos(a) sin(t)
> y(t) := combine(x(t))
            y := t -> combine(x(t))

>y(t);
            sin(a + t)

> x(t);
            sin(a) cos(t) + cos(a) sin(t)

# wanted but cannot do this: x(t) := combine(x(t))

Any workaround?

Thanks.

why does the command Matrix(2, 2, undefined) create the following

Matrix(2, 2, [[undefined(1, 1), undefined(1, 2)], [undefined(2, 1), undefined(2, 2)]])

instead of simple undefined for each item?