Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I'm attempting to plot the Nyquist plot for a complex system such that Maple cannot determine the frequency limits automatically, and suggested I use 'range' to specify the frequency limits. How do you do this? For example if you have a system G=1/(s+1), how do you plot the Nyquist plot in the frequency range 0 to 1rad/s?  Thanks

Helo friends. Hope you will be fine. I need the command for take the conjugate of exp(I*x). I need the result exp(-I*x), x treated as real number.

 

PhD (Scholar)
Department of Mathematics

In a plot command, how do I indicate not to show the legend ?

(I know about unchecking "show legend" with a right click, and I have searched mapleprimes.)

 

Thanks, all !

Write a Maple code that performs the Gaussian elimination for an nxn matrix, converting it to an upper triangular matrix. 

(Hint: you will need to use three for .. do loops.)

How to add values in vector? This is what I am trying in my code below. I know there are more than one values satifying my conditon.

What is the best way constructing such vectors. If there are other options, please let me know. Thanks.

Write a recursive Maple procedure, called “decToBin”, that converts an integer from decimal (base 10) to binary (base 2). Ensure that only integers are allowed as arguments to the procedure. Passing a negative integer should result in no output. Test your procedure by outputting the result of the following:
decToBin(10)
decToBin(163)
decToBin(213)
decToBin(-1)
decToBin(“A”) # should result in an “invalid input” error

the calculation is like the following command, the result in the picture

restart;
with(VectorCalculus);
SetCoordinates(spherical[r, theta, phi]);
Fv := rho*VectorField(`<,>`(v[r](r, theta, phi), v[theta](r, theta, phi), v[phi](r, theta, phi)));
Divergence();
Divergence(Fv);

divergence

 

1) when the Divergence act on the Fv, then it will be expanded, which is lengthy and not like most book's formulation , especially when I want to continue for a Conversation law like in fluid mechanics, this will be too long and a messy for later check.

could there be a way to not expand this result, just as the eq(3) like.

2) when I want to calculate the Divergence of Fv, I must construct a VectorField at first, but this is in components way, is there a quick way for Vector Field Function

 

 A'(t)=20+0.03*A(t)-0.00005*A(t)*Q(t)-0.0002*A(t)*Q(t)-0.0743*A(t)-12,

B'(t)=0.03*B(t)+0.00005*A(t)*Q(t)+0.0002*A(t)*Q(t)-0.0743*B(t)-0.01072*B(t),

C'(t)=0.01072*B(t)-(0.0743+0.0123)*C(t),

P'(t)=15+0.03*P(t)-0.0001*P(t)*B(t)-0.0743*A(t)-10,

Q'(t)=0.03*Q(t)+0.0001*P(t)*B(t)-0.0743*Q(t)-0.25*Q(t),

R'(t)=0.25*Q(t)-(0.0743+0.0123)*R(t),

X'(t)=25-0.0003*X(t)*B(t)-0.0743*X(t)-10,

Y'(t)=0.0003*X(t)*B(t)-0.0743*Y(t)-0.3*Y(t),

Z'(t)=0.3*Y(t)-(0.0743+0.0123)*Z(t).

Initial conditions

A(0)=150,B(0)=30, C(0)=10,P(0)=130,Q(0)=45,R(0)=10,X(0)=85,Y(0)=30, Z(0)=20.

Can I plot the graph between {A(t),B(t),C(t),P(t),Q(t),R(t),X(t),Y(t),Z(t)}and time t{0,100}.

Hi all,

I've been working on this problem for a while and only need to figure out why the function vector is returning zeros every time I call the procedure. The problem statement is in the file I attached to this.

If any of you could help me, I would be a happy guy.

 

File: 1.2.mw

Just as the title says: How do I avoid Maple Prompt to appear automatically?

