MaplePrimes Questions

I'm trying to be able to add bullet points to this Plot as well as where they cross. Download provided below. Thanks in Advance.

Example 7: Equations of Perpendicular Lines


NULLy = -(3/2)*x+2


2*x-3*y = 5; "_noterminate"

2, -1; "_noterminate"and m = 2/3; "_noterminate"

y = -5/3+(2/3)*x





Why the performance of Maple GUI is bad. It is really hard to type anything in thew gui. 

I have a relatively involved modeling application that tracks  a particle beam through a beam line for multiple passes. The beam is described by a Maple Record with a number of different entries; one of them being a table with a 6-vector of numbers for each particle in the beam.

Each pass this beam gets updated. I would like to be able to save the whole record into a filel (i.e. append the file each pass) and then be able to read it back in; regenerating the Record structures. Is this possible?

(I'd normally upload a mwe; except this code uses a large package that would have to be installed. This is something I do not consider reasonable to ask for.)

If there is an explicit answer like: call this function, give it your Record and it should work; and call this reading function to recreate the Record, that would be swell. Or maybe some pointers in the right direction...Maple has so many read, write, save and restore functions that a systematic test is daunting.


Mac Dude



Digits := 10:




v := .7:

Disp := 20:

esp := 800000:

k := 0:

E := proc (x, t) options operator, arrow; int(exp((-esp*w^4+Disp*w^2+k)*t)*cos(w*(x+v*t))/Pi, w = 0. .. infinity) end proc;

proc (x, t) options operator, arrow; int(exp((-esp*w^4+Disp*w^2+k)*t)*cos(w*(x+v*t))/Pi, w = 0. .. infinity) end proc


f := proc (x) options operator, arrow; 20*exp(-(1/2000000)*(x-10000)^2)+15*exp(-(1/2000000)*(x-13800)^2) end proc:



u := proc (x, t) options operator, arrow; int(E(x-xi, t)*f(xi), xi = 0. .. 20000) end proc;

proc (x, t) options operator, arrow; int(E(x-xi, t)*f(xi), xi = 0. .. 20000) end proc



plot(u(x, t), x = 1500, t = 0 .. 60000, numpoints = 100)

Error, (in plot) unexpected options: [x = 1500, t = 0 .. 60000]









I'm reading this book : Algebra and Trigonometry - Real Mathematics, Real People, 7th Edition

Chapter 1 Functions and Their Graphs, Example 7, Monthly Wage , Page 80.

I'm trying to plot for a monthly wage.

The inital equation is:

y = 2000 + 0.1*x

Then other changes:

a. Sales are $1480 in August. What are your wages for that month?
b. You receive $2225 for September. What are your sales for that month?

I'm lost. Which is better Numerical Solution or Graphical Solution?

Thanks in Advance,




I want to have the fractal of von koch as shown in the photo, and then animate it? the command Kochcurve does not give exactly the objective. some ideas?


f := proc (a, b) options operator, arrow; (1/2)*b+(1/2)*arccos(sin(2*a-b)/tan(b)) end proc; dis := proc (A, B) options operator, arrow; sqrt(inner(A-B, A-B)) end proc; bisA := proc (A, B, C) local b, c, M; b, c := dis(A, C), dis(A, B); M := (b*B+C*c)/(b+c); x*(A[2]-M[2])+y*(M[1]-A[1])+A[1]*M[2]-A[2]*M[1] end proc; P := proc (a0, b0) local a, b, c, p1, p2, p3, p4, r, II; a := evalf(a0); b := evalf(b0); c := f(a, b); if b0-0.1e-2


The VectorCalculus package offers a convenient way to compute multivariate integral over a non-square geometry. See help("VectorCalculus/int"). Here is an example taken from this help page:



int( x*y, [x,y] = Triangle( <0,0>, <1,0>, <0,1> ) );


The question is: is there a way to draw the area “Triangle( <0,0>, <1,0>, <0,1> )” in order to quickly check that it matches the expected region of integration that one needs?

I am trying to plot the solution of the DE in the interval [0.3,4]

  I have read the article attached in the following

I wrote the code mentioned in the article in order to find the Hirota formula This is the code that I wrote For Hirota Method But when running the program it does not work well 

Could you help me to fix the mistake and run the program


restart; with(PDEtools); with(DEtools)

alias(u = u(x, t)); declare(u(x, t)); alias(f = f(x, t)); declare(f(x, t))

u, f


` u`(x, t)*`will now be displayed as`*u


u, f


` f`(x, t)*`will now be displayed as`*f



BD := proc (FF, DD) local f, g, x, m, opt; if nargs = 1 then return `*`(FF[]) end if; f, g := FF[]; x, m := DD[]; opt := args[3 .. -1]; if m = 0 then return procname(FF, opt) end if; procname([diff(f, x), g], [x, m-1], opt)-procname([f, diff(g, x)], [x, m-1], opt) end proc

BD([f(x, t), g(x, t)], [x, 1])

(diff(f(x, t), x))*g(x, t)-f(x, t)*(diff(g(x, t), x))


BD([f(x, t), g(x, t)], [t, 1])

(diff(f(x, t), t))*g(x, t)-f(x, t)*(diff(g(x, t), t))


BD([f(x, t), g(x, t)])

f(x, t)*g(x, t)




print_BD := proc (FF, DD) local f, g, x, m, i; f, g := FF[]; f := cat(f, '', g); g := product(D[args[i][1]]^args[i][2], i = 2 .. nargs); if g <> 1 then f := g*(`*`(f)) end if; f end proc