hey can someone plzz help me asap.
I need to write a procedure to determine the lcm and gcd of 2 integers. But only 1 will do bcz they are similar.
I have already written 2 procedures, one to determine the prime divisors returning the set which contains them, and the other determining the multiplicity of each divisor.
also, what if each integer has the same prime divisors with different multiplicities how do i fix this problem?

Hi,
I've this problem with maple11;
when I try to assigne the elements of any array (Tent in this case) using a procedure inside a for loop:
Tent:= array(1..7);
for i from 1 to 7 do
Tent[i]:= proc(v)
global i;
christoffel[i] + v^i;
end proc;
end do;
Maple dont understand that the 'i' inside the procedure is the same of the for loop. So after I've execute the for, i is actually 8, and if I write
Tent[3](p);
maple return
christoffel[8]+ p^8
instead I would
christoffel[3]+p^3
I've undersand that when I invoke Tent[3] (or Tent[2], etc) maple control the actual value of i (that is 8), so I think this is a problem of scope of the i variable.

Hi, I am using Maple 10 on Windows XP, and every time I try to open the Classic worksheet, the program opens, but whenever I try to run code, open an old worksheet, save a new worksheet, or pretty much do anything, the program crashes and I get the standard Windows message that "cwmaple.exe has encountered a problem and needs to close. We are sorry for the inconvenience. If you were in the middle of something, the information you were working on might be lost." Then it asks me to send the error report to Microsoft.
I have reinstalled several times using my original disk, contacted Maple support for help and reinstalled using the link they gave me, turned off my firewall, and nothing has helped. What else should I do? I have 1G of RAM and 1.67Ghz processor speed. Thanks so much, I would really appreciate any help as I need the classic worksheet to do the work for my final exam.

I'm trying to import an audio file, a two channel .wav file that I saved onto my computer, but every time I try to import it, I just get an array full of zeroes. Are there any restrictions no the types of wav files that Maple can import? Is there any way to fix this? I know the file is good.

Maple Professionals:
I imported a 331*17 matrix from a csv file with the help of the import assistant. Now I wish to take the summ of all rows so that I should get a 331*1 matrix. I attempted this with the following code:
`for i to 331 do DCdemand[i] := sum(Demand[i, j], j = 1 .. 17) end do;`

However, I then get the Error, "Matrix index out of range.
I would be extremely grateful if someone could help me on this issue. Thanks in advance, Frederik

I all, I have a system of three ODEs. I want to give an upper boundary condition for one of my state variables, B. Can I do this on dsolve? I am pretty new to maple, but it seems like all I can do is set a range for the independent variable (time, t in this case). Any ideas?
dsolve([Eq1,Eq2,Eq3,B(0)=100,W(0)=1200,N(0)=500],B(t),W(t),N(t),numeric)
can I write something like
dsolve([Eq1,Eq2,Eq3,B(0)=100,W(0)=1200,N(0)=500],B(t),W(t),N(t),**B=0..400**,numeric)
Thanks for your input!
Regards,
Oun

I am trying to plot the output of an equation.
I want to vary one parameter in the equation from 0 to 2.5 in increments of 0.1
I want to plot the output as a voltage
vs the change in current density J0.
Do I need to use an array? And how do I plot
the output( 2D)?
Here is the Maple input:

```
View 2470_CHE260_HW_MODEL.mw on MapleNet or Download 2470_CHE260_HW_MODEL.mw
```

How do I get the Pi symbol to be a label on an axis. I know how to do labels, my only problem is getting symbols to appear as a label.
Thanks in advance,
-Greg

