Hi, I have an equation that i need to integrate. Howver, due to the absolute equations I have in it, an analytical integration is not possible. How could i use a numerical method instead?

In matlab, instead of using "int", i use "quad". Is there such a thing in maple?

Please advice. Thanks

I have a problem with inputting a new package in my version of maple (15) for mac. The packes makes it possible to do calculations with cos/sin/tan in degrees instead of radians. I need to use this for my Engeneering studies in Denmark.

In pc you just place the package in the labrary subfolder under the Maple 15 folder, but i can't find that subfolder anywhere on my mac??
Where should i put the package to make it work?

I have 2 files; "Trig.ind" and "Trig.lib" - whichone to use?

Hi everyone,

I have a problem with DEplot.

Here is the code:

restart;with(DEtools), with(plots), with(linalg):

z:=exp(I*((t/(2*Pi))+(2*Pi*j/3)));

f1:=Re((I/(2*Pi))((1/(subs(j=1,z)-(x(t)+I*y(t))))+(1/(subs(j=2,z)-(x(t)+I*y(t))))+(1/(subs(j=3,z)-(x(t)+I*y(t))))));

g1:=-1*Im((I/(2*Pi))((1/(subs(j=1,z)-(x(t)+I*y(t))))+(1/(subs(j=2,z)-(x(t)+I*y(t))))+(1/(subs(j=3,z)-(x(t)+I*y(t))))));

sys1:=[diff(x(t),t)=f1,diff(y(t),t)=g1]:

Is there a command (possibly in a package) which can help in reading parameters from ini text files like

    [section1]

    option1=1.0

    option2=2.0

    [section2]

    option3=3.0

Thanks in advance!

A := Matrix([[5, 4, 5, 4, 8, 4, 6, 1, 4], [2, 8, 4, 7, 5, 1, 9, 6, 1], [8, 4, 9, 8, 8, 8, 6, 2, 8], [7, 1, 8, 6, 6, 5, 3, 7, 7], [8, 7, 9, 4, 7, 8, 2, 7, 4], [5, 1, 7, 8, 5, 7, 2, 2, 8], [7, 4, 1, 7, 7, 1, 7, 8, 6], [2, 9, 9, 4, 7, 5, 2, 6, 8]]);

Create the vector b that is column 7 from A and the matrix C that is made from columns 1 to 6 and 8 to 9 of A (in the same order as the columns of A

E := proc (phi) options operator, arrow; int(1/(1-(1-1/kappa^2)*sin(phi)^2)^(1/2), phi = 0 .. phi) end proc

E(pi/2)

Is there a way to get this function to work properly without having to simplify the integral manually ?

kappa is a constant (with a value of 12.4)

Thanks in advance!

There is a utilde on the Latin pallette, among others. How do I get xtilde as a variable?

Consider a simple polymer , which consist of long chains of simple molecules attached to each other, e.g:

....-CH₂-CH₂-CH₂-...

If the interaction energy between consecutive sections ca be ignored or is low enough, the chain can be bend around at varios places to form varios configurations. One can start by  fixing the first section to a particular co-ordinate and then add additional components by use of random walk with...

Dear everyone,

  I have a question about the memory requirement of generalized eigenvalue solver in Maple 10 and further versions.

Suppose I have a code like

***

with(LinearAlgebra):
Digits:=34;
n:=30;
A:=Matrix(n,n);
B:=Matrix(n,n);

for i from 1 to n do
  for j from 1 to n do
     if i=j then
       A[i,j]:=i*j;
       B[i,j]:=1;...

I want to create a simple executable bat fil.

1) When I click on it maple should launch in the background (I dont want to see maple open)
2) Then Maple should open a txt file located on the desktop (containing a 1), read the file
and save  a new txt file on the desktop containing the number in the previous fil +1.

How can this be done?

i study a simple differential equations as follow and get the solution,but i faced two warning which puzzle me:

> restart;
> eq2 := diff(x(t), t) = x(t), diff(y(t), t) = x(t)^2+y(t)^2-1;
   ans1 := dsolve({eq2, x(0) = 1, y(0) = -1}, numeric, [x(t), y(t)], range = -10 .. 10)       ;

Warning, cannot evaluate the solution further left of -1.1646698, probably a singularity

Just a small digression.  It would be nice if the numbers on the x-axis could be rotated or angled so it does not look like they collide with each other.  When the numbers reach the 100's or 1000's there will be congestion. 

This is perhaps a suggestion for a future Maple version as the current versions do not support number rotations on the axis.

Hi everyone,

I have this differential equations:

gamma1:=1;gamma2:=1;

l:=t->sqrt((x1(t)-x2(t))**2+(y1(t)-y2(t))**2);
2 2
l := t -> sqrt((x1(t) - x2(t)) + (y1(t) - y2(t)) )

DE:=[diff(x1(t),t)=(-1/(2*Pi))*gamma2*(y1(t)-y2(t))/l(t),diff(y1(t),t)=(-1/(2*Pi))*gamma1*(x1(t)-x2(t))/l(t),diff(x2(t),t)=(-1/(2*Pi))*gamma2*(y2(t)-y1(t))/l(t),diff(y2(t),t)=(-1/(2*Pi))*gamma1*(x2(t)-x1(t))/l(t)]:

I want to draw the phase portraits (x1(t...