Advanced Search

Maple 13 Questions and Posts

These are Posts and Questions associated with the product, Maple 13
View unanswered posts
Maple 13 Questions and Posts Feed

second question about a simple two column matrix...

Asked by: Matt C Anderson 180
matrix export
February 06 2022
0  1

Thank you very much, Maple Community !

I have another (hopefully) easy question about a data set

seematrix_start_2_trim.xlsx

matrix_start_2_trim.pdf

I want to imitate Online Encyclopedia of Integer Sequences .b file

reference oeis.ort/A40

Cheers,

Matt

how to make a data file...

Asked by: Matt C Anderson 180
matrix export save
February 04 2022
0  3

Hi again,

Are there any Maple experts out there that can help me save a matrix full of data?

two_column_matrix.pdf

Regards,

Matt

Problem solving nonlinear equations...

Asked by: Delali Accolley 20
optimization rootfinding fsolve
January 26 2022
2  10

Hi, I would be grateful if anyone could help find inconsistency in a system of equations I have been solving. I can easily solve the system (that I have called ‘EQs’) when I set the parameter ‘tau\_l’ equal to .06, but for a different value, it will not solve. The same thing happens whenever I change the value of some other parameters like ‘l3’. Regards

ACCOLLEY_Delali_-_Demographics_3.mw

Numbering in set...

Asked by: Maria2212 25
union set
September 27 2021
1  2

При использовании «объединения» нумерация не сохраняется. Непонятно, по каким алгоритмам «объединение» применяется к тому или иному месту. Почему sistema [1] не равно GRAN [1]? Как я могу это исправить?

PS Прошу прощения за грамматические ошибки, английский плохо знаю

1. mws
 

> 

рестарт;
> 

 

> 

 

> 

B:=matrix([[1,        0,        1,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0],
[-1,        -1,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0],
[0,        0,        -1,        0,        -1,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0],
[0,        1,        0,        1,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0],
[0,        0,        0,        0,        1,        1,        0,        1,        0,        0,        0,        0,        0,        0,        0,        0],
[0,        0,        0,        -1,        0,        -1,        -1,        0,        0,        0,        0,        0,        0,        0,        0,        0],
[0,        0,        0,        0,        0,        0,        1,        0,        1,        0,        0,        0,        0,        0,        0,        0],
[0,        0,        0,        0,        0,        0,        0,        -1,        0,        -1,        0,        0,        0,        0,        0,        0],
[0,        0,        0,        0,        0,        0,        0,        0,        -1,        0,        -1,        -1,        0,        0,        0,        0],
[0,        0,        0,        0,        0,        0,        0,        0,        0,        1,        1,        0,        1,        0,        0,        0],
[0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        1,        0,        1,        0,        0],
[0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        -1,        0,        -1,        0],
[0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        -1,        0,        -1],
[0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        0,        1,        1]
]):
> 

E:=16:V:=14;

V: = 14

 (1)
> 

G:=2*E; êîëè÷åñòâî ãðàíè÷íûõ óñëîâèé

G: = 32

 (2)
> 

l:=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]:
> 

 

> 

for i to E do
De[i]:=-diff(phi[i](s[i]),s[i]$2)=He*phi[i](s[i]):  Äó âèäà -(diff(phi[i](s[i]), `$`(s[i], 2))) = H[e]*phi[i](s[i])
od:
> 

for i to E do
FSR[i]:=subs(s[i]=x,dsolve(De[i],phi[i](s[i]),output=basis)):  Îáùèå ðåøåíèÿ
od:
> 

for i to E do
OBCHEE_RECHIE[i]:=unapply(C[i+i-1]*FSR[i][2]+C[i+i]*FSR[i][1],x): od:
> 

 

> 

for i to E do
PROIZ:=diff(OBCHEE_RECHIE[i](x),x):
PROIZ_OB_RECH[i]:=unapply(PROIZ,x):od:
> 

 

×òåíèå ìàòðèöû èíöèäåíòíîñòè

