Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Dear All,

I want all calculations in Maple code to be done with double precision. Using Digits:=18 and evalf greatly affects the calculation speed. But using the evalhf performs calculations much faster. The problem is that the evalhf cannot work with expressions that contain symbolic terms. For example, the evalhf(1/3*x+1/5*exp(-x^2)) generates an error. How can I perform calculations on expression with double-precision floating point without significantly slowing down the calculation speed? Is there a method that can be used at the beginning of the code to use double precision (preferably hardware) in all computation steps without slowing down the computations, unlike software floating-point computations?
Best wishes
 

Digits:=18;
evalf(1/3.);
evalhf(1/3.);


str:=time[real]();
for i from 1 to 100000 do
evalf(sqrt(3)*i):
end do:
time[real]()-str;

 

str:=time[real]();
for i from 1 to 100000 do
evalhf(sqrt(3)*i):
end do:
time[real]()-str;

 

Dear All,

I have the following algebraic function. In vibrations, it is common to write the functions exp(alpha[i]*I*t) as separate terms of cos(alpha[i]*t) and sin(alpha[i]*t). While exp(-beta[i]*t) remains without converting to sinh and cosh. How can I find the solution given for T(t) as the sum of terms C1[i]*sin(alpha[i]*t)*exp(-beta[i]*t) and C2[i]*cos(alpha[i]*t)*exp(-beta[i]*t) with non-complex (real) coefficients C1[i] and C2[i].

Can anyone help me to achieve my goal in the following expression?

 

T(t):=(1.450761945*10^(-11) + (3.836655196*10^(-14))*I)*exp((-0.5000000000 + 222.6866468*I)*t) + (-3.770333746*10^(-9) + (2.179000257*10^(-6))*I)*exp((-0.5000000000 - 924.5904413*I)*t) + (-2.584086158*10^(-12) + (4.273321932*10^(-13))*I)*exp((-0.5000000000 + 326.7549627*I)*t) + (1.986287340*10^(-9) + (1.330623218*10^(-11))*I)*exp((-0.5000000000 - 74.63720909*I)*t) + (-5.980910367*10^(-12) + (5.816480027*10^(-11))*I)*exp((-0.5000000000 - 453.7574402*I)*t) + (1.450761945*10^(-11) - (3.836655196*10^(-14))*I)*exp((-0.5000000000 - 222.6866468*I)*t) + (8.923968224*10^(-10) - (8.844466162*10^(-9))*I)*exp((-0.5000000000 + 637.9999953*I)*t) - (1.217986141*10^(-10) + (4.431771836*10^(-13))*I)*exp((-0.5000000000 - 138.7904660*I)*t) + (-1.217986141*10^(-10) + (4.431771836*10^(-13))*I)*exp((-0.5000000000 + 138.7904660*I)*t) + (0.0002537882980 + 0.00002277791755*I)*exp((-0.5000000000 - 5.570928456*I)*t) - (3.770333746*10^(-9) + (2.179000257*10^(-6))*I)*exp((-0.5000000000 + 924.5904413*I)*t) + (-0.0001618723219 + 0.01288449595*I)*exp((-0.5000000000 - 1638.001654*I)*t) + (8.923968224*10^(-10) + (8.844466162*10^(-9))*I)*exp((-0.5000000000 - 637.9999953*I)*t) + (-1.153529195*10^(-7) + (1.908444485*10^(-9))*I)*exp((-0.5000000000 + 30.22171212*I)*t) - (1.153529195*10^(-7) + (1.908444485*10^(-9))*I)*exp((-0.5000000000 - 30.22171212*I)*t) - (0.0001618723219 + 0.01288449595*I)*exp((-0.5000000000 + 1638.001654*I)*t) + (1.986287340*10^(-9) - (1.330623218*10^(-11))*I)*exp((-0.5000000000 + 74.63720909*I)*t) + (0.0002537882980 - 0.00002277791755*I)*exp((-0.5000000000 + 5.570928456*I)*t) - (5.980910367*10^(-12) + (5.816480027*10^(-11))*I)*exp((-0.5000000000 + 453.7574402*I)*t) - (2.584086158*10^(-12) + (4.273321932*10^(-13))*I)*exp((-0.5000000000 - 326.7549627*I)*t);

Dear All,

