Hi,

Is there a way to export worksheet to pdf format but not in A4 size since some lines are going beyond the page?

I have a problem solving a system of PDEs.

The system of PDEs are

PDE01 := -(l^2+1^2)*(diff(v(l, t), t))+(l^2+1^2)*(diff(R(l, t), l, l))+4*l*(diff(R(l, t), l))+4*l*v(l, t)/(l^2+1^2)^(1/4)-6*R(l, t)/(l^2+1^2)+(l^2+1^2)^(1/2)*(-1.1+sqrt(.1))^2*sqrt(24)*u(l, t) = 0

PDE02 := diff(R(l, t), t) = v(l, t)

PDE03 := diff(u(l, t), t)-sqrt((1.1^2-1)/1.1^2)*(diff(u(l, t), l))-2*l*sqrt(1.1^2-1)*u(l, t)/(l^2+1^2) = 0

the initial condisions are

v(l, 0) = 0, R(l, 0) = 0, u(l, 0) = sqrt((l^2+1^2)^(1/2))*10^(-5)*exp(-(l-10)^2/.5^2)

and the BCs are

bdry00 := {((30^2+1^2)/30^2)^(1/4)*v(-30, t) = -((30^2+1^2)/30^2)^(1/2)*(D[1](R))(-30, t), ((30^2+1^2)/30^2)^(1/4)*v(30, t) = -((30^2+1^2)/30^2)^(1/2)*(D[1](R))(30, t), u(-30, t) = sqrt(30^2+1^2)*10^(-5)*exp(-40000), u(30, t) = sqrt(30^2+1^2)*10^(-5)*exp(-10000)}

to solve the system,

I enter

pde := pdsolve({PDE01, PDE02, PDE03}, {bdry00, init00}, time = t, numeric, range = -30 .. 30, timesstep = 1/60, spaceste = 1/254)

then, I failed to get the result constantly.

I tried several cases changing the initial conditions...

Can you let me know what I am doing wrong?

I have a problem with IsMatrixShape. I have in my part of formulation this matrix expression: Q^TI_bQ

While I_b is a symmetric matrix, this matrix expression is clearly symmetric. However, when I try to check this issue with IsMatrixShape command, it returns false. I am extremely confused. Can anyone help me? Thanks in advance.

prob1.mw    prob2.mw

Hi, I don't know why this happened...

If I write the command in new sheet, it runs correctly( prob2.mw).

Second, help me to evaluate N, please...  (prob1.mw)

Regards :-)

Another application for the study of rational numbers in operations, generating fraction, etc.

Numeros_Racionales.mw

(in spanish)

Atte.

L.AraujoC.

Download Question_algsubs_3.27.15.mw

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);

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