Good afternoon sir.

I request your kind suggestion to the above cited query.

With thanks & Regards

M.Anand

Assistant Professor in Mathematics

SR International Institute of Technology,

Hyderabad, Andhra Pradesh, INDIA.

Hi,

I am doing calculation with binomial numbers with some parameters. Usually the result is correct, but this time a strange notation "table" appears. I do not know how does this happen?

Thank you!

Yan

Hi,

1-Triying to plot a function divided by its maximum value,sometimes it works with some parameters that means, the max.value of the plot is 1.

But when i change the data the max. value in the plot in graeter than 1 which is wrong!! should be 1.

dont know why??

2- Changing different data in the parameters, the programme takes long long time then i stop it?

please help me with these two problems.

restart:
>
------------------------- Defining the nature of the variables used ----------------------
assume(T,real):Digits:=25:n:=1:tau:=Pi:
theta:=0:phi:=0:
lambda:=n;Omega:=1:Gamma:=0.01:
--------------------- Input---------------------------------
1

J1

term1:=(exp((Gamma+I*d)*tau)-1)/(2*(Gamma+I*d)):
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):
term2:=(evalf(-0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-(Omega*Gamma)/(4*n))*Ak2(d),k=0..50))):
J1:=(term1+term2):
J1mod:=(Re(J1))^2+(Im(J1))^2:
###### J2#########################
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))-(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))-(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):

J2:=(evalf(-0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-Omega*Gamma/(4*n))*Ak2(d),k=0..100))):
######################

J2mod:=(Re(J2))^2+(Im(J2))^2:
J3 same as J1differ in sign
term1:=(exp((Gamma+I*d)*tau)-1)/(2*(Gamma+I*d)):
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):
term2:=(evalf(0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-Omega*Gamma/(4*n))*Ak2(d),k=0..100))):
J3:=term1+term2:
J3mod:=(Re(J3))^2+(Im(J3))^2:
J4 same as J2 but -0.25-->2

J4:=-2*J2:
######################

J4mod:=(Re(J4))^2+(Im(J4))^2:

calculate the spectrum

Spec:=d->(exp(-2*Gamma*tau)*(J1mod*cos(theta/2)^2+J2mod+J3mod*sin(theta/2)^2-0.5*Re(J3*J4*sin(theta)*exp(I*phi))+0.5*Re(J1*J4*sin(theta)*exp(-I*phi)))):

with(plots):

tit:=sprintf("l=%g,W=%g,G=%g",lambda,Omega,Gamma):
Smax1:=max(seq(evalf(Spec(d)),d=-100..100)):
plot(evalf(Spec(d)/Smax1),d=-15..15,axes=boxed,title=tit,color=black,font=[2,3,18],thickness=2,tickmarks=[3,3],titlefont=[SYMBOL,14],font=[1,1,18],linestyle=1);

I have a multivariate polynomial P (written with the order "tdeg")in the (x_i)_i ; the coefficients are functions of the (a_i)_i. I'd want the set of the monomials that appear in the development of P. For instance P=a_1xy^2z^3+a_2x^4yz^2+a_3 is associated to {xy^2z^3,x^4yz^2,1}. Does there exist a command to do that ?

In a second time, I'd want the coefficients of {xy^2z^3,x^4yz^2,1} in the same order. Does there exist a command to do that ?

Thanks in advance.

Suppose one sets a system of differental equations in vector form, say 2 ODE's, like this:

restart;
with(LinearAlgebra):
ode:=Vector([diff(x(t),t),diff(y(t),t)])=Vector([2*x(t)+y(t),3*y(t)-x(t)]);

Then to solve these, what would be an easy way to do it, without having to rewrite them again manually as a set, as what one would normally do. Clearly one needs to map dsolve, and also convert the vectors to a set somewhere? I am not able to get the syntax right.

Is there an easy way to automatically convert/rewrite the above to

ode2:= diff(x(t),t)=2*x(t)+y(t)  ,   diff(y(t),t)=3*y(t)-x(t);

so that I can just do

dsolve({ode2},{x(t),y(t)});

Or, a way to map dsolve directly into the first from as shown?( the Vector = Vector form).

What is the Maple command that allows this:

restart;
set_:={a=3,b=5,c=9}:
magic_command( set_ );  %this will cause the display to show as below

to display as

          a=3
          b=5
          c=9