I have the given expression f as follows. I want to extract terms with sin, sinh, and exp in this expression. For example, for sin, cos, sinh and exp, I want the outputs to be -4*sin(x), -5*cos(x^3)*sin(y^2), 5*sinh(x^2) and 2*exp(y^2), respectively. Also, if I want to get terms having sin(x) and sin(y^2), I like to receive -4*sin(x) and -5*cos(x^3)*sin(y^2), respectively. And also, for terms including sinh(x^3), I want to get Void output. The above goals cannot be achieved with the following commands. Can anyone help me?

f := -4*sin(x) + 2*exp(y^2) + 5 - 5*cos(x^3)*sin(y^2) + 5*sinh(x^2);

indets(f, trig);

selectfun(f, exp);

selectremove(has, f, trig);

Multiplaction "dot" in Maple 2022 is way too small - causes errors.

e.g. two variables multiplied s*m ends up being sm a new variable as I cannot really see that there is a missing multiplication operator between the variables. This causes huge unnecessary errors.

Maple 9.x e.g had nice clear and big operators and this kind of error was avoided.

How can I undo this unfortunate regression in Maple 2022 to increase the size of multiplication operator and other operators, so that they actually becom visible and not just a little dot almost a pixel in size.

If I was a falcon (20x20)^infinity then this would have been ok, but I am not, I am human.

So how do I change this unfortunate regression so that these errors can be avoided.?

It was found on the social networks of the WM group. Written in Python. Perhaps someone would like to adopt it.
 

Deal All,

I have a system of linear differential equations with unknown functions T[1](t) to T[n](t). In the attached example, I considered the value of n equal to 10, but depending on the problem, the value of n may be higher. Maple is not able to solve this problem analytically with the ‘dsolve’. 

Does anyone have an idea to analytically solve for such a set of linear differential equations?

Best wishes

Set_of_Linear_DEs.mw

See attached worksheet in Maple 2023.

This example is taken from the Maple help page. I want to 'zoom in' on a plot3d object. The only way I have found was from responses [1] on the maple primes forum. It uses InlinePlot and the scale option to perform the 'zoom in'. Since InlinePlot generates the plot in terms of XML there is no graphic out, only a text based output. In order to reconstitute the InlinePlot as a plot object I can view visually I need to use some additional commands from the DocumentTool package. This is all great but the output, which in our case is P3, is not a plot object and therefore cannot be exported as a png. Is there a way to convert the InlinePlot with the scaling applied back to a typical plot object so I can export it as a .png, using Export("output_plot.png",P3,base=worksheetdir)?

can_I_convert_InlinePlot(P3)_back_to_a_regular_plot_object_so_I_can_export_it_as_a_png.mw

When I try to contract the tensor with connection, maple encounters such promble:

Error, (in DifferentialGeometry:-Tensor:-ContractIndices) expected 2nd argument to be a tensor. Received: _DG([["connection", O, [["con_bas", "cov_bas", "cov_bas"], []]], [`...`]])

The expression is Cacd Hab

I try to solve triple integraton in Maple with this code.

r := a + (b - a)*z/h;
x1 := sqrt(r^2 - y^2);
V := int(int(int(1, x = -x1 .. x1), y = -r .. r), z = 0 .. h);

but it leaves the last integral dz in the answer and warns: unable to determine if a*h/(-b+a) is between 0 and h; try to use assumptions or use the AllSolutions option
What is the problem?
and i need to get V = Pi*h(a^2 + ab +b^2)/3

Regards

Could someone help me understand why Maple hangs solving these two equations when the names of unknowns is c__1,c__2   or the old _C1, _C2?

This worksheets shows this. I had to put timelimit of 90 seconds, else it hangs may be forever.    

This example came from looking at why Maple seems to hang randomly when I run the same problem. I still do not know why, but is seems to have something to do with the use of  lower case c__1 or _C1. 

