Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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

I understand that matlab has function 'rfwrite' to generate touchstone file. This is a datafile used for high frequency modeling.

Is there any such package in Maple?

 

 

I´m new to Maple and I´m trying to find the minimum distance of an ECC from its parity check matrix or the generator matrix. Is this possible?

A fact that seems to be not documented. Probably it should be obvious.

 

Digits := d0
evalf[d](expr);

 

The toplevel float sub-expressions in expr  are computed with Digits=d0 but in procedures Digits is set to d.
Notice that the actual float parameters of the toplevel procedures are evalf-ed with Digits=d0.

 

restart;

g:=proc() convert(1/3., string) end:
h:=proc() 1/3. end:

Digits:=3;
evalf[10]([
  1/3. = h(),
  convert(1/3.,string) = g(),  
  fsolve(3*x=1) = add([1/3]),
  fsolve(x/3=1/9.) # 1/9. being at top level is evalf-ed with Digits=3
])

3

 

[.333 = .3333333333, ".333" = ".3333333333", .3333333333 = .3333333333, .3330000000]

(1)

 

k:=proc(x) convert(x,string) end:

kernelopts(floatPi);

true

(2)

4.+Pi;

7.14

(3)

evalf(k(1/3.+Pi));
# floatPi seems to be ignored inside actual parameters

".333+Pi"

(4)

evalf(k(4+evalf(Pi)));

"7.14"

(5)

evalf(k(4+Pi));
# 4 not being float (or "infected" by a float) is not evalf-ed

"4+Pi"

(6)

 

### (this is documented)

`evalf/h` := proc() 7.777 end:

evalf(h(1/3));

.333

(7)

evalf('h'(1/3));

7.777

(8)

 

 

Download digits.mw

(edited)

First 885 886 887 888 889 890 891 Last Page 887 of 2215