Hey all, I am new to the Maple software.

I have a question: how to find sum of all elements in a sequence?

For example , A:=[1,2,3,4]  the stupid way to find the sum is add them together by : A[1]+A[2]+A[3]+A[4]=10

Is there a better way to sum them together?

Hey all, I am new to the maple software. 

I have a question to create a binary sequence -consisting of m zeroes and n ones.  For example, m=n=2. then the number of combinations of the binary sequence is 6,nchoosek(4,2).  And the combinations could be {1,1,0,0}, {0,0,1,1}

{1,0,1,0} ,{0,1,0,1},{0,1,1,0} and {1,0,0,1}. How do I program the maple code ,that could print out all the combinations above?

hello dear freinds

im new comer in maple.

i want to find  particular solution of an ode by following code:

ode := diff(u[1](t), t, t)+u[1](t) = -(1/4)*a^3*cos(3*beta[0]+3*t)-(3/4)*a^3*cos(beta[0]+t)

m := combine(convert(particularsol(ode), trig))

but maple solution is : m := u[1](t) = (81/32)*a^3*cos(-3*beta[0]+t)-(81/16)*a^3*cos(3*beta[0]+t)-(3/8)*a^3*t*sin(beta[0]+t)+(3/16)*a^3*cos(-beta[0]+t)-(27/16)*a^3*cos(beta[0]+t)+(1/32)*a^3*cos(3*beta[0]+3*t)

but  particular solution is :

u[1](t) = -(3/8)*a^3*t*sin(beta[0]+t)+(1/32)*a^3*cos(3*beta[0]+3*t)

is there any idear for finding the solution?

thanks in advance

A 37 foot ladder is placed against a wall that is 9 feet away from its base. Will the top of the ladder reach a window ledge that is 35 feet above ground? Explain.

How do I get values out of solved augmented matric in 2D math input? I cannot simply assign the solutions to a value by x:=A[1,5}, because it's an augmented matrix.

The ability of Maple to solve differential equations is unsurpassed, but when the solutions appear in terms of Heun functions that result is disappointing because it is either difficult or impossible to convert those functions to other functions more commonly used and for which plots are readily generated.

Specifically, does any reader have a suggestion what to do with Heun C and Heun G functions?  In principle, they seem to be related to 1F1 and 2F1 hypergeometric functions, but the conversion seems not to succeed, and it is not obvious how to make it succeed.  In both cases of interest, the literature contains hints of solutions in other functions.

It seems that a solution of a differential equation in terms of Heun functions is not a solution at all.

Dear all,

I have one quetion. I have written maple code based on Dr.Robert J. Lopez's one.

Here is my code.

> q1 := diff(x(t),t) = -a*x(t)*y(t) + b*z(t);

q2 := diff(y(t),t) = -a*x(t)*y(t) + b*z(t);

q3 := diff(z(t),t) =  a*x(t)*y(t) - b*z(t);

q4:=k(t)=x(t)*y(t)/z(t);

ic := x(0)=10,y(0)=10,z(0)=0;

C := < I put Experimental data….>:

> SS := proc(a,b)

local F, V;

if not type([a,b],['numeric','numeric']) then return 'SS'(a,b);

elif a<0 or b<0 then return 1e100;

end if;

F := dsolve(eval({q1,q2,q3,q4,ic},{:-a=a,:-b=b}),[x(t),y(t),z(t),k(t)],numeric, output=Array([seq(k,k=0..N)]));

V := convert(Column(F[2,1],5),Vector);

Norm(V-C,2);

end proc:

 I just put one more equation -eqn4. it is not differential equation.

I currently have error massage from maple. Could you let me know how I solve the equations in this case?

Thanks.

Hello,

Since I was working in Matlab with Galerkin method which implies periodic boundary conditions I was wondering how to implement this in maple.

I tried this:

restart;

pde2 := diff(u(x, t), t)+3*(diff(u(x, t)^2, x))+diff(u(x, t),x$3) = 0