> 

i0:=1:i1:=1:z:=1:y:=1:
> 

for i from 1 to V do
 y:=1: z:=1:
  for j from 1 to E do
   if B[i,j]=1 then
    x[z]:=j: z:=z+1:
   fi:

   if B[i,j]=-1 then
    x[y]:=j:
    y:=y+1:
   fi:

  od;

  if z>=2 then
   v0[i0]:=[seq(x[s],s=1..z-1)]:
   i0:=i0+1:
  fi:

  if y>=2 then
   v1[i1]:=[seq(x[s],s=1..y-1)]:
   i1:=i1+1:
  fi:
od:
> 

v0:=convert(v0,list):
> 

v1:=convert(v1,list):
> 

 

Ïîäñ÷åò êîëè÷åñòâà ãðàíè÷íûõ óñëîâèé òèïà f(0)

> 

number0:=0:
> 

for i to i0-1 do
if nops(v0[i])=2 then number0:=number0+2; fi:
if nops(v0[i])=3 then number0:=number0+4; fi:
od:
> 

number0:=number0/2; #êîëè÷åñòâî ãðàíè÷íûõ óñëîâèé òèïà f(0)

number0 := 9

 (3)
> 

 

Ïîäñ÷åò êîëè÷åñòâà ãðàíè÷íûõ óñëîâèé òèïà f(1)

> 

number1:=0:
> 

for i to i1-1 do
if nops(v1[i])=2 then number1:=number1+2; fi:
if nops(v1[i])=3 then number1:=number1+4; fi:
od:
> 

number1:=number1/2; #êîëè÷åñòâî ãðàíè÷íûõ óñëîâèé òèïà f(1)

number1 := 9

 (4)
> 

 

Êîëè÷åñòâî ãðàíè÷íûõ óñëîâèé ñ ïðîèçâîäíûìè

> 

GP:=i0-1+i1-1;

GP := 14

 (5)
> 

 

ÍÓËÅÂÛÅ ÓÑËÎÂÈß

> 

i:=1:
> 

while i<number0+1 do
for j to i0-1 do
if nops(v0[j])=2 then
GRAN[i]:=OBCHEE_RECHIE[v0[j,1]](0)-OBCHEE_RECHIE[v0[j,2]](0)=0;
i:=i+1;
fi;
if nops(v0[j])=3 then
GRAN[i]:=OBCHEE_RECHIE[v0[j,1]](0)-OBCHEE_RECHIE[v0[j,2]](0)=0;
GRAN[i+1]:=OBCHEE_RECHIE[v0[j,2]](0)-OBCHEE_RECHIE[v0[j,3]](0)=0;
i:=i+2;
fi;
od; od;
> 

 

ÅÄÈÍÈ×ÍÛÅ ÓÑËÎÂÈß (ÈËÈ îò l)

> 

i:=number0+1:
> 

while i<number0+number1+1 do
for j to i1-1 do
if nops(v1[j])=2 then
GRAN[i]:=OBCHEE_RECHIE[v1[j,1]](l[v1[j,1]])-OBCHEE_RECHIE[v1[j,2]](l[v1[j,2]])=0;
i:=i+1;
fi;
if nops(v1[j])=3 then
GRAN[i]:=OBCHEE_RECHIE[v1[j,1]](l[v1[j,1]])-OBCHEE_RECHIE[v1[j,2]](l[v1[j,2]])=0;
GRAN[i+1]:=OBCHEE_RECHIE[v1[j,2]](l[v1[j,2]])-OBCHEE_RECHIE[v1[j,3]](l[v1[j,3]])=0;
i:=i+2;
fi;
od; od;
> 

 

ÍÓËÅÂÛÅ ÏÐÎÈÇÂÎÄÍÛÅ

> 

i:=number0+number1+1:
> 

