Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi, 

I'm new to maple and have a problem when using the dsolve in a for loop. I get the error message Error, (in dsolve/numeric/process_input) the independent variable (t in phi(t)) cannot be a procedure name and I dont now why. 

The code design is 

for k to kk do
vars := x(t), y(t), phi(t);
ic := x(0) = xx, y(0) = yy, phi(0) = zz;
sys := diff(x(t), t) = y(t),
diff(y(t), t) = -2*gamma*y(t)-alpha*x(t)-beta*x(t)^3+F*cos(phi(t)),
diff(phi(t), t) = 1;
sol := dsolve({ic, sys}, {vars}, numeric, output = listprocedure, numeric = -1);
assign(sol);
x := x(t);
y := y(t);
phi := phi(t);
numpts := 50;
pts := seq([x(2*Pi*i), y(2*Pi*i)], i = 1 .. numpts);
plt(k) := plot([pts], view = [-1.5 .. 1.5, -.5 .. 1], color = black);
xx := 'xx'; yy := 'yy'; zz := 'zz';
xx := x(2*Pi*numpts);
yy := y(2*Pi*numpts);
zz := phi(2*Pi*numpts);
unassign(x, y, phi);
c := x(0) = xx, y(0) = yy, phi(0) = zz
end do

Let's say you have a m times n matrix and you want to ask Maple to computes all of its p minors (p smaller than minimum of m and n). Is there any prepared command for it at Maple? I can write a procedure myself, I just want to know if there is a command in Maple or not. The minor command in LinearAlgebra package is just computing determinant of a square matrix after removing a given row and a given column which is not what I'm talking about.

Hi all,

 

I am totally stuck in the question:

I can do it for seperate cases, but how to do it, with somekind of sequence

 

Integral for x= 0 to infinity of [2 sin(x/2)/x] ^ 2n 

this has to be equal to:

2n* pi (sum of [(-1)^j * (n-j)^(2n-1)] / [j! * (2n-j)!]

I showed it for j=0 and n=1

But how to do it for n =1,2,3 and 4

 

Dear all

I need a help to finish just two lines in my code:

1)How substitute condition in a system of equations

2) How solve the obtained system ( not a linear system)

3) If its possible to plot the set of points (x_i,y_j, u_{ij})

Many thanks for your help

codetosee.mw

hi

how i can pdsolve these partial differentialequations?

thanks

10.mw
 

restart; R := 1; UB := 1; lambda := 1; pe1 := 1; pe2 := 1; L := 1; Gr := 1; Br := 1; p := 1; LinearAlgebra:-HermitianTranspose(L) := 1; Nb := 1; Nt := 1; a := lambda*pe1*u(r, z)*(diff(sigma(r, z), z))-(diff(r*(diff(sigma(r, z), r)), r))/r-lambda^2*(diff(sigma(r, z), z, z))-Nt*((diff(r*(diff(theta(r, z), r)), r))/r+lambda^2*(diff(theta(r, z), z, z)))/Nb

u(r, z)*(diff(sigma(r, z), z))-(diff(sigma(r, z), r)+r*(diff(diff(sigma(r, z), r), r)))/r-(diff(diff(sigma(r, z), z), z))-(diff(theta(r, z), r)+r*(diff(diff(theta(r, z), r), r)))/r-(diff(diff(theta(r, z), z), z))

(1)

b := lambda*pe2*u(r, z)*(diff(theta(r, z), z))-(diff(r*(diff(theta(r, z), r)), r))/r-lambda^2*(diff(theta(r, z), z, z))-Nb*((diff(sigma(r, z), r))*(diff(theta(r, z), r))+(diff(sigma(r, z), z))*(diff(theta(r, z), z))*lambda^2)-Nt*((diff(theta(r, z), r))^2+lambda^2*(diff(theta(r, z), z))^2)

u(r, z)*(diff(theta(r, z), z))-(diff(theta(r, z), r)+r*(diff(diff(theta(r, z), r), r)))/r-(diff(diff(theta(r, z), z), z))-(diff(sigma(r, z), r))*(diff(theta(r, z), r))-(diff(sigma(r, z), z))*(diff(theta(r, z), z))-(diff(theta(r, z), r))^2-(diff(theta(r, z), z))^2

(2)

c := -p+(diff(r*(diff(u(r, z), r)), r))/r-L^2*(diff(r*(diff((diff(r*(diff(u(r, z), r)), r))/r, r)), r))/r+Gr*theta(r, z)+Br*sigma(r, z)

-1+(diff(u(r, z), r)+r*(diff(diff(u(r, z), r), r)))/r-(-(diff(u(r, z), r)+r*(diff(diff(u(r, z), r), r)))/r^2+(2*(diff(diff(u(r, z), r), r))+r*(diff(diff(diff(u(r, z), r), r), r)))/r+r*(2*(diff(u(r, z), r)+r*(diff(diff(u(r, z), r), r)))/r^3-2*(2*(diff(diff(u(r, z), r), r))+r*(diff(diff(diff(u(r, z), r), r), r)))/r^2+(3*(diff(diff(diff(u(r, z), r), r), r))+r*(diff(diff(diff(diff(u(r, z), r), r), r), r)))/r))/r+theta(r, z)+sigma(r, z)

