## Row sum...

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

## boundary condition for dependent variable...

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

## How to call an equation 25 times, then plot result...

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```

I'm study in Thailand's University. Now I would like to make seminar about Maple T.A. But I don't have knowledge about that. Thank you for all information.

## Pi as an axis label...

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

## ODE BVP...

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.

## help with numerical answer to maple expression...

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

## Stripping units...

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

## Liblip problem - Loss of Kernel connection...

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):

## Serious Help needed !!!.....

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]))

## Manipulate plot3d polygons...

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?