I am evaluating a difficult integral of an oscillating function.

f:=k->evalhf(piecewise(k=0,1,k>0,cos(k*(dz))*BK0(k*rho1)/BK0(k*a)));

evalf(Int(f, kmin..kmax, method = _d01akc,epsilon=eps,maxintervals=1000000))

with dz>0 and possibly very large
kmin is Pi/(2*dz)
the number of cycles is nc = ceil((1/2)*kl*dz/Pi) with kl = 91.563
kmax=  kmin+nc*2*Pi/dz

Each evaluation of this integral using a call to the NAG routine (_d01akc ) eats up 

I'm using the printf command and i want my output a decimal number to get it published

s1=0.5555552487 and f(s1)=0.2542365647 what should i do?

hi;

when run file , maple give wrong result, plz help me

my file is attached

thx

 EKM3-cylindrical_pan.mw

I have solved  the following BVP system in Maple

restart; e := 12; Sc := .1; Pr := .7; Gr := 1; M := .1; sys := diff(f(t), `$`(t, 3))+2*f(t)*(diff(f(t), `$`(t, 2)))-(diff(f(t), t))^2-M^2*(diff(f(t), t))+Gr*theta(t) = 0, diff(theta(t), `$`(t, 2))+2*Pr*f(t)*(diff(theta(t), t)) = 0, diff(phi(t), `$`(t, 2))+Sc*(2*f(t)*(diff(phi(t), t))-phi(t)) = 0; bcs := f(0) = .1, (D(f))(0) = 0, theta(0) = .5, phi(0) = 1, (D(f))(e) = 0, theta(e) = 0, phi(e) = 0; sol := dsolve({bcs, sys}, numeric); plots[odeplot...

For educational purposes, I would like to hide some code (procedures) in a maple file. I put these procedures in a subsection which I can collapse. But everytime you execute the worksheet, the subsection is automatically expanded and I have to collapse the subsection again.
Is there a way to avoid the automatic expansion of subsection while executing the worksheet? 

If we have a matrix with n rows and 1 column of n-variable function  how to use the subs?

How to implement in maple, the differentiation matrix of second order ?
Is it possible to implement it with the command LinearAlgebra[BandMatrix]?

I can not get the matrix:



Lets look at an example using the following set of 8 discrete points (x,y)

(2,5),(4,15),(6,35),(8,65),(10,100),(12,145),(14,195), and (16,255).

Hi,

I need to extract the extremes of a cyclical function, actually I have all the data (x,y) and I was wondering whether extracting these extremes as I have shown in the attached figure will be easy with Maple or not? NOTICE that the attached figure is actually 2 series of data overlaid, you just consider the blue one and ignore the red one

TIA

Hi,

I wanted to know if it possible to get maple to skip all operations to the bottom if it receives an error anyway through my program? I have a number of seperate lines and don't want to put it into one big if statement.

Thanks

So, i'm using Maple 13 in my school, and our homework was to make two procedures which could calculate something.

I'm interested in making a procedure which can calculate out the angles of my triangle with the law of cosines, and this is how my procedure looks so far:

findA := proc (a, b, c)

local solA;

solA := arccos((b^2+c^2-a^2)/(2*b+2*c));

return evalf(solA);

end proc

findA(67, 50, 34)

After many calculations tried to run eval(exp). Was  in shock :( Big table used. Is there smth like forget(all) to forget all functions. After forget(exp) memory of course lowers enought.

Hello..

I need to write a function in maple which creates matrices from 1x1 to 20x20.. And the filling in the matrix needs to be 1 in the diagonal. But not like the diagonal in the identity matrix, but the opposite..: see picture:

So the amounts of 1 increase depending on the size of the matrix...

I know that i can put the function in a for loop to do the steps 20 times...

CTmodel-1.mw

Can anybody spot the problem in the second solve() in the attached file?  I am using the same approach that I did for the first solve(), but Maple is returning any empty list. 

Thanks!

I've set up part of my procedure to calculate the values of i from 1 to 15 using for loops.  I only need the values from 1 to 15, but a later part of my procedure if failing because it says it cannot evaluate the value for i=16.  How do I tell Maple that i>=16 doesn't matter/do nothing for i>= 16?