(3)

bc := {sigma(R, z) = 0, sigma(r, 0) = 1, theta(R, z) = 0, theta(r, 0) = 1, u(R, z) = UB, (D[1](sigma))(0, z) = 0, (D[1](sigma))(r, LinearAlgebra:-HermitianTranspose(L)) = 0, (D[1](theta))(0, z) = 0, (D[1](theta))(r, LinearAlgebra:-HermitianTranspose(L)) = 0, (D[1](u))(0, z) = 0, (D[2](u))(R, z) = 0}

{sigma(1, z) = 0, sigma(r, 0) = 1, theta(1, z) = 0, theta(r, 0) = 1, u(1, z) = 1, (D[1](sigma))(0, z) = 0, (D[1](sigma))(r, 1) = 0, (D[1](theta))(0, z) = 0, (D[1](theta))(r, 1) = 0, (D[1](u))(0, z) = 0, (D[2](u))(1, z) = 0}

(4)

``


 

Download 10.mw

 

I have this situation:

We have sol[1] = "A=1", etc. How do I pick up just the solution values 1, 2 and 3 from sol? This is a generated system in a larger program, so I don't have access to the individual numbers coming out of a single step.
Thanks!
mapleatha

Good day,

I have a problem with inequalities plotting (trsh_con.mw). I need to plot them(eq5,eq7) in the same box. But I get the wrong solution. The plot not giving me correct region for 2 inequalities. I attach the mw file. Are there any options to overcome such problems?

Thanks in advance.

 

I need help solving Translation Step4 and Dilation Step5

 

Translation

Ta,b(x,y)=(x+a,y+b)

 

Dilation

Dk(x,y)=(kx,ky)

 

 

In the partial fraction decomposition Maple often changes the order of the terms.
For example,

How do we remedy this?
Thank you!

mapleatha

Dear Maple  users,

I am facing the memory allocation problem while calculating the inverse Laplace transformation of a big expression. I have many such terms (401) and I am importing them from a text file.

For an individual term itself, Maple consumes huge memory but at the end, it is able to calculate it in about an hour. I want to automate the process of reading and writing the terms one by one from the different text files. While doing so, maple eats a lot of memory and I see the error msg

"Error, (in expand/bigprod) Maple was unable to allocate enough memory to complete this computation.
Please see ?alloc "

Is there a way to do so? The maple file is attached.

InvLaplace_ex.mw

Thanks,

Dear all

I have invalid operator parameter name, i don't know where is the mistake.

Many thanks for your help

 

fdfd.mw

I'm producing multiple matrices via a for-loop (n = + 100) and I want to export each matrix into the same Excel-sheet, however I want to insert all matrices in cells beneath each other. I'm wondering if anyone knows how to do this? I've tried a loop variable but that doesn't work inside the "  " to specify the cells...

Note: Matrices are all 11x6. That means I want them inserted in the same sheet, in cells like this:
1st matrix: A1:F11
2nd matrix: A12:F23
3rd matrix: A24:F35
etc.

In other words: how to execute a loop-variable inside the " ", or to bypass it???
Please help! It's for my MSc Thesis:)

Helo dearz!

I want to find the inverse laplace transfrom (numerically) of the attached file of solution u and v for any values of M. I am waiting for your positive response. 

Help_v2.mw

how to download  wkptest from http://cpc.cs.qub.ac.uk/summaries/ADTY_v1_0.html and install on maple?

As a powerlifter, I constantly find myself calculating my total between my competition lifts, bench, squat, and deadlift. Following that, I always end up calculating my wilks score. The wilks score is a score used to compare lifters between weight classes (https://en.wikipedia.org/wiki/Wilks_Coefficient ), as comparing someone who weighs 59kg in competition like me, to someone who weighs 120kg+, the other end of the spectrum; obviously the 120kg+ lifter is going to massively out-lift me.

So I decided to program a wilks calculator for quick use, rather than needing to go search for one on the internet. For anyone curious about specific scoring, a score of 300+ is very strong for the average gym goer, and is about normal for a powerlifter. A score of 400+ makes you strong for a powerlifter, putting you in the top 250 powerlifters, while 500+ is the very top, as far as unequipped powerlifting goes, including the top 30. For anyone wondering, my best score at my best meet was 390, although given the progress I’ve made in the gym, should be above 400 by my next meet.

Hope you all enjoy!

 Find it here: https://maple.cloud#doc=5687076260413440&key=301A440EFD2C4EDD8480D60B5E3147BF40CA460F842942449C939AB8D2E7D679

First 894 895 896 897 898 899 900 Last Page 896 of 2223