I meet a problem in the  Lipschitz–Hankel integral as follows,

Int(x^-2*BesselJ(1,x)*BesselJ(1,0.6*X)*e^(-10*x/7),x=0..infinity)

The output is always the same as shown above rather than a float value. How can I deal with such a integral?

Many Thanks!

I am displaying in 2D, multiple arrows with amplitudes oscillating sinusoidally.  I wish each arrow to have its own color. I have included the commands I am using. Please provide the command that will provide the colors. Thanks.

Ax:=t->cos(5*t);Ay:=t->0;Bx:=t->-cos(5*t)/2;By:=3*cos(5*t)/2;

animate(arrow,[[[0.0],[0.0]],[[Ax(t),Ay(t)],[Bx(t),By(t)]]],t=0..Pi,view=[-2..2,-2..2],color=[blue,green])

Is there a way to create an interactive plot, and put it in a Maplet, so I could get onclick coordinates and display some points there on my plot ? I checked Maple help and I didn't find any onclick event for plot or Maplets,Elements,Plotter but I found something called PlotComponent and there are some options like clickx, clicky. Can I use it somehow in Maplet, or maybe there's another solution ?

Hi all

Im trying to calculate scientific error propagation with data from Excel. Unfortunately Im having troubles making the package work. What I want to do is import a list from Excel with values, and another list of the error of each value.

The problem comes when using the Combine function on the imported arrays from excel.
Could anyone help me work this out? I feel it should be possible for Maple to import tables from Excel, calculate the propagated error. ...

Hi, below is what I have. I am trying to change elenum:=Array([2,4,6,8,16])

Then I changed my E0:=evalf((map(c->c*((map(int,((B^%T).(B1).J),eta=0..1))),L))

I also changed my E2:=evalf(map(c->c*map(int,((B2^%T).(B2).J),eta=0..1)),(1/(2*L))))

My E0M:=Matrix(nodes), which I am not sure if it is correct.

Then I tried to use the for do statements, it would not allow me to do so. What i'm trying to get is the first one to be a 5x5 matrix, followed...

Hi, I have a array with 5 numbers in it. I would like to multiply it into a 3x3 matrix, how do i do that?

for example A:=([3,5,7,9,15])

B:=Matrix([3x+4,2x-9,42],[8x+2,93x^2,6y],[5x/y,34y,32y^2])

how do i get 5 sets of output when i multiply it together? And can i get pick out each outpur as well? for example, I want the 3rd output which is the matrix multiply with 7, can I pick it out seperately if I wanted?

Please advice. Thanks

Dear all,

Some time ago I asked a question concerning the algorythm for numerical dsolve and related issues here.

I got a nice answer from pagan and implemented the algorythm.

However, since my system is big, the iterations take long time.

Here is the code itself:

sys1:=[TRIC[1],TRIC[2],TRIC[3],TRIC[7],TRIC[8],TRIC[9]]:
> sys2:=[TRIC[4],TRIC[5],TRIC[6],TRIC[10],TRIC[11],TRIC[12]]:
> GRAD:=subs(u=FF[1](t),a=FF[2](t),k=k(t),tau=tau(t),psi[1...

Hi,

I have a 1x3 matrix, with s0,s1,s2 terms in it. For example :
[(s-s1)*(s-s2). (s-s0)*(s-s2), (s-s0)*(s-s1)]

Is there a simple way that i can integrate the first one by s, second one by s0 and third one by s1?

to get a matrix for example:

[(int((s-s1)*(s-s2),s). int((s-s0)*(s-s2),s0), int((s-s0)*(s-s1),s1)]

Please advice,

Thanks.

Let S1 be the set of solutions of the inequalities x^2+ 3*x + 2 >0 and S2  be the set of solutions of the inequalities x^2 - 9>=0. I want to the set S1 union S2. What must i do? Please help me. 

I converted the package,module to a worksheet then executed the example before realizing it all executes at once. 

MyPackage := module()
        export  f1, f2;
        local   loc1;
        option package;
        f1 := proc() loc1 end proc;
        f2 := proc( v ...

I have four equations:

eqs:={8x-5z=0, -5y-2x=0, 2z-8y=0} and 13x-4y+9z=1

I've been trying for the past hour to solve both equations for x, y, and z and have not been successful.

What commands do I use to solve this problem?

I would like to save a simple plot to jpeg file on my hard disk. I've got a Maplet, and a Button which looks like that

Button['SIM']
 (
     "Save image",
     onclick=Evaluate('function'='saveImage()')
)

and a procedure:

saveImage := proc()

   plotsetup('jpeg', 'plotoutput'="plot.jpeg", 'plotoptions'="portrait,noborder");
   plot3d(sin(x)+cos(y), x=0..2*Pi, y=0..2*Pi, 'axes'='boxed');

I was solving one Math puzzle and found out a strange behaviour of the Optimization package commands in the Classic Worksheet.

Let us minimize the function 7*x+3*y under the constraint 6*x+2*y >= 49 for non-negative integers. When I type

Optimization:-Minimize(7*x+3*y, {6*x+2*y >= 49}, assume = nonnegint);

Let function f(x)=x(2011+√(2013-x²)).Find maximum and minimum value of function f(x)?

Why won't isolve show all the solutions here?

> a := 6*x+2*y >= 48:
> b := 6*x+2*y <= 60:
>
> with(plots):
> aa := implicitplot(a, x = 0 .. 10, y = 0 .. 5, filledregions = true, transparency = .5):
> bb := implicitplot(b, x = 0 .. 10, y = 0 .. 5, filledregions = true, transparency = .5):

display(aa,bb)

We can see the region of solutions.  But if we use isolve ...