Maple 2024 Questions and Posts

What "style" are the icons?...

Asked by: Scot Gould 1059

December 01 2024

1 0

Question: Is there a style for the icons from the palettes?

(The rest of this post is a screenshot because MaplePrimes can't show the change in fonts.)s a

DataFrame display alpha as a Greek letter...

Asked by: Ronan 1381

December 01 2024

1 4

How do I get alpha to display as the Greek letter in DataFrame Row column?

>

restart

>	#varp:=alpha

>

QQFProj := proc(q12::algebraic, q23::algebraic,
                q34::algebraic, q14::algebraic,{varp:=:-alpha},
                  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);
  sub2:=-4*(q12*q23*q34);
  sub3:=64*q12*q23;
  qqf:=(sub1+sub2)^2-sub3;
  q13:=(q12-q23)^2;
  q24:=varp*(q23-q34)^2;
  if prnt then

   values:=<qqf,q13,q24>;
   DF:=DataFrame(<values>, columns=[`"Values Equations"`],rows=[`#1 QQF`,`#2 Q13`,cat(`#3 Q24 (`,varp,`)`)]);
   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:=(r^2+t^2)^2/10:q34:=((a+b+c)^4/26):q41:=sqrt(17+b)/130:

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

Download 2024-12-01_Q_Data_Table_alpha_as_Greek_Letter.mw

A Christmas package needs to be wrapped......

Asked by: Alfred_F 440

December 01 2024

1 2

An oddly shaped Christmas package is a collection of centrally stacked cube-shaped boxes with ever smaller sides. From bottom to top, the edge lengths of the stacked cubes run like the sequence of the reciprocals of the sequence of the square roots of the natural numbers. How high is the package? How much paper is needed to wrap the package? What is the total volume enclosed by all the packages?

Why when i take integral from whole equation is n...

Asked by: salim-barzani 1695

November 30 2024

0 1

why i take intgral from equation it must apply for each term but just by sign is show up how it will work?

>

(1)

>

(2)

>

(3)

>

(4)

>

w(x, t) = 2*(diff(diff(f(x, t), x), x))/f(x, t)-2*(diff(f(x, t), x))^2/f(x, t)^2

(5)

>

(6)

>

(7)

>

(8)

(9)

Download int.mw

Unit conversion of results in the form of equations...

by: C_R 3562 Product: Maple 2024

November 29 2024

1 0

As an enthusiast and frequent user of unit functionalities, I would like to suggest some improvements.

Example 1:

Assume a quantity equation where we substitute (and probably other parameters with units) to obtain a numerical result.
After some simplifications the result is

(1)

A result in this form is the most frequent physical quantity equation with units in use:

A physical quantity denoted by the name that is equal to a numerical value times a physical unit.

In this form it is also the shortest way to define/name scientific constants, scientific and engineering parameters, and results for variables/unknowns.

To better interpret the above result a conversion to micron is desirable. This unfortunately does not work with the context panel. Alternatively one could think of mapping a command onto the equation.
(Mapping a command is a good practice to threat the left-hand side and the right-hand side of an equation equally.)

Error, (in convert/units) unable to convert `1` to `micron`

This however does not work, since the left-hand side does not contain the dimension length, but only the name of the physical quantity.

A way to convert is:

(2)

For a single result this might be acceptable, but for a list of results this not practical.

In any case not being capable to apply to an equation the same operation to both sides bares the risks of introducing errors.

Example 2:

More complicated are results of this kind

(3)

A way to convert this to :

(4)

This is too complicated.

Instead, eval with a manually entered conversion can be used for this

(5)

However, this is old style: The way we had to work before the introduction of the unit package.

In other areas Maples unit conversion functionalities are more advanced.

Plotting with units for example has been improved to the point where unit conversion errors of quantity expressions are almost impossible (in contrast to numerical-value expressions): The user can specify the units to display and the numerical values of the physical quantities are converted accordingly.

S1: To catch up with that improvement, conversion of the results as in Example 1 in the form of quantity equations would be desirable as well.

Other desirable unit conversions of quantity expressions are

•	S2: Right click (or select and right click) a unit in Maple output and conversion to a new unit with the context panel (try it with one of the above outputs: nothing happens).

•	S3: Replacement of a unit by a desired unit of the same dimension (Example 2).

•	S4: Replacement of a unit by conversion to a new equivalent unit (e.g.: ) matching the pattern of the unit but not necessarily the pattern of the underlying Unit commands (for this subs can be used already).

For all suggestions, it is important that the physical quantity to be converted does not change, i.e. if a unit is changed, the numerical value of the physical quantity is changed accordingly.

Download Unit_conversion_in_quantity_expressions.mw

I don't like indenting input text...

Asked by: michele 90

November 29 2024

2 11

With the latest version of Maple I get an indenting space when I click "Shift+Enter".
For example

But I don't like it: suggestion? Thanks

How to get this recursive sequence of A[n] (warnin...

Asked by: Aixleft math 120

November 26 2024

1 4

Hi, how to rsolve this recursive sequence auto?

The manual calculation is taking the reciprocal of both sides of the equation, and then do some simple calculations to get the answer:1/(3·2^n-5).

Download ask_recrusive_sequence_of_A[n].mw

Why 2d input sometimes becomes roman...

Asked by: C_R 3562

November 25 2024

4 10

Something has changed in the past weeks (not sure for how long) and I have not figured out what exactly.

Probably it is me doing something differently while editing an input line. Below is an example where the name data is not italic any more

Anybody else experiencing the same?

My interface seetings:

Error in NLPSolve and recursion for parameter fitt...

Asked by: countryplaza 10

November 24 2024

2 4

Hello to whomever this may concern. I am curently trying to estimate parameters and initial conditions for a system of ODEs. I have taken inspiration from the code given by Robert J. Lopez. I am having issues with minimizing the SSE (sum of squared errors) which I named as a procedure using the NLPSolve command from the Optimization package. It returns an error saying there are too many levels of recursion, which I do not understand where it is coming from or how to solve the issue and I was hoping to get an explination on how to solve it.

Below is the code up until the error. Any help or guidance would be greatly appreciated.

Parameter_estimation.mw

problem in replacing mu in equation ?...

Asked by: salim-barzani 1695

November 24 2024

0 4

when i replace mu=0 in my equation mu[3] and mu[4] also become zero why ?

>

(1)

>

(2)

>

U1 := lambda*(1+sqrt(-lambda)*(cosh(xi*sqrt(-lambda))*B[2]+sinh(xi*sqrt(-lambda))*B[1])*`&+-`(sqrt(lambda/(lambda^2*B[1]^2-lambda^2*B[2]^2-mu^2))))*beta[0]*exp(I*(d*x+e*y+f*t^beta/beta))/(B[1]*cosh(xi*sqrt(-lambda))*lambda+B[2]*sinh(xi*sqrt(-lambda))*lambda+mu)

(1+(-lambda)^(1/2)*cosh(xi*(-lambda)^(1/2))*B[2]*`&+-`((-1/(lambda*B[2]^2))^(1/2)))*beta[0]*exp(I*(d*x+e*y+f*t^beta/beta))/(B[2]*sinh(xi*(-lambda)^(1/2)))

(3)

>

V1 := -mu[4]*b*lambda^2*(1+sqrt(-lambda)*(cosh(xi*sqrt(-lambda))*B[2]+sinh(xi*sqrt(-lambda))*B[1])*`&+-`(sqrt(lambda/(lambda^2*B[1]^2-lambda^2*B[2]^2-mu^2))))^2*beta[0]^2/((B[1]*cosh(xi*sqrt(-lambda))*lambda+B[2]*sinh(xi*sqrt(-lambda))*lambda+mu)^2*a*mu[3])

Error, numeric exception: division by zero

Download tester.mw

Find ODE solution with condition parameter ?...

Asked by: salim-barzani 1695

November 23 2024

0 12

How we can calculate solution of this ODE by give assumption to the equation , we have otehr case too, like lambda>0&mu<0

why is latex of Laplace is empty?...

Asked by: nm 11693

November 21 2024

1 3

Any idea why Maple returns empty string when asked for the latex of the Laplace of x(t)? Am I doing anything wrong here? I do not see it

Download latex_of_laplace_nov_21_2024.mw

I was expecting something like this using another software

Did Maple always behave this way for Laplace? I do not have earlier version now to check. Any workaround?

Error, (in PDEtools:-Solve) found functions with s...

Asked by: nm 11693

November 21 2024

1 1

Is this a known error when using PDEtools:-Solve? I get no error using solve on same input, so thought to ask, just in case it should not happen.

For now, I will change my code to use solve for this.

I am basically solving two equations in Laplace domain for Y1(s) and Y2(s). But since there are initial conditions x(0) and y(0) in the equations, and Laplace has L(x(t),t,s) then PDEtools:-Solve is not happy, as it sees x(t) and x(0) in same input.

But solve has no problem with this. Who is correct? solve or Solve?

>	interface(version);

>	Physics:-Version();

>

libname;

$"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"$

>

restart;

>	sys:=[slaplace(x(t),t,s)-x(0) = -3laplace(x(t),t,s)+4laplace(y(t),t,s)+1/(s^2+1), slaplace(y(t),t,s)-y(0) = -2laplace(x(t),t,s)+3laplace(y(t),t,s)+1/s^2]; Ys:=[laplace(x(t),t,s), laplace(y(t),t,s)];

[s*laplace(x(t), t, s)-x(0) = -3*laplace(x(t), t, s)+4*laplace(y(t), t, s)+1/(s^2+1), s*laplace(y(t), t, s)-y(0) = -2*laplace(x(t), t, s)+3*laplace(y(t), t, s)+1/s^2]

>	solve(sys,Ys);

[[laplace(x(t), t, s) = (x(0)*s^5+4*y(0)*s^4-3*x(0)*s^4+x(0)*s^3+4*y(0)*s^2-3*x(0)*s^2+s^3+s^2+4)/((s^4-1)*s^2), laplace(y(t), t, s) = (y(0)*s^5+3*y(0)*s^4-2*x(0)*s^4+y(0)*s^3+3*y(0)*s^2-2*x(0)*s^2+s^3+s^2+s+3)/((s^4-1)*s^2)]]

>	PDEtools:-Solve(sys,Ys)

Error, (in PDEtools:-Solve) found functions with same name but depending on different arguments in the given DE system: x(0), x(t). Specification of the dependent variables is required

Download differenece_between_solve_and_Solve_nov_20_2024.mw

Screen display on Maple 2020 vs Maple 2024...

Asked by: Jean-Marc Roy 20

November 20 2024

1 1

Any reason why the display of a fairly large plot in Maple 2024 (only show a print screen here)

be worst then the display of the same plot in Maple 2020

In 2024, it seams to be rasterized, while in 2020 it is still in vector form. Also the plot does not resize as well in 2024 compare to 2020. Any hint would help!

Maybe large plot are displayed in raster image, there is probably a setting somewhere in the documentation. When I export both are in vector format...

Thanks!

Experimental format for projective vectors...

Asked by: Ronan 1381

November 20 2024

1 8

I am experimenting using the this format of Vector( [Vector] ) to make projective vectors a different data type to Vectors. I don't want to use 1 x 3 or 3 x 1 matrices. The format holds some promise.
I would like to be able to copy the Maple format of Vector or Vector[column] and Vector[row] for my varaition.

ProjVectoC and ProjVectorR so ProjVector or ProjVector[column] and ProjVector[row]
A secondary question is on type checking (see previous question How to setup special type check in a procedure? - MaplePrimes ). Would it be possible to have the type check return ProjVector[column] or ProjVector[row]?
The attached worksheet contains a procedure for factor reducing the vectors to to a minimal format of <x,y,z>. Also Cross product and Dot product procedures to suit.

I am open to any efficiency improvements.

>

restart

>	interface(rtablesize=50)

(1)

>	with(LinearAlgebra):

>

FactReduce:=overload([
     proc(v::{list,Vector})
          option overload;
          description " removes linear factor from",
                      " a list, vector, matrix or expression";
          uses LinearAlgebra;
          local i, num,tgdc,dnm, V1;
          num:=`ifelse`(type(v,Vector),numelems(v),nops(v));
          dnm:=frontend(lcm, [seq(denom(v[i]),i=1..num)]);
          V1:=radnormal(v*~dnm);
          tgdc:=V1[1];

          for i from 2 to num do
               tgdc:=frontend(gcd, [tgdc, V1[i]]);
          end do;

          return simplify(V1/~tgdc);
     end proc,

     proc(M::{Matrix})
          option overload;
          uses LinearAlgebra;
          local i, num,r,c, tgdc,dnm, V1, Ml;
          r,c:=Dimension(M);
          num:=r*c;
          V1:=convert(M,list);
          dnm:=frontend(lcm, [seq(denom(V1[i]),i=1..num)]);
          Ml:=radnormal(dnm*~M);
          V1:=convert(Ml,list);#print((dnm,V1));
          tgdc:=V1[1];#print("xx")

          for i from 2 to num do
               tgdc:=frontend(gcd, [tgdc, V1[i]])
          end do;

          return simplify(Ml/~tgdc);
     end proc,

     proc(l::{`+`,`*`,`=`, `symbol`,procedure}, {vars::list:=[:-x,:-y]})
          option overload;
          uses LinearAlgebra;
          local i, num,f1,f1a,lv,lr, tgdc,dnm, V1,Vs;
          f1 := `if`(l::procedure, l(vars[]), l);
               f1a:=`if`(f1::`=`,lhs(f1)-rhs(f1),f1) ; # Remequal(f1);
          lr:=primpart(f1a,vars);
          return lr
end proc

]):

>	ProjVectorC := proc(a, b, c) local cfs, vectr; description " A Projective Column (Line) Vector in Reduced format"; cfs := FactReduce([a, b, c]); vectr := <[<cfs>]>; end proc:

>

>	ProjVectorR := proc(a, b, c) local cfs, vectr; description " A Projective Row (Point) Vector"; cfs := sign(c)*FactReduce([a, b, c]); vectr := <[<cfs>^%T]>^%T; end proc:

>

>	`&otimes;` := proc(A, B) local cp; description "Cross Product of Projective Vectors in Reduced format"; cp :=sign(c)* FactReduce(LinearAlgebra:-`&x`(A[1], B[1]))^%T; cp := ifelse(cp[3] <> 0, <[sign(cp[3]) *~ cp]>, cp); #makes sure format is [x,y,z] and not [x,y-z] end proc: