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!

Hello,

I am quite new on Maple and Maplesim. A co-worker and I are evaluating a demo version (Maple 18 with Maplesim 7) for a possible purchase in 2015. For the project my collegue and I are working on, we need to detect when a signal becomes zero in order to switch between two options. We were trying to do this using the "Change" or "Boolean Change" components (the difference between those two is not really clear to us). Correct me if I am wrong, but those components should detect when a boolean signal has a change in value (both rising or decreasing) and in that case the output should become TRUE, while otherwise is FALSE.

The problem is that the component work if a "Terminate" is connected, but the output is always FALSE in any other case.

At the bottom we are attaching a file with two very easy examples.

In the first one the input of the "Boolean Change" component is a "Boolean Pulse". We compared also the output of the "Change", "Rising Edge", and "Edge with Initial Value" components. We were expecting istantaneous pulse with frequency equal to that of the boolean pulse, or its half for the Edge block. On the contrary the output of all the components stays zero (thus FALSE) for the whole simulation.

The second example use the output of a "Greater Threshold" component to enter the "Boolean Change", whose output is used to trigger a "Switch", whose output is in turn used as an input for the "Greater Threshold" component. The expected result is to have the "Switch" output to alternate the two inputs, instead it keeps the second one (thus the input for the switch is always FALSE).
The feature we need is more similar to this second output, but we don't understand what we are doing wrong or if there is another way to do it.
Thank you in advance
Irene

test_change.msim