Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

A few people have asked me how I created the sections in the Maple application in this video: https://youtu.be/voohdmfTRn0?t=572

Here's the worksheet (Maple 2016 only). As you can see, the “sections” look different what you would normally expect (I often like to experiment with small changes in presentation!)

These aren't, however, sections in the traditional Maple sense; they're a demonstration of Maple 2016's new tools for programmatically changing the properties of a table (including the visibility of its rows and columns). @dskoog gets the credit for showing me the technique.

Each "section" consists of a table with two rows.

  • The table has a name, specified in its properties.
  • The first row (colored blue) contains (1) a toggle button and (2) the title of each section (with the text in white)
  • The second row (colored white) is visible or invisible based upon the state of the toggle button, and contains the content of my section.

Each toggle button has

  • a name, specified in its properties
  • + and - images associated with its on and off states (with the image background color matching the color of the first table row)
  • Click action code that enables or disables the visibility of the second row

The Click action code for the toggle button in the "Pure Fluid Properties" section is, for example,

tableName:="PureFluidProperties_tb":
buttonName:="PureFluidProperties_tbt":
if DocumentTools:-GetProperty(buttonName, 'value') = "false" then   
     DocumentTools:-SetProperty([tableName, 'visible[2..]', true]);
else
     DocumentTools:-SetProperty([tableName, 'visible[2..]', false]);
end if;

As I said at the start, I often try to make worksheets look different to the out-of-the-box defaults. Programmatic table properties have simply given me one more option to play about with.

using worksheet mode, is there a way to automatically color any comments after (#) with a different color than the default red?

I have been using Maple for years and did not need to significantly document or comment on my worksheets before because I only needed to share my worksheets with close collegues who have a lot of experience with Maple too. However, now I need to share my code with a general audience that might not be a Maple user. So I need to add enough comments after each line, and I have been manually changing the color of comments after the # to green to give the reader the indication that this is a comment and not part of the code, like this for example:

 

> 

restart;

> 

f:=(x,y)->sin(sqrt(x^2+y^2))/sqrt(x^2+y^2); # Define f as a function of the variables x and y

proc (x, y) options operator, arrow; sin(sqrt(x^2+y^2))/sqrt(x^2+y^2) end proc

(1)
> 

df:=(x,y)->eval(diff(f(a,y),a),a=x); # Define df as the partial derivative of f with respect to x

proc (x, y) options operator, arrow; eval(diff(f(a, y), a), a = x) end proc

(2)
> 

df(1,3); # Evaluate df at x=1, y=3

(1/10)*cos(10^(1/2))-(1/100)*sin(10^(1/2))*10^(1/2)

(3)
> 

 

 

Download Worksheet-coloring-comments.mw



It would be great if Maple can automatically color comments like typical editors do. Is this possible? if not, does anyone have an advice for a nice easy way to add proper documentation to worksheets?

Hello everybody.

I'm trying to obtain the numerical solution of a differential equation. Unfortunately, this prove to be quite challenging. I was able to obtain a rough solution using mathematica, but nothing more. The function is strictly increasing (for sure).

Any help is really REALLY appreciated, thanks!

 

``

deq1 := 1/(b-f(b)) = (2*(3-(1-f(b)*(diff(f(b), b, b)))/((diff(f(b), b))*(diff(f(b), b)))))/(1-2*(b-(1-f(b))/(diff(f(b), b))))

1/(b-f(b)) = 2*(3-(1-f(b)*(diff(diff(f(b), b), b)))/(diff(f(b), b))^2)/(1-2*b+2*(1-f(b))/(diff(f(b), b)))

(1)

ic1 := eval(f(b), b = 3/8) = 0, eval(f(b), b = 1/2) = 1/2

f(3/8) = 0, f(1/2) = 1/2

(2)

digits := 3

3

(3)

dsol1 := dsolve({deq1, ic1}, method = bvp[middefer], numeric, range = 3/8 .. 1/2)

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

 

``

 

Download diffeqn.mw

Hi everybody

In the following attached file that is a simple code, 2 problems occur:

1- The values of Lc1 and Lc2 are specified in the third and fourth lines of the code. when I execute Pcl10, Maple does not replace Lc1 with its value that is 0.5*L. This problem does not happen for Tcl10 and Tcl21. Why Maple does not replace Lc1 with its value that is specified at the start of the worksheet?

2- In the last line, when I try to simplify the Tcl10, Maple returns an error, while for example the first element i.e. Tcl20(1,1) can be simplified as: cos(thet1+thet2). What is the source of error?

 

Thanks in advance

P1.mw

I had another problem whereby I need to shaded a region between y=ln(x+2) , y-axis and y=-1 to y= 2. What is the suitable comand I shall apply ( Maple 18)


Trying to build a block matrix. Having a problem getting the syntax correct. Can't add in the predefind matrices to the lower band.

restart

with(LinearAlgebra):

interface(displayprecision = 5)

5

(1)

interface(rtablesize = 81)

10

(2)

S := 2

2

(3)

dmax := 7

7

(4)

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

(5)

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

(6)

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

(7)

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

(8)

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

(9)

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

(10)

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

(11)

S*dmax

14

(12)

M1 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), Matrix(S, S*dmax, [seq(-CCnew || n, n = 0 .. dmax-1)])])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -336750255587/3769550688757, (1, 2) = 14853552191797/1696297809940650, (1, 3) = -49655436033349/56543260331355, (1, 4) = 75647656451147/1413581508283875, (1, 5) = -299962512141959/80776086187650, (1, 6) = 1231816081155781/8481489049703250, (1, 7) = -50445725001719/5769720441975, (1, 8) = 9065291388901/40388043093825, (1, 9) = -142685068141037/11539440883950, (1, 10) = 116560067351321/565432603313550, (1, 11) = -60560690824604/5769720441975, (1, 12) = 88774498025543/848148904970325, (1, 13) = -20398649489879/4121228887125, (1, 14) = 5446073753219/242328258562950, ...

 