16020

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1744 and is the same as the version installed in this computer, created 2024, April 17, 19:33 hours Pacific Time.`

restart;

dsolve(diff(y(x),x$9)=1,arbitraryconstants=subscripted);

y(x) = (1/362880)*x^9+(1/40320)*c__1*x^8+(1/5040)*c__2*x^7+(1/720)*c__3*x^6+(1/120)*c__4*x^5+(1/24)*c__5*x^4+(1/6)*c__6*x^3+(1/2)*c__7*x^2+c__8*x+c__9

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*c__1)*sqrt(2)/sqrt(c__1))*sqrt(2)/sqrt(c__1) = c__2, 1 = (-exp(2*c__2*sqrt(c__1)*sqrt(2)) + 1)*sqrt(c__1)*sqrt(2)/(exp(2*c__2*sqrt(c__1)*sqrt(2)) + 1)];
unknowns:=[c__1,c__2];
time();
C_sol:=[timelimit(90,solve(eqs,unknowns))];
time()

[-(1/2)*arctanh((1/2)*(1+2*c__1)^(1/2)*2^(1/2)/c__1^(1/2))*2^(1/2)/c__1^(1/2) = c__2, 1 = (-exp(2*c__2*c__1^(1/2)*2^(1/2))+1)*c__1^(1/2)*2^(1/2)/(exp(2*c__2*c__1^(1/2)*2^(1/2))+1)]

[c__1, c__2]

.609

Error, (in evalf/hypergeom) time expired

106.125

restart; #try again but do not do the arbitraryconstants=subscripted now.

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*c__1)*sqrt(2)/sqrt(c__1))*sqrt(2)/sqrt(c__1) = c__2, 1 = (-exp(2*c__2*sqrt(c__1)*sqrt(2)) + 1)*sqrt(c__1)*sqrt(2)/(exp(2*c__2*sqrt(c__1)*sqrt(2)) + 1)];
unknowns:=[c__1,c__2];
time();
C_sol:=[timelimit(90,solve(eqs,unknowns))];
time()

[-(1/2)*arctanh((1/2)*(1+2*c__1)^(1/2)*2^(1/2)/c__1^(1/2))*2^(1/2)/c__1^(1/2) = c__2, 1 = (-exp(2*c__2*c__1^(1/2)*2^(1/2))+1)*c__1^(1/2)*2^(1/2)/(exp(2*c__2*c__1^(1/2)*2^(1/2))+1)]

[c__1, c__2]

106.234

[[]]

106.390

#see? it finisghed instantly now.


Download why_solve_hangs_with_subscripted_may_15_2024.mw

The same thing happens If I use the old _C1 and _C2 instead of c__1 and c__2. It also hangs. 

The following worksheet shows this.  If I change _C1 and _C2 to other symbols, say C1 and C2, then it does not hang. 

Why the names of the unknowns makes difference to solve?
 

18792

restart;

18792

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*_C1)*sqrt(2)/sqrt(_C1))*sqrt(2)/sqrt(_C1) = _C2, 1 = (-exp(2*_C2*sqrt(_C1)*sqrt(2)) + 1)*sqrt(_C1)*sqrt(2)/(exp(2*_C2*sqrt(_C1)*sqrt(2)) + 1)];
unknowns:=[_C1,_C2];
time();
C_sol:=[timelimit(90,solve(eqs,unknowns))];
time()

[-(1/2)*arctanh((1/2)*(1+2*_C1)^(1/2)*2^(1/2)/_C1^(1/2))*2^(1/2)/_C1^(1/2) = _C2, 1 = (-exp(2*_C2*_C1^(1/2)*2^(1/2))+1)*_C1^(1/2)*2^(1/2)/(exp(2*_C2*_C1^(1/2)*2^(1/2))+1)]

[_C1, _C2]

.125

Error, (in evalf/cos) time expired

105.218

restart;

18792

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*C1)*sqrt(2)/sqrt(C1))*sqrt(2)/sqrt(C1) = C2, 1 = (-exp(2*C2*sqrt(C1)*sqrt(2)) + 1)*sqrt(C1)*sqrt(2)/(exp(2*C2*sqrt(C1)*sqrt(2)) + 1)];
unknowns:=[C1,C2];
time();
C_sol:=[timelimit(90,solve(eqs,unknowns))];
time()

[-(1/2)*arctanh((1/2)*(1+2*C1)^(1/2)*2^(1/2)/C1^(1/2))*2^(1/2)/C1^(1/2) = C2, 1 = (-exp(2*C2*C1^(1/2)*2^(1/2))+1)*C1^(1/2)*2^(1/2)/(exp(2*C2*C1^(1/2)*2^(1/2))+1)]

[C1, C2]

105.312

[[]]

105.468

 

 

Download why_solve_hangs_with_OLD_C_also_may_15_2024.mw

Only thing I see in help related to name of symbols to solve for is this:

The solve command solves one or more equations or inequalities for the specified unknowns. The unknowns may be names, including indexed names (though for efficiency reasons, indexed names should be avoided when possible), 

ps. This looks like a bug to me. So I send bug report to Maplesoft support also.

 

pps. I tried this in Maple 2023 and Maple 2022 and same behavior. Could someone with earlier version of Maple try to see if this behavior was there also?   It can possibly be correct that the choice of letter used makes difference for solving equations. I have to use c__1 and c__2 etc.. since these equations come from differential equations and this is what I use for constants of integrations., Otherwise I have to make lots of changes now to use different letters.

Update

did trace on solve using _C1 and _C2 and then using A,B for variables to solve for. Code flow is different. This expalins why it hangs. The flow starts the same until it gets to 

TriangularDecomposition: something went wrong during backsubstitution - trying a different variable order

Then when using _C1 and _C2 the code goes into different path than when using A,B. 

Here is flow when using _C1,_C2.

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*_C1)*sqrt(2)/sqrt(_C1))*sqrt(2)/sqrt(_C1) = _C2, 1 = (-exp(2*_C2*sqrt(_C1)*sqrt(2)) + 1)*sqrt(_C1)*sqrt(2)/(exp(2*_C2*sqrt(_C1)*sqrt(2)) + 1)];
unknowns:={_C1,_C2};
time();
infolevel[solve]:=5;
C_sol:=[timelimit(60,solve(eqs,unknowns))]; 
time()


Main: Entering solver with 2 equations in 2 variables
Main: attempting to solve as a linear system
Dispatch: dispatching to Radicals handler
Recurse: recursively solving 2 equations and 2 inequations in 2 variables
Dispatch: dispatching to Radicals handler
Transformer:   solving for linear equation in _X000001
Recurse: recursively solving 2 equations and 2 inequations in 2 variables
Dispatch: dispatching to Radicals handler
Recurse: recursively solving 3 equations and 2 inequations in 3 variables
Dispatch: dispatching to Exponentials handler
Transformer:   solving for linear equation in _S000004
Recurse: recursively solving 3 equations and 3 inequations in 4 variables
Dispatch: dispatching to Rename handler
Dispatch: renaming _S000005 = arctanh(1/2*_S000002*RootOf(_Z^2-2,index = 1)/_S000001)
Recurse: recursively solving 3 equations and 3 inequations in 5 variables
Dispatch: handling polynomials of the form a*x^n-b
Dispatch: dispatching to PolynomialSystem handler
Main: polynomial system split into 1 parts under preprocessing
Main: using RegularChains based methods
SolverVariableOrder: using the variable order  _S000003 > _X000002 > _S000005 > _S000002 > _S000001
TriangularDecomposition: using deterministic algorithm for decomposition
TriangularDecomposition: decomposition successfully found 1 components
TriangularDecomposition: backsubstituting to form solutions
TriangularDecomposition: something went wrong during backsubstitution - trying a different variable order
Transformer:   solving the uncoupled linear subsystem in {_S000003, _X000002}
Linear: solving 2 linear equations
Polynomial: # of equations is: 2
Polynomial: best equation / unknown _S000005*_z1 _X000002 2*_S000001
Polynomial: # of equations is: 1
Polynomial: best equation / unknown -_z1*_S000001+1 _S000003 _z1*_S000001+1
Polynomial: backsubstitution at: 2
Polynomial: backsubstitution at: 1
Main: polynomial system split into 1 parts under preprocessing
Main: subsystem is essentially univariate
UnivariateHandler: subsystem has only one equation
UnivariateHandler: solving as if univariate in _S000002
Recurse: recursively solving 1 equations and 0 inequations in 1 variables
Dispatch: dispatching to OnlyIn handler
Transformer:   solving for linear equation in _S000006
Recurse: recursively solving 1 equations and 0 inequations in 1 variables
Transformer:   solving the uncoupled linear subsystem in _S000006
Recurse: recursively solving 1 equations and 1 inequations in 1 variables
Transformer:   solving the uncoupled linear subsystem in t
Error, (in evalf/hypergeom) time expired

 

This is trace when using A,B

 

eqs:=[-1/2*arctanh(1/2*sqrt(1 + 2*A)*sqrt(2)/sqrt(A))*sqrt(2)/sqrt(A) =B, 1 = (-exp(2*B*sqrt(A)*sqrt(2)) + 1)*sqrt(A)*sqrt(2)/(exp(2*B*sqrt(A)*sqrt(2)) + 1)];
unknowns:={A,B};
time();
infolevel[solve]:=5;
C_sol:=[timelimit(90,solve(eqs,unknowns))]; 
time()


Main: Entering solver with 2 equations in 2 variables
Main: attempting to solve as a linear system
Dispatch: dispatching to Radicals handler
Recurse: recursively solving 2 equations and 2 inequations in 2 variables
Dispatch: dispatching to Radicals handler
Transformer:   solving for linear equation in A
Recurse: recursively solving 2 equations and 2 inequations in 2 variables
Dispatch: dispatching to Radicals handler
Recurse: recursively solving 3 equations and 2 inequations in 3 variables
Dispatch: dispatching to Exponentials handler
Transformer:   solving for linear equation in _S000004
Recurse: recursively solving 3 equations and 3 inequations in 4 variables
Dispatch: dispatching to Rename handler
Dispatch: renaming _S000005 = arctanh(1/2*_S000002*RootOf(_Z^2-2,index = 1)/_S000001)
Recurse: recursively solving 3 equations and 3 inequations in 5 variables
Dispatch: handling polynomials of the form a*x^n-b
Dispatch: dispatching to PolynomialSystem handler
Main: polynomial system split into 1 parts under preprocessing
Main: using RegularChains based methods
SolverVariableOrder: using the variable order  _S000003 > B > _S000005 > _S000002 > _S000001
TriangularDecomposition: using deterministic algorithm for decomposition
TriangularDecomposition: decomposition successfully found 1 components
TriangularDecomposition: backsubstituting to form solutions
TriangularDecomposition: something went wrong during backsubstitution - trying a different variable order
Main: polynomial system split into 1 parts under preprocessing
Main: applying the solver for domain=absolute, engine=traditional
PseudoResultant: 225530 [1 200002087 _S000001] 3 3 151 2 45 0
PseudoResultant: 120059 [2 200004772 _S000003] 2 2 131 0 3 0
PseudoResultant: 134507 [1 700002396 _S000002] 1 1 53 0 3 0
PseudoResultant: -10 [] 0 0 3 0 3 0
PseudoResultant: 1 solutions found, now doing backsubstitution
PseudoResultant: backsubstitution of _S000002
PseudoResultant: backsubstitution of _S000003
PseudoResultant: backsubstitution of _S000001
SolutionsLost: setting solutions lost flag
Main: solving successful - now forming solutions
Main: Exiting solver returning 0 solutions

solve: Warning: no solutions found

 

Why does Maple 2022 still have only a history of ten files in the Recents ?

It makes it very difficult to search for files you opened but only realised later you need, but then it is long out of the silly 10 history list.

Isnt there a way to make a perpetual list so that all files opened are saved chronologiacally against date, or... at least be able to increase the silly 10 recent files to something like 30

Thanks

This is a follow on from and earlier question on tabulation. Some procdures returs up tp eight values. So using a Table make it easier to see what is what. Sometimes the equations returned are very long. Is there a way to truncate what is display in the table using something along the lines of term elision? Say the 1st 300 chatacters for example?  

Edit: Corrected an error in the worksheet.

restart

QQFProj := proc(q12::algebraic, q23::algebraic,
                q34::algebraic, q14::algebraic,
                  prnt::boolean:=true) #{columns:=[QQFproj,Q13proj,Q24proj]}
  description "Projective quadruple quad formula and intermediate 13 and 24 quads. Useful for cyclic quadrilaterals";
  local qqf,q13,q24, sub1,sub2,sub3, R,values,DF,lens;
  uses   DocumentTools;
  sub1:= (q12 + q23 + q34 + q14)^2 - 2*(q12^2 + q23^2 + q34^2 + q14^2) ;
  sub2:=-4*(q12*q23*q34+q12*q23*q14+q12*q34*q14+q23*q34*q14)+8*q12*q23*q34*q14;
  sub3:=64*q12*q23*q34*q14*(1-q12)*(1-q23)*(1-q34)*(1-q14);
  qqf:=((sub1+sub2)^2=sub3);
  q13:=((q12-q23)^2-(q34-q14)^2)/(2*(q12+q23-q34-q14-2*q12*q23+2*q34*q14));#check this
  q24:=((q23-q34)^2-(q12-q14)^2)/(2*(q23+q34-q12-q14-2*q23*q34+2*q12*q14));#check this
  if prnt then
  
   values:=<qqf,q13,q24>;
   DF:=DataFrame(<values>, columns=[`"Values Equations"`],rows=[`#1  QQF`,`#2  Q13`,`#3  Q24`]);
   lens := [4 +8* max(op(length~(RowLabels(DF)))),4+ min(max( 10*(length~(values))),1000)];#op(length~(ColumnLabels(DF)0)
   Tabulate(DF,width=add(lens),widthmode = pixels,weights = lens);
  return qqf,q13,q24
  end if;
  return qqf,q13,q24