(just for display purposes, so it is easier to see each element on its own line. I'll use this to format result from say dsolve and solve and such)

HI,

Here is my question,

I want to plot dependent variable or say function f(x) on x-axis and independent variable x on x-axis. But dont know how to do it. Please help me for this.

Regards

Sunit

Hello,

I look for solving a DAE system i obtain after having determined the equations of a mechanical system with kinematic closed loops.

For that, thanks to the partitioning method, i could transform my DAE in ODE system.

But now, i don't manage to solve my ODE system.

The first issue was the calculation of big matrix with trigonometric functions. With your help in the post "Resolution of a big product of matrix with trigonometric function" (http://www.mapleprimes.com/questions/200012-Rsolution-Of-A-Big-Product-Of-Matrix) , i could calculate the different matrix involved are calculated.

Now, the resolution with my differential system is very long and never finished.

@Carl Love 3670 gives me good advices. He asked me to try to avoid symbolic calculations of the matrix.

Do not use the output option for dsolve.

Then the following procedure evaluates the Matrix AA:

AA:= proc(S,t)
local Cu_inv:= eval(Cu, S(t))^(-1), tCu_inv:= eval(tCu, S(t))^(1);
     eval(Avv - Avu.Cu_inv.Cv + tCv.tCu_inv.Auu.Cu_inv.Cv - tCv.tCu_inv.Auv, S(t))
end proc:

           To use it, invoke AA(Sol, t) where t is an actual numeric value.

Unfortunately, i don't see how i can use this method for the moment? But, i think that there is a step i don't understand. For me, i have to make the calculation of the big matrix AA and QQ before solving my differentiel system because my differentiel system is composed with AA and QQ.

Here you can find the system i try to solve.

and here you can find the maple file without the steps leading to the setting of the equations

calcul_des_matrices_.mw

I am so thankful if there are any hints about this question:

http://math.stackexchange.com/q/523048/8581

Generally, can we check the existence of a valid limit by using epsilon-delta method in Maple?

Thanks

I thought to ask here on this, may be someone has this book.

According to Amazon, it shows this book

http://www.amazon.com/Introduction-Maple-Andre-HECK/dp/1461265053/ref=sr_1_27?s=books&ie=UTF8&qid=1383016657&sr=1-27

(rest of message was deleted by this amazing Maple prime software, will write it again)

The book above is supposed to new?

"Book Description
Publication Date: November 6, 2012 | ISBN-10: 1461265053 | ISBN-13: 978-1461265054 | Edition: 3rd ed. 2003"

but 3rd edition was published in 2003. Is this really a new book for Maple?

Find a loop do define these matrices, even for large values of n. Record the Matrix Mn 

all i,j=1,2,13,..,n for example 

N := `<|>`(`<,>`(1, 2), `<,>`(2, 2))

P := `<|>`(`<,>`(1, 2, 3), `<,>`(2, 2, 3), `<,>`(3, 3, 3))

Q := `<|>`(`<,>`(1, 2, 3, 4), `<,>`(2, 2, 3, 4), `<,>`(3, 3, 3, 4), `<,>`(4, 4, 4, 4))

Please help!!!

How do you type  x bar in Maple? I did not find !

after composite

test5 := subs(y=-q/p, subs(x=-q/p, x+y));
test6 := subs(p=x,subs(q=y,test5));

how to recover test6 := -2*y/x back to x+y with y/x or x/y;

hi,

How can i solve Volterra integral equation in maple 13?

please help me.

thanks and nice job.

Hi there,

How can I get gridlines in my plots in Maple TA? Se eg. the source of the following question:

Thanks alot in advance!

mode=Multiple Choice@
name=Grafen for den afledede@
comment=@
editing=useHTML@
solution=@
algorithm=$r = rint(3);
$f1 = rint(3);
$f2 = rint(3);
condition: ne($r,$f1);
condition: ne($r,$f2);
condition: ne($f1,$f2);
$rf = switch($r,x^2,cos(x),x);
$ff1 = switch($f1,x^2,cos(x),x);
$ff2 = switch($f2,x^2,cos(x),x);
$fm = mathml("$rf");
$opgaveplot = plotmaple("plot($rf,x=-10..10), plotdevice='gif', plotoptions='height=250,
width=250'");
$rplot=plotmaple("plot(diff($rf,x),x=-10..10), plotdevice='gif', plotoptions='height=250,
width=250'");
$fplot1=plotmaple("plot(diff($ff1,x),x=-10..10), plotdevice='gif', plotoptions='height=250,
width=250'");
$fplot2=plotmaple("plot(diff($ff2,x),x=-10..10), plotdevice='gif', plotoptions='height=250,
width=250'");@
uid=61d45fb9-b1cc-4e0b-a986-9bbb24985734@
privacy=10@
allowRepublish=false@
description=@
difficulty=0.0@
modifiedBy=669aa1fb-4bde-408c-ae4d-6fa47454164c@
school=d1f34b37-d475-429f-b82f-5a0479c33756@
attributeAuthor=true@
numberOfAttempts=1@
numberOfAttemptsLeft=1@
numberOfTryAnother=0@
numberOfTryAnotherLeft=0@
question=<p>
Betragt funktion f(x) = $fm og dens graf herunder.&nbsp;</p>
<p style="text-align: center; ">
$opgaveplot</p>
<p style="text-align: center; ">
&nbsp;</p>
<p>
Angiv hvilken af følgende figurer repræsenterer grafen for funktionen?</p>
@
answer=1@
choice.1=$rplot@
choice.2=$fplot1@
choice.3=$fplot2@
choice.4=Ingen af dem@
choice.5=Der er flere rigtige@
fixed=@