M2 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity)])

M2 := Matrix(14, 14, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 1, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (1, 10) = 0, (1, 11) = 0, (1, 12) = 0, (1, 13) = 0, (1, 14) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 1, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (2, 10) = 0, (2, 11) = 0, (2, 12) = 0, (2, 13) = 0, (2, 14) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (3, 5) = 1, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (3, 10) = 0, (3, 11) = 0, (3, 12) = 0, (3, 13) = 0, (3, 14) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 0, (4, 5) = 0, (4, 6) = 1, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (4, 10) = 0, (4, 11) = 0, (4, 12) = 0, (4, 13) = 0, (4, 14) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (5, 7) = 1, (5, 8) = 0, (5, 9) = 0, (5, 10) = 0, (5, 11) = 0, (5, 12) = 0, (5, 13) = 0, (5, 14) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0, (6, 7) = 0, (6, 8) = 1, (6, 9) = 0, (6, 10) = 0, (6, 11) = 0, (6, 12) = 0, (6, 13) = 0, (6, 14) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 0, (7, 5) = 0, (7, 6) = 0, (7, 7) = 0, (7, 8) = 0, (7, 9) = 1, (7, 10) = 0, (7, 11) = 0, (7, 12) = 0, (7, 13) = 0, (7, 14) = 0, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 0, (8, 5) = 0, (8, 6) = 0, (8, 7) = 0, (8, 8) = 0, (8, 9) = 0, (8, 10) = 1, (8, 11) = 0, (8, 12) = 0, (8, 13) = 0, (8, 14) = 0, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = 0, (9, 6) = 0, (9, 7) = 0, (9, 8) = 0, (9, 9) = 0, (9, 10) = 0, (9, 11) = 1, (9, 12) = 0, (9, 13) = 0, (9, 14) = 0, (10, 1) = 0, (10, 2) = 0, (10, 3) = 0, (10, 4) = 0, (10, 5) = 0, (10, 6) = 0, (10, 7) = 0, (10, 8) = 0, (10, 9) = 0, (10, 10) = 0, (10, 11) = 0, (10, 12) = 1, (10, 13) = 0, (10, 14) = 0, (11, 1) = 0, (11, 2) = 0, (11, 3) = 0, (11, 4) = 0, (11, 5) = 0, (11, 6) = 0, (11, 7) = 0, (11, 8) = 0, (11, 9) = 0, (11, 10) = 0, (11, 11) = 0, (11, 12) = 0, (11, 13) = 1, (11, 14) = 0, (12, 1) = 0, (12, 2) = 0, (12, 3) = 0, (12, 4) = 0, (12, 5) = 0, (12, 6) = 0, (12, 7) = 0, (12, 8) = 0, (12, 9) = 0, (12, 10) = 0, (12, 11) = 0, (12, 12) = 0, (12, 13) = 0, (12, 14) = 1, (13, 1) = 0, (13, 2) = 0, (13, 3) = 0, (13, 4) = 0, (13, 5) = 0, (13, 6) = 0, (13, 7) = 0, (13, 8) = 0, (13, 9) = 0, (13, 10) = 0, (13, 11) = 0, (13, 12) = 0, (13, 13) = 0, (13, 14) = 0, (14, 1) = 0, (14, 2) = 0, (14, 3) = 0, (14, 4) = 0, (14, 5) = 0, (14, 6) = 0, (14, 7) = 0, (14, 8) = 0, (14, 9) = 0, (14, 10) = 0, (14, 11) = 0, (14, 12) = 0, (14, 13) = 0, (14, 14) = 0})