end proc:

 q12:=1/2:q23:=9/10:q34:=25/26:q41:=9/130:#Cyclic quadrilateral
q12:=sqrt(17+a)/2:q23:=expand(r^2+t^2)^2/10:q34:=expand((a+b+c)^4/26):q41:=sqrt(17+b)/130:

Q:=QQFProj(q12,q23,q34,q41,true):

#Q[1]

#(Q[2])

#Q[3]

length(Q[1])

2669

(1)

length(Q[2])

1024

(2)

length(Q[3])

1024

(3)
 

 

Download 2024-05-14_Q_Data_Table_Limit_Equation_lentgh.mw

In my code, why  does InitState(5,5) return only one element of the vector while Coherent(5) returns 5 elements of the vector?

restart;

with(LinearAlgebra):

 

Coherent := proc({zeta:=1,phi:=Pi/2},n_max) local alpha,i,ICs:
 #alpha := sqrt(zeta)*exp(I*phi):
 ICs := Vector(n_max,i->evalc(evalf(exp(-zeta/2)/sqrt((i-1)!))*(sqrt(zeta)*exp(I*phi))^(i-1)),datatype=complex[8]);
end proc:

Projectile := proc({L:=1,sigma:=0.01,beta:=0.02,k0:=-5*Pi},n_max) local x0,g,c,v,ICs,j:
  x0 := evalf(beta*L);
  g := unapply(exp(-(x-x0)^2/2/sigma^2),x);
  c := evalf(int(g(x)^2,x=-L/2..L/2,numeric=true));
 
  v:=Vector(n_max,datatype=complex[8]);

  for j from 1 to n_max do:
     if (is(j,odd)) then  v[j]:= evalc[8](evalf(Int(cos(Pi*j*x/L)*g(x)*cos(k0*x),x=-L/2..L/2,method=_d01akc))+ I*evalf(Int(cos(Pi*j*x/L)*g(x)*sin(k0*x),x=-L/2..L/2,method=_d01akc)));
      else v[j]:= evalc[8](evalf(Int(sin(Pi*j*x/L)*g(x)*cos(k0*x),x=-L/2..L/2,method=_d01akc))+ I*evalf(Int(sin(Pi*j*x/L)*g(x)*sin(k0*x),x=-L/2..L/2,method=_d01akc)));
     end if:
  end do:

 ICs :=evalf[8](sqrt(2/L*c))*v;
