Maple 13 Questions and Posts

These are Posts and Questions associated with the product, Maple 13

 

restart

A := Omega; assume(Omega > 0); B = Omega

B = Omega

(1)

``

Download assumetest.mw

By printing each processing step for a parametric surface entered in 3-dimensional space,

i- Finding whether the surface is a ruled(regle) surface and printing it in the ruled surface form,

ii- Calculating the dispersion parameter of the surface and the throat line,

iii- Draws the surface and the throat line on the surface together,

Simple question

Got stuck

tup_10_g_trying.mw

Want output to look like - 

tuple10(50,12,14,0,0,0,0,0,0,0,0)

1 5 17 19
2 17 29 31
3 29 41 43
4 47 59 61

Please assist.
Regards,
Matt A

Hi again everybody,
I hope this can help in some classrooms.

This procedure has some warnings.  Someone else can clean it up.
Also, the ithprime() command can make it a bit slicker.
I did work with prime constelations project, and you can see my web page at mattanderson.fun  .

Refference oeis.org/A40 .
Some prime constelations, and k-tuples, are not of 'general intrest' yet, so they are not in the OEIS encyclopedia YET!

But prime numbers are exciting.  I like doing calculations.

Hopefully someone in the next generation will take up an intrest and do what I did.

pairs_procedure.mw

pairs_procedure.pdf

I wish the best on all of you.

Regards,

Matt C. Anderson

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

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

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

При использовании «объединения» нумерация не сохраняется. Непонятно, по каким алгоритмам «объединение» применяется к тому или иному месту. Почему 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

 

 

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

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

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

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

 

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

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

 

 

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

 

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