while i<number0+number1+nops(v0) do
 for j to nops(v0) do
  for y from 1 to nops(v0[j]) do
   x[y]:=PROIZ_OB_RECH[v0[j,y]](0);
  od;
  m:='m';
  GRAN[i]:=sum(x[m],m=1..nops(v0[j]));
  i:=i+1;
 od;
od;
> 

 

ÍÅÍÓËÅÂÛÅ ÏÐÎÈÇÂÎÄÍÛÅ

> 

i: = число0 + число1 + nops (v0) +1:
> 

в то время как i <G + 1 сделать
 для j до nops (v1) сделать
  для y от 1 до nops (v1 [j]) сделать
   x [y]: = PROIZ_OB_RECH [v1 [j, y]] (l [v1 [j, y]]);
  od;
  m: = 'm';
  ГРАН [i]: = сумма (x [m], m = 1..nops (v1 [j]));
  я: = я + 1;
 od;
od;
> 

 

> 

 

> 

 

новогодний

> 

sistema: = {}:

новенькие

> 

время: = {}:

> 

 

> 

для я от 2 до G делаю

> 

  sistema: = sistema union {GRAN [i]}:
> 

od:
> 

для я от 2 до G делаю

> 

  peremenie: = peremenie union {C [i]}:
> 

od:
> 

 

> 

коэфф: = решить (система, изменение):
> 

коэфф:
> 

 

> 

система [1];

C[2]*He^(1/2)+C[6]*He^(1/2)

 (6)
> 

ГРАН [1];

C[1]-C[5] = 0

 (7)
> 

 

> 

urav: = subs (koeff, GRAN [1]):
> 

uravnenie: = collect (numer (lhs (normal (urav))), {He, C [1]}):
> 

uravn: = remove (has,%, C [1]):
> 

uravnenie_sobs_chisla: = evalf (algsubs (He = x, uravn)):
> 

ВЕСЕЛЫЕ: = неприменить (uravnenie_sobs_chisla, x):
> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 

> 

 


 

Скачать 1.mws

 

 

I cannot use series function. How can I solve this...

Asked by: Tolibjon1996 5
syntax series
September 27 2021
1  1

Good day everyone. I need to calculate series of the function given above. But I can not find result. It writes "Error, (in series/function) improper op or subscript selector". How can I fix this error? Please help me.

The detailed file has been attached above.
 

> 

Ek4:=sqrt((Ea^2*Eb^2-Eab^2*X5^2)*Ds);

((((1/2)*k((1/2)*k^2/m2+X2)^2/m1((1/2)*k^2/m2+X2)+X1((1/2)*k^2/m2+X2))*((1/2)*k((1/2)*k^2/m2+X1)^2/m2((1/2)*k^2/m2+X1)+X2((1/2)*k^2/m2+X1))-((1/2)*k((1/2)*k^2/m2+X1)^2/m1((1/2)*k^2/m2+X1)+X2((1/2)*k^2/m2+X1))*X5^2)*(((1/2)*k((1/2)*k^2/m2+X2)^2/m1((1/2)*k^2/m2+X2)+X1((1/2)*k^2/m2+X2)-(1/2)*k((1/2)*k^2/m2+X1)^2/m2((1/2)*k^2/m2+X1)-X2((1/2)*k^2/m2+X1))^2+4*((1/2)*k((1/2)*k^2/m2+X1)^2/m1((1/2)*k^2/m2+X1)+X2((1/2)*k^2/m2+X1))^(1/2)*X5^2))^(1/2)

 (5)
> 

series(Ek4,k=0,4);

Error, (in series/function) improper op or subscript selector

  
> 

series(sqrt(Ek4),k=0,4);

Error, (in series/function) improper op or subscript selector

  


 

Download function.mwfunction.mw

prime 2 tuple with difference of 36...

Asked by: Matt C Anderson 180
export filetools incomplete-question
September 02 2021
1  1

Maybe some expert can help me with this.

I am trying to improve oeis.org/A156104. It is pairs of prime numbers p and p+36.
According to the k-tuple conjecture, this should be potentially an infinite list. I want to make the database .b file larger - from 1,000 to 10,000.