There are two equations f'''+1/2*f*f''=0 and g'''+1/2*g*g''=0
boundary value
f(0)=0
g(0)=0
f'(0)=g'(0)(is not zero)
f''(0)=g''(0)
f'(4)=1
g'(-4)=0.5
Because of two coupled boundary value (f'(0)=g'(0, f''(0)=g''(0)), it
is difficult to solve.
Could you tell me how to solve this equation using matlab?
Thanks.

Hello, I have the following maple expression (sorry, it is somewhat lengthy)

e_x0 :=

4*ep^2*sin(m*Pi/b*y)*sin(n*Pi/c*z)*J0*(x1-x0)*(-cos(n*Pi/c*z1)+cos(n*Pi/c*z0))*(-cos(m*Pi/b*y1)+cos(m*Pi/b*

y0))/((2*m^2*Pi^4*c^2*n^2*b^2+n^4*Pi^4*b^4+8*Pi^4*m^4*c^4*Q^2+16*Pi^4*m^4*c^4*Q^4+16*k^4*b^4*c^4*Q^4+8*m^2*

Pi^2*c^4*k^2*b^2*Q^2-32*Pi^2*n^2*b^4*Q^4*k^2*c^2+8*n^2*Pi^2*b^4*k^2*c^2*Q^2+m^4*Pi^4*c^4-32*Pi^2*m^2*c^4*Q^

4*k^2*b^2+16*Pi^4*m^2*c^2*Q^2*n^2*b^2+8*Pi^4*n^4*b^4*Q^2+32*Pi^4*m^2*c^2*Q^4*n^2*b^2+16*Pi^4*n^4*b^4*Q^4)/b

^4/c^4/Q^4)^(1/2)/a/m/Pi^2/n;

and the input parameters are

Can someone please explain the rational behind why this works.
b := seq(x^2*Unit('m'), x = 1 .. 5);
1[m],4[m],9[m],16[m],25[m]
eval(b, Unit = 1);
1, 4, 9, 16, 25
I use this method but never understood why it works.
why does using Unit=1 strip the units?
Thanks,
DG

Hello all.
I am have recently gotten help in here, where it was suggested to me that i used Liblip to interpolate my 3d datasets. This has worked fine after i figured out how to make it work. I have however now encountered a problem with the LipIntSmoothLipschitz command, where on the latest dataset i am running, it makes a memory drain which ultimately after a long data accumulation, results in Kernel connection loss. I hope someone has an idea of what is wrong.
My loaded data is assigned like this:
First line: 2 (number of dimensions), 1852 (number of lines)
Second line and below: X Y Z (In first part sorted by increasin Y)
Here is my setup (Have made notes in, and at the bottom of worksheet):

Hi everyone - I have two problems Im really hoping someone can help me understand how to overcome (I have spent hours and hours going round and round in circles!)
Firstly the topic is square pyramid numbers i.e (1), (1+4), (1+4+9), (1+4+9+16), .. etc
Im trying to come up with a program that will see a procedure, nPyrLayers(), which says which "layer" of a pyramid a number is.. i.e 5 is the 2nd layer, 14 is the 3rd and then draw a pyramid.
I can build the pyramids with blocks (cubes) myself, I can find the nPyrLayers number, but I cant come up with a program that will draw ANY pyramid for any layer USING nPyrLayers... Is there a program I could write that could see "nPyrLayers" and then draw the pyramid with blocks itself?

I am trying to create a plot for two differential equations but I'm not getting any graphs with the commands below. I'm not sure what's causing the problem and I've tried everything I can think of. Please help if you can.
a := 2*ln(2)
b := 1/5*ln(2)
dose := proc (t) options operator, arrow; sum(2*Heaviside(t-6*n)-2*Heaviside(t-6*n-1/2), n = 0 .. 10) end proc
J:=DEplot([diff(x(t), t) = dose(t)-a*x, diff(y(t), t) = a*x-b*y], [x, y], 0 .. 50, {[0, 0, 0]}, stepsize = .5, scene = ([t, x]))
K:=DEplot([diff(x(t), t) = dose(t)-a*x, diff(y(t), t) = a*x-b*y], [x, y], 0 .. 50, {[0, 0, 0]}, stepsize = .5, scene = ([t, y]))

Suppose I generate a PLOT3D structure by
p:= plot3d(f(x,y),x=-1..1,y=-1..1,grid=[3,3]):
p2 := convert( p, POLYGONS );
Now I want to go through these polygons and throw away those that have a vertex far away from the origin (say more than 10 units), and then redisplay the trimmed structure.
Is there a snazzy way to do this?