IBC := {u(0, t) = u(2, t), u(x, 0) = sech(50*(x-1/2))^2+2*sech(30*(x-1))^2, (D[1](u))(0, t) = (D[1](u))(2, t), (D[2](u))(0, t) = (D[2](u))(2, t)}

pds := pdsolve(pde2, IBC, numeric, time = t, range = 0 .. 2)

But it's telling me: 

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions can only contain derivatives which are normal to the boundary, got (D[2](u))(0, t)

So what's wrong?

following commands on my computer got an error.

restart;
with(LinearAlgebra):
A:=Matrix([[1,3],[2,5]],datatype=float):
b:=Vector([1,1],datatype=float):
LinearSolve(A,b,method=hybrid);
Error, (in SWcallhybrid[1]) param 4 should be an rtable

any suggestion is appreciated.

win7, 12.02

Hi,

I'm sorry, a few das ago I asked about it. Problamy my question was unclear.
I've decieded asking again.

I need a prcocedure (or contitional loops/sequences) which allow me to put values which are equivalent
to symbols into the array. The value are taken from eight diferent list. 

I enclose maple worksheet and screenshoot with explanation what I need.

I wil be appreciated for any help.

DA_-_help.mw

I have a new question here about using maple, welcome to answer

http://math.stackexchange.com/questions/1049655/how-to-calculate-a-in-this-description

Dear all, I have been trying to use Runge-Kutta method to plot an approximate solution with the following code. However, although I can get the numerical approximation the plot would not show.

h := .1;

x[0] := 0;

y[0] := 1;

xf := 3;

n := floor(xf/h)

f:= (x,y)->1/(3 y-x-2)

x := x[0]

y := y[0]

for i to n do

k1 := f(x, y);

k2 := f(x+(1/2)*h, y+(1/2)*h*k1);

k3 := f(x+(1/2)*h, y+(1/2)*h*k2);

k4 := f(x+h, h*k3+y);

k := (k1+2*k2+2*k3+k4)*(1/6);

y := h*k+y;

x := x+h

end do;

y[n]

data := [seq([x[n], y[n]], n = 0 .. 30)];

p[2] := plot(data, style = point, color = blue);

p[3] := plot(data, style = line, color = blue);
display(seq(p[n], n = 2 .. 3));

Does anyone in the community know example worksheets that analyze "billiards", viz. trajectory sets of free, specularly reflected particles constrained to bounce around in a particular 2d domain?

Are there any examples using the Bunimovich Stadium Billiard?

Thanks,

Bob Terry

declare(W(x, y), Z(x, y));

sys := [-A*kappa3-`&PartialD;`(`&PartialD;`(W(x, y))/`&PartialD;`(x))*(2*G-A)/`&PartialD;`(x)-2*G*(`&PartialD;`(`&PartialD;`(W(x, y))/`&PartialD;`(y))/`&PartialD;`(y)+`&PartialD;`(`&PartialD;`(Z(x, y))/`&PartialD;`(x))/`&PartialD;`(y))+A*`&PartialD;`(`&PartialD;`(Z(x, y))/`&PartialD;`(x))/`&PartialD;`(y) = 0, `&PartialD;`(`&PartialD;`(Z(x, y))/`&PartialD;`(y))*(A-4*G)/`&PartialD;`(y)+`&PartialD;`(`&PartialD;`(W(x, y))/`&PartialD;`(x))*(A-2*G)/`&PartialD;`(y)-2*G*`&PartialD;`(`&PartialD;`(Z(x, y))/`&PartialD;`(x))/`&PartialD;`(x) = 0];

I have this system of coupled PDE and I wish to solve it using Maple.

It gives me error of this kind:

pdsolve(sys, [[W(x, y)], [Z(x, y)]]);

Error, (in pdsolve/sys) found functions depending on different variables in the given DE system: [`&PartialD;`(x), `&PartialD;`(y)]

Thanks a lot for help