When writing a code in Maple Prompt and hitting enter, a new [> is created and I can't delete it. Which is unfortunate, for most of the time I want to write some TEXT afterwards, not more code.

Sometimes it doesn't happen, but I can't find the logic in it. Any ideas? Much obliged!

For an example, see this post:

http://www.mapleprimes.com/questions/203735-Print-To-Pdf-Results-In-Blank-Space

2nd picture. I asked it as a BTW-question.

 

Thanks in advance

hi friends

i have a problem in maple with an error

  

Read'dsnumsort.map';
dsnumsort(Ns(0),[x,y,z]);
for i from 0 to 1000 do;
T:=i/25;
NsT:=Ns(T):
X[i]:=rhs(NsT[C1]); Vx[i]:=rhs(NsT[V1]);
Y[i]:=rhs(NsT[C2]); Vy[i]:=rhs(NsT[V2]);
Z[i]:=rhs(NsT[C3]); Vz[i]:=rhs(NsT[V3]);
KepVec[i]:=convert(crossprod([X[i],Y[i], Z[i]],[Vx[i],Vy[i], Vz[i]]),list);
KepAbs[i]:=norm(KepVec[i],2);
od:

pacecurve({[seq(X[i],Y[i],Z[i],i=0..1000)],[[-1/(2),0,0],[1/(2),0,0]]},labels=['x','y','z']);

but i see:

Error, invalid input: seq expects 2 arguments, but received 4

can you helpe me?

Thanks

Hello!
Please help solve this equation!

restart;

Dd:=4.9*10^(-10):
n:=4.5:
m:=3.5:
ny:=20:
B:=1.47*10^(-11):
H1:=70:
H2:=20:
H:=100:
b1:=40:
b2:=80:
delt_y:=H/ny:
A_max:=b1*H1+b2*H2+10*b1:
M:=1000:




b:=y->piecewise(`and`(y>=0,y<=H1),b1,`and`(y>H1,y<H1+H2),b2,`and`(y>=H1+H2,y<=H),b1):
y:=0: eq:=0:

for j from 1 to ny do
y:=(j-1)*delt_y+delt_y/2;

eq:=evalf(eq+b(y)*sign(y-y0)*abs((y-y0))^(1/n)*delt_y);
end do;


#eq=0, y0-?
sol1:=fsolve({eq=0},{y0});

Hi

I'm trying to use the command: DocumentTools[Retrieve](path,label), but i recieve the message "unable to retrieve label reference" no mather what i do. When i use the GUI: Insert - reference and get the value it works fine, but in that way I can't assign the value to a variable in the document.

here is an example:

The document i want to retrieve a value from:

restart;

a := 5;

#a gets the label (1)

 

The document i want to retrieve the value to:

restart;
with(DocumentTools);

path := FileTools:-JoinPath(["test.mw"], base = worksheetdir); = "C:\Users\Nicolai\CloudStation\Nicolai\Skole stof\Sem 7\test.mw"

 

Retrieve(path, "());
Retrieve(path, "(1)");
Retrieve(path, "a");

#None of the above works

 

Is it even possible to do such thing?

with best regards

Nicolai

with(Physics) :
t:=Intc(Dirac(k1+k2+k3)*phi(k1)*phi(k2)*phi(k3),k1,k2,k3) ;


# how to force
Fundiff(t,phi(-k));

# to find
3*Intc(Dirac(k2+k3-k)*phi(k2)*phi(k3),k2,k3) ;

 

hi friends

i have a problem in maple with an error

 restart:
 with(plots):with(linalg):
 r:=sqrt((xi-x)^2+y^2+z^2):
 sigma:=m/l:#linear mass density
 U:=Int(G*sigma/r,xi=-L/2..L/2);
 U:=normal(value(U));
 Us:=subs(G=1,m=1,L=1,z=1.5,U);
 plot3d(Us,x=-2..2,y=-2..2,axes=boxed,style=hidden,color=red,orientation=[-80,-130],numpoints=35^2,labels=['x','y','U']);

after this i see :

 Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct.

can you help me?

Thanks for your time

Best regards

First 1287 1288 1289 1290 1291 1292 1293 Last Page 1289 of 2256