(13)

lowerband := Matrix(S, S*dmax, [seq(-evalf(CCnew || n), n = 0 .. dmax-1)])

lowerband := Matrix(2, 14, {(1, 1) = -0.893343e-1, (1, 2) = 0.875645e-2, (1, 3) = -.878185, (1, 4) = 0.535149e-1, (1, 5) = -3.71351, (1, 6) = .145236, (1, 7) = -8.74318, (1, 8) = .224455, (1, 9) = -12.36499, (1, 10) = .206143, (1, 11) = -10.49630, (1, 12) = .104669, (1, 13) = -4.94965, (1, 14) = 0.224740e-1, (2, 1) = -0.864552e-2, (2, 2) = -.115162, (2, 3) = -0.776011e-1, (2, 4) = -1.10141, (2, 5) = -.282396, (2, 6) = -4.51408, (2, 7) = -.537873, (2, 8) = -10.26209, (2, 9) = -.568850, (2, 10) = -13.97156, (2, 11) = -.318088, (2, 12) = -11.39493, (2, 13) = -0.737097e-1, (2, 14) = -5.15750})

(14)

M3 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), lowerband])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -0.8933432215e-1, (1, 2) = 0.8756453086e-2, (1, 3) = -.8781848755, (1, 4) = 0.5351488825e-1, (1, 5) = -3.713506389, (1, 6) = .1452358276, (1, 7) = -8.743183575, (1, 8) = .224454831, (1, 9) = -12.36498974, (1, 10) = .2061431666, (1, 11) = -10.49629552, (1, 12) = .1046685287, (1, 13) = -4.949652167, (1, 14) = 0.2247395242e-1, (2, 1) = -0.8645515381e-2, (2, 2) = -.1151624586, (2, 3) = -0.7760109089e-1, (2, 4) = -1.101411529, (2, 5) = -.2823958356, (2, 6) = -4.514076892, (2, 7) = -.537872817, (2, 8) = -10.26209174, (2, 9) = -.5688498906, (2, 10) = -13.97155838, (2, 1...

 

``


Download Test_Block_matrix.mw

restart

with(LinearAlgebra):

interface(displayprecision = 5)

5

(1)

interface(rtablesize = 81)

10

(2)

S := 2

2

(3)

dmax := 7

7

(4)

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

(5)

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

(6)

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

(7)

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

(8)

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

(9)

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

(10)

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

(11)

S*dmax

14

(12)

M1 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), Matrix(S, S*dmax, [seq(-CCnew || n, n = 0 .. dmax-1)])])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -336750255587/3769550688757, (1, 2) = 14853552191797/1696297809940650, (1, 3) = -49655436033349/56543260331355, (1, 4) = 75647656451147/1413581508283875, (1, 5) = -299962512141959/80776086187650, (1, 6) = 1231816081155781/8481489049703250, (1, 7) = -50445725001719/5769720441975, (1, 8) = 9065291388901/40388043093825, (1, 9) = -142685068141037/11539440883950, (1, 10) = 116560067351321/565432603313550, (1, 11) = -60560690824604/5769720441975, (1, 12) = 88774498025543/848148904970325, (1, 13) = -20398649489879/4121228887125, (1, 14) = 5446073753219/242328258562950, ...

 