end proc:

InitState := proc({zeta:=1,phi:=Pi/2,L:=1,sigma:=0.01,beta:=0.02,k0:=-5*Pi},d1,d2) local Z0:
z1 := Coherent(zeta,phi,d1);
#Z0:= MatrixMatrixMultiply(Coherent(zeta,phi,d1),Transpose(Projectile(L,sigma,beta,k0,d2))):
end proc:

Warning, (in InitState) `z1` is implicitly declared local

 

InitState(5,5); Coherent(5);

Vector(1, {(1) = .6065306597+0.*I})

 

Vector[column](%id = 36893490076814315516)

(1)
 

 

Download test3.mw

Can I export a Table as an image like .png, where Table is defined (with DocumentTools)? See Maple worksheet for example.

Why would anyone want to do this? It all started because I wanted to include a color bar(with a specific color range) in my 3dplot. There is no native way to do this with plot3d so I searched Maple Primes for alternative strategies. One strategy is to generate the 3d plot and the color bar(with plot3d) seperately, then combine the plots in a table so they sit side by side, using with(DocumentTools). I have been almost successful with this strategy. There remain two outstanding problems. 1. I can't re-size the table cells since there are no such options with(DocumentTools). The color bar should have a smaller cell because the figure itself is tall and thin. 2. I need to export this combined object as an image(.png) but its a table with plot objects inside, and not itself a plot object and therefore one can't simply export it as a .png like one would a typical plot object. Is there a way to export this table as a png? I am beginning to think that my idea of combining  plots with document tools and attempting the export the resulting table is not feasible. How does a normal person add a custom-color-range color bar to their plot? I'm not trying to move the earth here but it certainly feels like it.

with_document_tools_how_to_i_resize_the_cells_and_export_the_table_as_png.mw

 I can't achieve this using StringTools:-RegSplit. The StringTools,Regular_Expressions documentation doesn't seem to address these functionalities.

My goal is to split the string at the occurences of "." that are somewhere between "]" and "["* yet not enclosed between the two digits, and to ignore all other occurences of "."

1 2 3 4 5 6 7 Last Page 2 of 2112