In Maple 11 I need to write define_external procedure with 23 arguments.

For example:
   test1 := define_external('myproc',
                                               N1::(integer[4]),

A few weeks ago I mentioned the ncrunch comparison of "mathematical programs for data analysis" in a comment in another thread.  There is now a new, 5th release of that review. The systems reviewed are:

• GAUSS
• Maple
• Mathematica
• Matlab
• O-Matrix
• Ox
• SciLab

The review is skewed towards statistical computation and data manipulation, but it includes several interesting comparisons of the major computer algebra systems (CAS).

There is a comparative performance section, and the worksheets used for that benchmarking are available for download. Here is the Maple worksheet, which was used with Maple 11.

Hi

 

how can i do this in maple

 

if(a>b>c) then 4

 

or maybe if(a>b) and if(b>c) then 4

All:

This should be simple but I am clearly missing some of the nuances - can you please help me see what I am doing wrong?  I would like to solve 2 eqns in 2 unknowns:

eq15 := Kzeff15 = (Kvr*Kz15)/(Kvr+Kz15);

eq33 := Kzeff33 = (Kvr*Kz33)/(Kvr+Kz33);

Kz33 := (3.3/1.5)*Kz15;

Kzeff15 := 150612;

Kzeff33 := 253629;

soln_z15 := fsolve ({eq15,eq33},{Kvr,Kz15});

 

 

Greetings all

I would like to show the x y values next to the point on the plot does maple have a command to do this?

data := [[0, 0], [.5, -1], [1, 0], [2, 2], [3, 0], [5, -4], [7, 0], [11, 8], [15, 0], [18.5, -7], [22, 0], [24.5, 5], [27, 0]];
 with(CurveFitting);
cubfit := BSplineCurve(data, x, order = 3);
plots[display](plot(data, style = point, color = black, symbol = BOX),
plot(cubfit, x = 0 .. 27, thickness = 2))

 

tia sal

What is the best way to print a maple file?  It seems that a lot of white space is added when making an html file.  Also it scrambles in a ton of strange symbols when its in a pdf, Im just trying to figure out how to print a 3 or 4 page maple worksheet without it printing on 20 pages!

 

Thanks!

How can I plot sin(y)>0, x>0 ?

I've tried using inequal:
with(plots):
inequal({sin(y)>0,x>0}, x=-2*pi..2pi, y=-2*pi..2*pi);

but I get the message 'unable to compute coeff' .

Thanks in advance.

I'm interested in finding a command that could stretch slope fields for me. Or perhaps I need to make a procedure of my own?

If that is the case, I would like to know how to draw a segment of a graph by varying the slope.

I think it is possible for me to use

with(geometry);

segment([y=y'*x+C, [x,y]],a,b);

But I don't have any idea how to make a program out of it yet.

 

Thanks in advance :)

Can I have maple show me the steps it takes to integrate?  Is this possible?

 

thanks!

Using the procedure permsPosNeg as given in the blog entry Positive and negative permutations, including the improvement obtained in the blog entry Refactoring Maple code, consider the following procedure:

Hi, I can't find any information what does FI mean, for example: f := proc(x) if x>0 then x else -x fi; end proc: f(-3); or: if iseuler(G) then continue else RETURN (print ("The graph does not have an Euler circuit")) fi; I'm tired of googling already, Maple Help doesn't have any information about fi, what might it be? perhaps some closing declaration instead of 'end if'? I would be grateful for any comments, thanks in advance, Sven

is there any way of creating a polynomial in maple from a list of coefficients (i could write my own function, but this is annoying, no?). specifically i would want an inverse to the PolynomialTools:-CoefficientList function, f([1, 2, 3]) would yield me 1 + 2*x + 3*x^2. i cant understand the apparent lack of such an elementary function, i must be missing someting? thanks

Hiya,

I have two questions...
Firstly when plotting a graph in Maple using plots[display]({A1,A2,A3}), is there anyway to switch the horizontal and vertical axes such that f(z) is on the horizontal axis and z along the vertical?

Secondly, I'm trying to plot a sum at various times t,

c(z,t):= Sum(A[n]*Phi[n](z)*exp(-lambda[n]*t), n = 1 .. 10)

I several peaks from emission spectra which I am trying to fit wit a Gaussian distribution. 

y = y_0 + A*exp[-{(x-x_0)/(width)}^2

 

The problem I am running into is I am getting fits which are wrong.  Maple returns:

fit_658 := 1.+1.*exp(.50000000000000*(z-1.)^2)

The correct answer is:

-7e05+2.7e06*exp(((z-656.1)/9)^2)

Also, I am interested in finding the deviation in the parameters which are found.

John Fredsted posted some interesting code dealing with permutations, and I suggested a small improvement.  Here, I want to continue the story of that improvement.  First, let us focus one particular line of code:

  posMaps,negMaps := seq(map((perm::list) ->
		(x::list) -> [seq(op(perm[i],x),i=1..nops(perm))]
	,perms),perms in [posPerms,negPerms]);

which uses a lexically scoped procedure to perform the permutations. The first thing to notice is that op(perm[i],x) is really equivalent to x[perm[i]]. Now that we have that perky op gone, we see that the resulting code expression will return unevaluated if x is unknown, unlike op which throws an error message (correctly so!). So now instead of using scoping, we can let Maple actually evaluate the inner perm[i] calls, and use unapply to recover a routine. Putting that together gives us my suggestion:

posMaps,negMaps := seq(map((perm::list) -> unapply([seq(x[perm[i]],i=1..nops(perm))],x),perms), 
perms in [posPerms,negPerms]);