M2 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity)])

M2 := Matrix(14, 14, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 1, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (1, 10) = 0, (1, 11) = 0, (1, 12) = 0, (1, 13) = 0, (1, 14) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 1, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (2, 10) = 0, (2, 11) = 0, (2, 12) = 0, (2, 13) = 0, (2, 14) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (3, 5) = 1, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (3, 10) = 0, (3, 11) = 0, (3, 12) = 0, (3, 13) = 0, (3, 14) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 0, (4, 5) = 0, (4, 6) = 1, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (4, 10) = 0, (4, 11) = 0, (4, 12) = 0, (4, 13) = 0, (4, 14) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (5, 7) = 1, (5, 8) = 0, (5, 9) = 0, (5, 10) = 0, (5, 11) = 0, (5, 12) = 0, (5, 13) = 0, (5, 14) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0, (6, 7) = 0, (6, 8) = 1, (6, 9) = 0, (6, 10) = 0, (6, 11) = 0, (6, 12) = 0, (6, 13) = 0, (6, 14) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 0, (7, 5) = 0, (7, 6) = 0, (7, 7) = 0, (7, 8) = 0, (7, 9) = 1, (7, 10) = 0, (7, 11) = 0, (7, 12) = 0, (7, 13) = 0, (7, 14) = 0, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 0, (8, 5) = 0, (8, 6) = 0, (8, 7) = 0, (8, 8) = 0, (8, 9) = 0, (8, 10) = 1, (8, 11) = 0, (8, 12) = 0, (8, 13) = 0, (8, 14) = 0, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = 0, (9, 6) = 0, (9, 7) = 0, (9, 8) = 0, (9, 9) = 0, (9, 10) = 0, (9, 11) = 1, (9, 12) = 0, (9, 13) = 0, (9, 14) = 0, (10, 1) = 0, (10, 2) = 0, (10, 3) = 0, (10, 4) = 0, (10, 5) = 0, (10, 6) = 0, (10, 7) = 0, (10, 8) = 0, (10, 9) = 0, (10, 10) = 0, (10, 11) = 0, (10, 12) = 1, (10, 13) = 0, (10, 14) = 0, (11, 1) = 0, (11, 2) = 0, (11, 3) = 0, (11, 4) = 0, (11, 5) = 0, (11, 6) = 0, (11, 7) = 0, (11, 8) = 0, (11, 9) = 0, (11, 10) = 0, (11, 11) = 0, (11, 12) = 0, (11, 13) = 1, (11, 14) = 0, (12, 1) = 0, (12, 2) = 0, (12, 3) = 0, (12, 4) = 0, (12, 5) = 0, (12, 6) = 0, (12, 7) = 0, (12, 8) = 0, (12, 9) = 0, (12, 10) = 0, (12, 11) = 0, (12, 12) = 0, (12, 13) = 0, (12, 14) = 1, (13, 1) = 0, (13, 2) = 0, (13, 3) = 0, (13, 4) = 0, (13, 5) = 0, (13, 6) = 0, (13, 7) = 0, (13, 8) = 0, (13, 9) = 0, (13, 10) = 0, (13, 11) = 0, (13, 12) = 0, (13, 13) = 0, (13, 14) = 0, (14, 1) = 0, (14, 2) = 0, (14, 3) = 0, (14, 4) = 0, (14, 5) = 0, (14, 6) = 0, (14, 7) = 0, (14, 8) = 0, (14, 9) = 0, (14, 10) = 0, (14, 11) = 0, (14, 12) = 0, (14, 13) = 0, (14, 14) = 0})

(13)

