A colleague is trying to evaluate some contour integrals numerically using evalf(Int(f,t=0..1)). Unfortunately f involves square and higher roots and Maple is getting an incorrect answer because it is not choosing the branch for these roots continuously through the range of integration.
Apparantly Mathematica has an option in numerical integration which forces a continuous choice of branch for the integrand. Is there any reasonably simple way to achieve the same with Maple?
Thanks.

Hello,
I am trying to do a 3d plot from 3 datapoints from a text file. I have the file DATA with 4000 data points in the following format:
-1 -1 2
-1 -0.976522 2
-1 -0.952609 2
-1 -0.915507 2
...
Where columns respectivly are the x, y, and z values that I want plotted. I figured that this command should plot the points:
listplot3d(readdata("DATA",3), axes=boxed, labels=[x, y, z]);
I do get a graph here, but instead of plotting the points x, y, z, I get the x value correspoindng to the line number, a y value that seems to be the value of the y column + 2, and the actual z value that I want. So, my question is, how do I just plot the three values that I am inputting?

Hi!
I have a lot of programs made by fortran77, can i use this programs under maple10?
thanks.

Hi,
How to solve 12 nonlinear equations with 6 unknowns with Maple 10 or other software but analytic, not numeric?
I was tried with solve(). After waiting 24h computation I was interrupt it. And I don't think that Maple will solve that problem.
In begining CPU usage was 100% but latter was 0% and HD always read or write something without stoping.

I want a function depending on n1 and n2, and a1 and a2, more or less in the following way:
sum( sum( ((if ((j1*k1/n1 +j2*k2/n2+k2/n2 > 1/2) &and
(j1*k1/n1 +j2*k2/n2 >1/2)) 1 else 0), j1=0..n1), j2=0..n2).
How do I write in maple the previous expresion?.
That is.. how may I write
MYfunction( n1, n2) = sum( sum( BOOLEAN , j1=0..n1),j2=0..n2);
where BOOLEAN stands for terms involving j1,j2,n1,n2 and other "constants"?
Thanks in advance!
Jacques

several beams are connected by gemel one after another, the 1st is fixed at one end on the roof, and the last is driven by a force(fx,fy vector summation) at the end point.
i am using lagrange method for the virtual displacement of angle theta[i] , and momentum balance in 2 direction (x,and y)to build the model.
in the last step, i got an ODE system of 2nd order, but i cant code in maple to solve it. i think it is solveable, (6 equations, 6 unknowns).
Kindly let me know how to write the code to solve it!
below is the code.

> restart:

interface(warnlevel=0):

with( plots ): with( plottools ):

with( DEtools ): with( PDEtools ):

g:=9.81:

x[0][2]:=0:

y[0][2]:=0:

#initialization function of each beam,in this function, we should give 3 variables, n=the number of this #beam, la=the length of this beam, ma=the mass of this beam

beams:=proc(n,la,ma)

global l,m,x,y:

l[n]:=la: #length of beam n

m[n]:=ma: #mass of beam n

#initialize the positions of each beam

x[n][2]:=0: #the x position of the end point of each beam

y[n][2]:=0: #the x position of the end point of each beam

x[n][1]:=0: #the y position of the center point of each beam

y[n][1]:=0: #the y position of the center point of each beam

end proc;

#define the forces that applied on the end of beam n

Force:=proc(n,fm,fn)

global fx,fy:

fx[n]:=fm:

fy[n]:=fn:

end proc:

Consider the task of sorting a list of complex floating-point numbers by magnitude.

**First Attempt**
The usual method to do this in Maple is with the

**sort** procedure. By passing a boolean-valued function that computes then compares the magnitudes of two complex numbers, we can sort the list. The following procedure shows how this is accomplished.

sort1 := proc(L)
return sort(L, proc(z1,z2) abs(z1) <= abs(z2) end proc);
end proc:

A disadvantage of this approach is that the absolute-value procedure is called twice every time a pair of numbers is compared. For a long list, the time spent in the absolute value routine dominates the computation time.

Is there a way to do a Symbolic Computations under Max-Plus Algebra? For instance, linear algebra computations under under Max-Plus Algebra. I would like to appreciate any comments on it.

I need to search for primitive polynomials of degree n over finite fields GF(q), where q is a PRIME POWER. I have major difficulties doing that, especially when I wish to fix one of the coefficients of the polynomial. On the other hand, polynomials over GF(p) pose no problem. Does anybody have any experience in this? PLEASE HELP!

Hi there everyone.
I have recently been trying to plot fourier transforms.
After taking the fourier transform of a particular function of the time domain for example a sinusoid. I get a sum of terms each with a translated Dirac function.
When I try and plot the series of Dirac functions it doesnt work in the sense that there is nothing displayed. The plot initializes but there are no Dirac functions.
Is there a way to plot Dirac Functions?
thanks

Dear Maple Users,
I have the following problem:
How to find the best-fitted function P(x,y) = ax+by+cxy+d to the following data:
i= 0, 1, 2, 3, 4 , 5 , 6 , 7
xi= 150, 152, 154, 156, 156, 157, 158, 158
yi= 134, 136, 106, 61, 90, 124, 41, 135
I know that in Maple functions that can solve interpolation and approximation problems with one variable are implemented. But how to solve the problem with two variables?
Best regards,
Wieslaw Kotarski
kotarski@gate.math.us.edu.pl

Hi to everybody
I'm a new user of Maple
How can I solve this system?
> eq1:=1039.44*Diff(Diff(ck(r),r),r)+80960*Diff(ck(r),r)*Diff(p(r),r)+80960*ck(r)*Diff(Diff(p(r),r),r)=0;
/ 2 \
| d | / d \ / d \
eq1 := 1039.44 |---- ck(r)| + 80960 |--- ck(r)| |--- p(r)|
| 2 | \ dr / \ dr /
\ dr /
/ 2 \
| d |
+ 80960 ck(r) |---- p(r)| = 0

hi everyone ,
please help me to sort out this problem ...
suppose i have a matrix
1 2 3
? 2 ?
2 1 6
where ? are unknown values .
what i want to do is to "Fill in the missing values by the respective attribute means. " How i will do it matlab?
please help me as soon as possible ...

I have been experimenting with overloading existing Maple functions by creating an appropriate package, as described in the "overload" help. It works fine in general.
I was hoping to change the definition of `abs` on type Vector to be sqrt(V . V) in this way, but the change does not work because `abs` already has an alternative definition on this datatype and it seems that built in functions take precedence over overload versions that one has added. Is there any way of achieving this?
Grateful for any pointers.