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

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.

Today science professionals in engineering software used to only work on the desktop and even just looking to download and use mobile apps math; but they are not able to design their own applications.Maplesoft to set the solution to it through its Maple package; software supports desktop and mobile; solves problems of analysis and calculation with Embedded Components. To show this we have taken the area of different mathematical topics; fixed horizontally to a certain range of parameters and not just a constant as it is customary to develop. This paper shows how the Embedded Components allow us to develop mathematics in all areas. Achieving build applications that are interactive in mobile devices such as tablets; which are used at any time. Maple gives us design according to our university or research need, based on contemporary and modern mathematics.With this method we encourage students, teachers and researchers to use graphics algorithms.

CSMP_PUCP_2014.pdf

Coloquio_PUCP.mw

Lenin Araujo Castillo

Physics Pure

Computer Science

Download FermiCreationAnnihilation.mw

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

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

I don't know what mean is T and V0x?

Please.