The virtual clipboard was too large to copy to notepad.  I could not figure out how to write this to a file.  I tried fopen() and fclose().

See my efforts

pairs_p_and_p+36_2.mw

pairs_p_and_p+36_2.pdf

Let me know.

Matt

more proper divisors from Matt...

Asked by: Matt C Anderson 180
factorization prime number-theory
June 11 2021
1  4

Hi all again,

Tried this.

some_proper_divisor_examples_8.pdf

some_proper_divisor_examples_8.mw

 

Anyone know a theorem related to this?

My guess is (p1 -1) * (p2 -1) * ... * (pb -1)  .

For p1, p2, ... pb all prime numbers.

Regards,

Matt

Proper Divisors Procedure in Maple code

by: Matt C Anderson 180 Maple 13
number-theory divisors
June 09 2021
0  2

Hi again all

You can enjoy this simple Maple code to find the proper divisors of a given positive integer (whole number).

 

Download proper_divisor.mw

here  some_proper_divisor_examples_3.pdf   file

Hope this helps

Matt Anderson

 

Why is this example from Maple not working?...

Asked by: mapleatha 160
animation plotting plots
June 07 2021
1  14

I came across this example from Maple's Online Help:

F := proc(t)
plottools[line]([-2,0], [cos(t)-2, sin(t)], color=blue),
plottools[line]([cos(t)-2, sin(t)], [t, sin(t)], color=blue),
plot(sin(x), x=0..t, view=[-3..7, -5..5]);
end proc:

animate(F,[theta],theta=0..2*Pi,background=plot([cos(t)-2,sin(t),t=0..2*Pi]),
scaling=constrained,axes=none);


Why is it not working?

Thank you!

mapleatha

How far can Maple go, in terms of displaying Ferma...

Asked by: Matt C Anderson 180
June 02 2021
1  1

Hi all,

Wrote some Maple code, and want to share.

I put less than eight hours into this exploration.

I have been running a computer to find Fermat factors for over ten years now.  Want support.

This is an exploration in some mostly composite whole numbers.

see 2_to_the_exponent_p_3.pdf

also see 

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/2_to_the_exponent_p_3.mw .
 

Download 2_to_the_exponent_p_3.mw

Have a look.

Regards

 

Matt

 

 

form of primes greater than four...

Asked by: Matt C Anderson 180
June 01 2021
1  4

Hi all,

Do some of you know the form of all primes greater than five?

Is it 6*p +/- 1?

maple_about_primes.pdf



Maple Worksheet - Error

et/convert/                                  Maple_prime_expository_number_one_for_sharing.mw .
 

Download Maple_prime_expository_number_one_for_sharing.mwmaple_about_primes.pdf

 

Regards,

Matt

 

plots showing sin(30°) without italics?...

Asked by: mapleatha 160
typesetting degree plotting
May 31 2021
2  8

The command

plots[textplot]({[4.2,1.9,typeset(sin(30°))]},
symbol=point,color=black,font=["Arial","Bold",25]);

will output

Can we have the number 30 in a non-italic format?
Thank you!
mapleatha


 

Integral of partial derivatives by combining the t...

Asked by: guras 100
derivative int multivariate-calculus
April 24 2021
0  2

Hi everyone!

Is it possible to evaluate int(diff(u(x,y,t),x)*v(x,y,t)+diff(v(x,y,t),x)*u(x,y,t),x) directly as u(x,y,t)*v(x,y,t) in Maple?

Thanks.

Do not give such terms...

Asked by: Maria2212 25
programming assumptions assume
April 01 2021
2  6

Good afternoon! I can not give such terms in the procedure, simplify and evalf do not help. Why does it write three lambdas instead of a lambda in a cube?
What can I do about it?

help.mw
 

 

Error, (in content/content) invalid arguments...

Asked by: panke 145
radical simplify simplification
March 25 2021
2  5

temp.mw

1 2 3 4 5 6 7 Last Page 2 of 54