M:= Matrix ([[a,b], [c,d]])

If a.d-b.c=0 then print (Matrix(2)); else print (Matrix ([[(d/a.d-b.c), (-b/a.d-b.c)], [-c/a.d-b.c, a/a.d-b.c]]); end if;

I don't understand where I have gone wrong. If anybody could help I would be greatly appreciative

Thank you

Has anyone used Maple for Empirical Mode Decomposition ? Any tips, Maple docs ? I can't seem to find anything when searching Maple sites.

I have recently written a maple program to deconvolute gamma-ray spectra using the Richardson-Lucy algrithm. Although this method works well I would prefer to use a method based on the Maximum Entropy algorithm, and would like to know if anyone has tried to write a Maple program to deconvolute 1 dimensional data?

Hi.

I'm having a problem with this little guy "cos(Phi1(x)-psi)". When printed, it becomes "cos(-Phi1(x)+psi)", when Maple does not know what Phi1(x) or psi are. No matter what sign i set for psi, it will only display +psi when printed. Anyone know what causes this?

In ode solve command i generated a large array data. The output shows a large order matrix of this form

[110001x6 Matrix

Datatype:Anything

Storage:rectangular

order:Fortran_order].

I want to export this matrix into a notepad. Which can then be used for plotting in TecPlot. 

Looking for good response

I downloaded Maple 17 student eddition for linux and when I click on the bin file to install it I get an error that says "This file is of unknown type". It gives the option to find an aplication local or online and it says no aplications available or the "MIME" file type canot be installed on this system. Im not sure what to do, do I need somthing to unpack the file or is it just not aplicable on my linux operating system (version: Cinnamon version 1.8.8 LinuxMint: olivia(X86-64))?

In physics as in math, we often use delta to indicate a (small but finite) variation in a parameter (e.g. delta x). Sometimes Delta (capital delta) is used for a larger variation, and there are even constructs like delta Delta x (a small variation of a larger variation in x). (Don't laugh or frown, this is being done & makes sense in certain situations.)

How do I write this best in Maple, esp. when I use 2-d math and am working on a (live-) script that may be handed out to others incl. students,? I.e. the appearance matters. Just writing delta x is interpreted (in 2-d input) as delta*x which is not what is meant (and leads to funny effects after the first simplify). If I write deltax (one word) then delta does not typeset in its proper greek form, which becomes an issue if you have many of these. I can write delta(x) and that often works (as long as you don't try to differentiate by x), but the typeset equation does not look right either (delta is not a function of x).

Is there a god way of doing this, or am I hitting a wall here?

TIA,

Mac Dude

question 2 of my homework is actually abit more simple imo. i think i performed all the operations correctly but however my 3dplot shows just a flat plane. i dont think thats right. if someone could take a look at my code in reference to this question. it would be great.

Suppose that a baseball is hit from 5 feet off of the ground with an initial velocity of 40 feet per second at an angle of above the horizontal. Let the horizontal and vertical components of the position, velocity, and acceleration be denoted respectively by and . Assume that the only acceleration is due to gravity (i.e, ).

a) If , plot the parametric equations over the time period that the ball is in the air.

b) Clear in maple by typing theta:='theta'; Now, suppose that we write in Cartesian form by solving for t: t:=solve(rx=x,t); Typing in ry should now yield a function that involves only x and . Make a 3-dimensional plot of for , 

ax:=0;ay:=-32;
0
-32
vx:=int(ax,t)+C1;
C1
vy:=int(ax,t)+C2;
C2
v0:=40;
40
theta:=Pi/3;
1
- Pi
3
C1:=solve(subs(t=0,vx)=v0*cos(theta),C1);
20
C2:=solve(subs(t=0,vy)=v0*sin(theta),C2);
(1/2)
20 3
vx; vy;
20
(1/2)
20 3
rx:=int(vx,t)+C3;
20 t
ry:=int(vy,t)+C4;
(1/2)
20 3 t + 5
C3:=solve(subs(t=0,rx)=0,C3);
0
C4:=solve(subs(t=0,ry)=5,C4);
5

ry;
(1/2)
20 3 t + 5
plot([rx, ry, t = 0 .. infinity]);

Hi,

I have an ode like this

ODE:=(diff(T(x), x, x))+P*(S+a*(1-exp(-L*x))/L)*(diff(T(x), x))=0;

bcs:=T(0)=1,T(infinit)=0;

where P, S, L a, are all constants.

let assume that 

z=P/L^2*exp(-L*x);

subing z into the ode, we can have

ode1:=diff(T(z), z$2)+(1+z*a-P)*(diff(T(z), z)) = 0;

bcs1:=T(P/L^2)=1,T(0)=0;

Is it possible to find a closed form solution (T(x) in a compact form)?

I'm trying to display some flow lines for a vector field. The vector field I have is:

What I put into Maple is:

I've tried this with a number of points, but I always get an error message along the lines of:

"Error, (in Student:-VectorCalculus:-FlowLine) cannot determine if this expression is true or false: (2.*(1)(0)-4.)^2+(-2.+4.*(1)(0))^2 < 0.9e-9"

Could someone please shed some light as to what's going wrong for me? I would be very grateful for any help.

Thanks.

For the recursive equation

x_n − 10x_n−1 + 50x_n−2 = 0
n = 2, 3, ...
x₀ = 0
x₁ = 10

How do I calculate x_n in Maple? I'm wondering because n doesn't start at 0 but at 2, how will that make a difference?

Thanks in advance!

Hi, the display function in the simplex package gives us a nice breakdown of the form Ax<=b, but how do we extract the individual elements, i.e. the matrix A and vectors x and b? Normal indexing into the expression doesn't work, i.e. if for some set of constraints of the form 2x+3y<=21 we enter

res:=display(constraintSet);

this gives us a nice expression, but res[1] etc. doesn't work as we might hope (i.e. by returning the matrix A).

a := matrix(...);
b := Convert to hermit matrix from a;
Norm(MatrixMatrixMultiply(b, a), 2) = Norm(a,2)^2

i have tried, but not exactly equal
Norm(MatrixMatrixMultiply(b, a), 2) = 0.01
Norm(a,2)^2 = 0.1

Guess what matrix a is?

How to calculate colimit with only substitution and solve? Any simple example to show the steps.

run a command string in C# by calling maple

it can run in maple if copy into maple

however return input string was not in correct format

String commandstring = "restart;with(LinearAlgebra):with(ExcelTools): filename := "0257.HK";open3 := Import(cat(cat("C://Temp//HK//Transportation//",filename),".xls"), filename, "B2:B100");high3 := Import(cat(cat("C://Temp//HK//Transportation//",filename),".xls"), filename, "C2:C100");low3 := Import(cat(cat("C://Temp//HK//Transportation//",filename),".xls"), filename, "D2:D100");close3 := Import(cat(cat("C://Temp//HK//Transportation//",filename),".xls"), filename, "E2:E100");n := 30;Round := proc(x,n::integer:=1) parse~(sprintf~(cat("%.",n,"f"),x)); end proc: t:=1; gg :=Matrix(n+1,1); ggg :=Matrix(n+1,1); for k from 0 to n do InputMatrix3 := Matrix([[close3[t+1+k] , close3[t+k], close3[t+2+k]],[close3[t+k], close3[t+2+k],0],[close3[t+2+k],0 , 0]]): InputMatrix3b := Matrix([[close3[t+2+k], close3[t+1+k] , close3[t+3+k]],[close3[t+1+k] , close3[t+3+k],0],[close3[t+3+k],0 , 0]]): InputMatrix3c := Matrix([[close3[t+3+k] , close3[t+2+k], close3[t+4+k]],[close3[t+2+k], close3[t+4+k],0],[close3[t+4+k],0 , 0]]): Old_Asso_eigenvector := Eigenvectors(MatrixMatrixMultiply(Transpose(InputMatrix3), InputMatrix3)): Old_Asso_eigenvector2 := Eigenvectors(MatrixMatrixMultiply(Transpose(InputMatrix3b), InputMatrix3b)): Old_Asso_eigenvector3 := Eigenvectors(MatrixMatrixMultiply(Transpose(InputMatrix3c), InputMatrix3c)): gg[k+1,1] :=Old_Asso_eigenvector[2][1,1]; od;Round(Re(gg[1,1][1,1]));";