lowerband := Matrix(S, S*dmax, [seq(-evalf(CCnew || n), n = 0 .. dmax-1)])

lowerband := Matrix(2, 14, {(1, 1) = -0.893343e-1, (1, 2) = 0.875645e-2, (1, 3) = -.878185, (1, 4) = 0.535149e-1, (1, 5) = -3.71351, (1, 6) = .145236, (1, 7) = -8.74318, (1, 8) = .224455, (1, 9) = -12.36499, (1, 10) = .206143, (1, 11) = -10.49630, (1, 12) = .104669, (1, 13) = -4.94965, (1, 14) = 0.224740e-1, (2, 1) = -0.864552e-2, (2, 2) = -.115162, (2, 3) = -0.776011e-1, (2, 4) = -1.10141, (2, 5) = -.282396, (2, 6) = -4.51408, (2, 7) = -.537873, (2, 8) = -10.26209, (2, 9) = -.568850, (2, 10) = -13.97156, (2, 11) = -.318088, (2, 12) = -11.39493, (2, 13) = -0.737097e-1, (2, 14) = -5.15750})

(14)

M3 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), lowerband])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -0.8933432215e-1, (1, 2) = 0.8756453086e-2, (1, 3) = -.8781848755, (1, 4) = 0.5351488825e-1, (1, 5) = -3.713506389, (1, 6) = .1452358276, (1, 7) = -8.743183575, (1, 8) = .224454831, (1, 9) = -12.36498974, (1, 10) = .2061431666, (1, 11) = -10.49629552, (1, 12) = .1046685287, (1, 13) = -4.949652167, (1, 14) = 0.2247395242e-1, (2, 1) = -0.8645515381e-2, (2, 2) = -.1151624586, (2, 3) = -0.7760109089e-1, (2, 4) = -1.101411529, (2, 5) = -.2823958356, (2, 6) = -4.514076892, (2, 7) = -.537872817, (2, 8) = -10.26209174, (2, 9) = -.5688498906, (2, 10) = -13.97155838, (2, 1...

 

``


Download Test_Block_matrix.mw


Dear all,

I wold like to find the solution of the next system of two equations with three unknowns but we assume that the unknows are positive integers. How the following code can work. Many thanks

 

 

 

> restart;
> assume(J, integer, J >= 0);
> assume(A, integer, A >= 0);
> assume(T, integer, T >= 0);
> eq1 := J+10*A+50*T=500;
   eq2 := J+A+T = 100;
  solve( {eq1,eq2},{J,A,T});

 

For many years I succesfully used the command:

Lambda := ImportMatrix(LambdaFile, delimiter = "", datatype = string);

 

but for some txt files I now got the following error in Maple 2016:
Error, (in ImportMatrix) cannot interpret file

The same code still runs without problems at all in previous versions of Maple.

The strange thing is that other txt files can be still be imported in Maple 2016 without problems.

I will upload two files:

ULSlx3.txt returns an error message

ULStd3.txt no problems

ULSlx3.txt

ULStd3.txt

complete code:

restart;
with(LinearAlgebra); with(Statistics);
interface(rtablesize = infinity);
with(Typesetting):
interface(typesetting=extended):

LambdaFile := "C:/path/ULSlx3.txt";
Lambda := Matrix(5, 1);

in the whatsnew there is no information on changes with respect to the ImportMatrix command, although there is some (unrelated) information on Importing Data:

Importing Data
By default, the Import command now returns a DataFrame when importing from Excel, CSV, DIF, and TSV file formats, and returns a DataSeries when importing from ODS, SXC, and TSV file formats. Use the output option to specify a different format. For example, Import(origin, output=Matrix).

 

I have no clue what is wrong with my code and why some text files can be imported as before whereas other files returns an error message.

greetings Harry

In using Maple there are many ways to complete a command on a function. Using shortcuts in the worksheet/document you can enter the command on the line rather than using the drop down menus or clickable math menu.  I would like to see what the shortcut command are for the packages in Maple.  There are so many, I am using the LinearAlgebra package and have seen Dr. Lopez use " ||A||2 to determine the Euclidean Norm of A.  But where does this use get described in the help or examples of other shortcut usage for Maple. I've seen the shortcut commands for Documents in General, but the use of shortcut commands for other packages  don't seem to discuss the usage of such shortcuts.

Is there a document that documents these shortcuts.

 

Hi my dear friends, sorry for boring you. I am hauted by a problem about 'mtalyor'.

I failed to expand the following equation by using the following command:

u:=x(t)/sqrt(x(t)^2+y(t)^2-2*x(t)+1)*diff(y(t),t);
indets(u,name);
w:=evalindets(u,function(identical~({x(t),y(t)})),s->op(0,s)(freeze(op(s))));
var:=indets(w,name) minus {t};
mtaylor(w,var,4);
thaw(%)

 Freeze is applied to the part in which the expansion will take place, obviously it didn't work well. 

Thank you in advance for taking a look. 

To gererate a random initail cofigration -1 or 1

spin = (-1).^(round(rand(N)));

 

for i=1:1000,

 

neighbours = circshift(spin, [ 0 1]) + ...

circshift(spin, [ 0 -1]) + ...

circshift(spin, [ 1 0]) + ...

circshift(spin, [-1 0]);

how to do this in maple?

Hello

I was looking for information on how to solve a differential equation in Maple, but I do not know how to start.

I have yet equations in symbolic form and I want to apply for the calculation of such a procedure ODE4 (Runge Kutta 4) with fixed step integration Tc and time calculations ongoing T.

Where to start?
1. This should be done in another sheet, or I could write command in next line under the symbolic transformations?
2. To enter numeric parameters that certain symbols were they equal?
3. How to perform iterative calculations?

Thank you in advance for your help

I am curious, can simplify/siderels be executed in mod p by some equivalent Maple function call?

Hi, I have the equations (below) with different parameters. I'd like to find out if is it possitble to:

(1) Create a loop that solves for values of on of the parameters, T, in the interval -0.3<T<0.3 and obtain the solutions allowing for a grid of 0.01.

That is, if one can first solve for T= -0.3 and then use the loop with increments that raise T by 0.01 each time until T=0.3. Other parameters are held constant at their assigned values.

(2) Generates the corresponding vectors for T and all the endogenous variables.

eq1 := PI_T = alpha*M*(kappa/(1+kappa-sigma)-1):
eq2 := l = alpha*((sigma-1)*kappa/(1+kappa-sigma)+1):
eq3 := M = phi_c^(-kappa)*F:
eq4 := e*F = PI_T*v+T:
eq5 := M*l+L_A = L_s:
eq6 := F*e+A = A_s:
eq7 := A = (1-beta)*(L_s+(1-v)*PI_T-T):
eq8 := A_s = L_A:
eq9 := p = sigma/((sigma-1)*phi):
eq10 := P_Y = M^(lambda/(1-sigma))*p:
eq11 := L_s = N*(theta*P_Y^beta)^(-1/delta):
eq12 := phi = (kappa/(1+kappa-sigma))^(1/(sigma-1))*phi_c:
eq13 := U = delta*theta*L_s/((1+delta)*P_Y^beta)+((1-v)*PI_T-T)/P_Y^beta:

Params:= [ T=0, N = 1, v = 1, beta = .75, lambda = 1, sigma = 5, kappa = 4.8, delta = 2, theta = 2.591350635, alpha = 0.3998699153e-1, e = 0.8333333332e-1]:

Init_Values:= {A_s = .2500000000, l = .9996747882, PI_T = 0.8333333332e-1, L_A = .2500000000, phi = 1.878101496, M = .4168022157, A = .1666666667, p = .6655657336, P_Y = .8283396488, L_s = 2/3, phi_c = 1.2, F = 1, U = 1.326439132 }:


SOL:= fsolve(eval({eq||(1..13)}, Params), Init_Values);

 

Many thanks!!

Hello there,

 

Suppose I have a parametric experssion like   P(ε)=1+ε+ε2+ε3 where ε is very small. How can I get P(ε)≈1+ε or P(ε)≈1+ε+ε2? 

 

Thanks

First 1136 1137 1138 1139 1140 1141 1142 Last Page 1138 of 2261