AoA. How are you? I want to solve the following equation in one step and convenient way please help

psi(x):=Matrix(1,3,[x,x-1,x^(2)-x-1])

P := Matrix(3, 3, [1, 2, 1, 0, 1, 1, -1, 2, 3])

U := Matrix(3, 3, [a[1, 1], a[1, 2], a[1, 3], a[2, 1], a[2, 2], a[2, 3], a[3, 1], a[3, 2], a[3, 3]])

u := psi(x)^T*(P^2)^T*U*P^2*psi(t)-t*psi(x)^T*(P^2)^T*U*P^2*psi(1)-x*psi(1)^T*(P^2)^T*U*P^2*psi(t)+t*x*psi(1)^T*(P^2)^T*U*P^2*psi(1)

If you reduce my computational work or have any suggestion to solve please give.

PhD (Scholar)
Department of Mathematics

AoA. How are you? Hope you will be fine. I want write the following in Loop for any M

F[r,s]:=2^(k+2) m sqrt((c[r-1])/(c[s-1])) for r=2...(M+1), s=1..r-1 and (r+s)=odd;

F[r,s]:=0  otherwise

Novice question

I have two sets of inequalties:

s1:= (0<x, x<100)

s2:= (5<x, x<95)

I want each set to define a "range", so I can do set operations on the ranges. For example I want to determine that S2 is a subset of S1, S1 - S2 would give me (0<x, x<=5, 95<=x, x<100), S1 Union S2 = (0<x, x<100), etc

I assume these are simple things that can be done with Maple - but I am not sure how! Any suggestions?

I under stand how to do the calculus that's easy enough, but Maple is stil giving me a hard time.  The math par t is let the radius of the cylinder be r cms. 0<r<10 Then the height is 2 * Sqrt 100-r^2, the volume is V(r)=2PIr^2* sqrt 100-r^2 cm^3

Differentiante V(r) with respect to r we get V' (r) =2PIr^2(1/(2*sqrt(100-r^2)))(-2r) + (2PIsqrt(100-r^2))(2r)

V'(r)=(-2PIr^3 +4PIr(100-r^2))/sqrt (100-r^2)      V'(r)=2PIr(200-3r^2)/sqrt 100r^2         to find the critical point find the value of r when V(r) =0     2PI(200-3r^2)/sqrt 100-r^2 =0   2PIr(200-3r^2) = 0   r(200-3r^2)=0 

r=0 or 200-3r^2 =0     r=0 or r^2=200/3 r=sqrt(200/3) =10sqrt(2/3)   r is not in the interval so ignore it.  The critical point for 0<r<10 occurs at r = 10sqrt(2/3)       Since V'(r) >0 for 0<r<10sqroot(2/3) and V'(r) <0<10sqrt(2/3)<10

The dimensions are r=10sqrt(2/3) =8.16 cm and h=20/sqrt3 = 11.55 cm the volume is 4000PI/3*sqrt3=2418.4 cm^3

So how do I translate this into maple? Any help is greatly appreciated.

I am really baffled by this. I was starting to think Maple control system is good, only to find I get an error when I try to obtain the controllability Gammian for a state space system which is time varying.  Maple only seem to support a non-time varying A,B,C,D for this function.  Can someone please confirm this, or may be I am doing something wrong?

I was trying to verify my hand solution for a HW using Maple.

----------------------
restart:
alias(DS=DynamicSystems):
A:=Matrix( [[2,-exp(t)],[exp(-t),1]]);
B:=Matrix([[0],[1]]);
sys:=DS:-StateSpace(A,B);  #created OK
m:=DS:-Grammians(sys);

Error, (in DynamicSystems:-Grammians) symbolic value(s) in model: t, exp(t), exp(-t)
---------------------------------------------------

If the system is time invariant, then it works. Example from help:

-------------------------------------------------
restart:
alias(DS=DynamicSystems):
aSys := DS:-StateSpace( <<-5,3>|<3,-4>>, <<2,3>>, <<1,0>|<0,1>>, <<0,0>> ):
m:=DS:-Grammians(aSys);
---------------------------------

The main advantage of using computer algebra system for control system is to be able to do things using symbolic which is much more powerful than just numerical computation as with Matlab. There is nothing in the above linked to help page that says Maple does not support time varying state space for this function (unless I missed it).

Is it possible to add support for LTV in Maple 19?

ps. note that the ControllabilityMatrix works on LTV:

m:=DS:-ControllabilityMatrix(sys); #OK

But I wanted the Grammian.

Hello,

I have a question: Why gives Maple no result for:

expr2 := sin(k*(t-tau))*sin(k*tau);
Int(simplify(expr2, symbolic), tau = 0 .. t);

It is possible with Wolfram Alpha: http://www.wolframalpha.com/input/?i=Int%28sin%28sqrt%28k%2Fm%29*%28t-tau%29%29*sin%28sqrt%28k%2Fm%29*tau%29%2C+tau+%3D+0+..+t%29

Thanks!

Download mmgf_part1.mw

Any help???

Suppose I have installed maple on a directory in Windows 7 64 Bit. Now I use the command mode in windows, by typing maplew "D:\xx.mw" > xx.txt, it doesn't work. 

Excuse me, how to run maplew on command mode?

Hello,

I need your help to display the matrix generated by a lenear system.

The linear system is okay, well written, but  there is an error in generate matrix, I get the wrong matrix.

Thank you very much for your help.

test.mw

I came across this website www.fuelsonline.ca that has their price for furnace oil and Propane.  I thought it would be fairly easy using HTTP[Get] to fish out the prices, but it turns out those values on the page are images in of course a format (PNG) that cannot be read by Maple.  I don't know where that web page is pulling the prices from but it looks like web surgery with Maple will be unable to (at least not easily) wean the numbers out. 

If anyone has any ideas I would be interested in knowing.

Hi,

I have two all 1.-1  square matrices and want to test whether they are equivalent or not. Can anyone help me with details?Thanks in advance

Dear all

I need your help to create a band matrix of n block,   and each block of the matrix is given by:

A=[[1,2,3];[2,1,3]; [3,2,1]]:

Thanks for your help

I'm trying to find the partial sum of the function. Now I need to plot the first couple partial sums onto 1 graph.

I'm not really sure how to input the plot function. I was able to graph it by inputting each partial sum function but I would prefer an easier solution.

Hi every body:

i will earn the cofficients of q(T),q(T)^2 and (diff(q(T), T, T)) with commonds maple,how?

eq := 324.6463527*(diff(q(T), T, T))+4.012505275*10^11*q(T)+3.589858529*10^12*q(T)^2 = 0;

Hey guys,

I have a problem with this code:

d := Array(1 .. 15);

for i to s1 do

for c2 to 15 do

if `mod`(c2+(J+H)*(`mod`(1/(H+1), A[i])), A[i]) = 0 then d[c2] := d[c2]+l(A[i]) else d[c2] := d[c2]end if

end do

end do;

d;

The problem is that my code stops running when I encounter that `mod`(1/(H+1), A[i]) doesn`t have an modular inverse. What I would like to have is that the loop just continues and does nothing when my array when it encounters an error.

Hope you can help me!