MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

hello . i have a Partial differential equation i need some help with 

 

This is how it goes 

Hi everyone.

I always used simple command to plot graphs with parametric function:

> plot([seq([i, t, t=0..11-i], i=1..10)], thickness=3, color=blue); 

And I always obtained this:

But yesterday I installed second update for Maple 18, and now I get only error:

"Error, (in plot) invalid input: ToInert expects 1 or 2 arguments, but received 0"

What command I should use now?

I'd like to implement this in Maple, but being a newbie having hard time getting the syntax right:



This is what I tried:

a:='a';
f:= i->product( `if`(evalb(i<>j),(a-z(j))/(z(i)-z(j)),NULL),j=1..3);

but it is not skipping the i=j case, since when I call it with f(1) I get division by zero.

f(1);
Error, (in product) numeric exception: division by zero

Tried also

f:= i->product(`if`(evalb(not(i=j)),(a-z(j))/(z(i)-z(j)),NULL),j=1..3);
and few other things. I think the reason it fails is deeper than my level of Maple undertstanding.

How would you do this in Maple?

hello, I have a list of numbers, for example [1.2,5.6,7.2,0.5,-0.25,-4,6]. I would like to find the position of it in the list where the number is closest to zero. May I know if there is a function for this please? Thanks.

Abel equation solution...

November 13 2014 Hamada 15

Hi,

I am interested to know if Abel equation of the form:

 

y'=a*y^3+(b*x+c)*y^2+d*y+(g*x+f)    where: a,b,c,d,g,f are constants. 

 

is belonging to any solvable classes? I noticed that most of the classes solved in the literature using invariant analysis has a rational coefficients in x and they do not have the full form as the one abov. Any suggestions for this ?

Regards

 

Hamada

 

AoA. How are you? Sir want to generate the attached square matrix for any value of M. Please help.

 

PhD (Scholar)
Department of Mathematics

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!


restart:

m:=2; k:=10.1; a:=k-m; b:=k+m-1;

m := 2

k := 10.1

a := 8.1

b := 11.1

z:=(k*m)/10^(0.1*15);

z := .6387800873

yo:=10^(0.1*10);

yo := 10.

N1:=evalf(sum(((((yo/z)^((b-a+2*p-1)/2))*GAMMA((1-(b-a+2*p))/2))/(p!*GAMMA(p-a+1)*GAMMA(1+((1-(b-a+2*p))/2)))),p=0..10));

Error, (in NumericRange) summand is singular in the interval of summation

N2:=evalf(sum(((((yo/z)^((b+a+2*p-1)/2))*GAMMA((1-(b+a+2*p))/2))/(p!*GAMMA(p+a+1)*GAMMA(1+((1-(b+a+2*p))/2)))),p=0..10));

N2 := -701860.9759

 

 

 

 

 

 

 


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

5 6 7 8 9 10 11 Last Page 7 of 1083