## How to find abundant/deficient numbers...

Hi MaplePrimes,

I have decided to look for an odd perfect number.
I want Maple code to determine if there is a perfect number beween bounds on a factor.
For definitions of deficient and perfect and abundant numbers, see
https://mathworld.wolfram.com/DeficientNumber.html
https://mathworld.wolfram.com/PerfectNumber.html
https://mathworld.wolfram.com/AbundantNumber.html

I need some sort of loop that closes in on the edge of abundant/deficient numbers.

Specifically, an IF() statement is needed about wheather the function

sigma(a)-2*a is positive or negative.

Regards,

Matt

abundant_edge_30.mw

## Sum of proper divisors procedure...

Hi all,

I wrote a little simple minded procedure for finding the sum of divisors for a positive integer.  I thought I would share.

sum_of_proper_divisors.mw

sum_of_proper_divisors.pdf

Has someone written better code for this task?

Regards,
Matt

## idk what this line should do...

H:= Fit(a*B+b*sinh(c*B), Bdata, Hdata, B)

can someone help me ? i have never used maple13 and i have to describe a code but i cant understand what this line does. can someone help me, please?!
thank u in advance.

## Maple unter Window 11...

My Maple Program run perfectly  with Window 10.  My new Laptop has Window 11 and I get system  crash wenn I applied

plot3d(...) .

So I think It is somethink with the Grapics Redering ?

Can somebody help ?

Betriebsystemname    Microsoft Windows 11 Home
Version    10.0.22621 Build 22621
Weitere Betriebsystembeschreibung     Nicht verfügbar
Betriebsystemhersteller    Microsoft Corporation
Systemname    ELIPAN
Systemhersteller    SAMSUNG ELECTRONICS CO., LTD.
Systemmodell    750XED
Systemtyp    x64-basierter PC
Prozessor    12th Gen Intel(R) Core(TM) i5-1235U, 1300 MHz, 10 Kern(e), 12 logische(r) Prozessor(en)
BIOS-Version/-Datum    American Megatrends International, LLC. P08CFG.033.220913.HQ, 13.09.2022
SMBIOS-Version    3.4
Version des eingebetteten Controllers    255.255
BIOS-Modus    UEFI
BaseBoard-Hersteller    SAMSUNG ELECTRONICS CO., LTD.
BaseBoard-Produkt    NP750XED-KC5DE
BaseBoard-Version    SAMSUNG_SW_REVISION_12345+10.0.22000
Plattformrolle    Mobil
Sicherer Startzustand    Ein
PCR7-Konfiguration    Erweiterung zum Anzeigen erforderlich
Windows-Verzeichnis    C:\WINDOWS
Systemverzeichnis    C:\WINDOWS\system32
Startgerät    \Device\HarddiskVolume1
Gebietsschema    Deutschland
Hardwareabstraktionsebene    Version = "10.0.22621.819"
Benutzername    ELIPAN\Pan
Zeitzone    Mitteleuropäische Zeit
Installierter physischer Speicher (RAM)    16,0 GB
Gesamter physischer Speicher    15,7 GB
Verfügbarer physischer Speicher    8,28 GB
Gesamter virtueller Speicher    16,7 GB
Verfügbarer virtueller Speicher    6,98 GB
Größe der Auslagerungsdatei    1,00 GB
Auslagerungsdatei    C:\pagefile.sys
Kernel-DMA-Schutz    Ein
Virtualisierungsbasierte Sicherheit    Wird ausgeführt...
Virtualisierungsbasierte Sicherheit – erforderliche Sicherheitseigenschaften
Virtualisierungsbasierte Sicherheit – verfügbare Sicherheitseigenschaften    Allgemeine Virtualisierungsunterstützung, Sicherer Start, DMA-Schutz, UEFI-Code Readonly, SMM Security Mitigations 1.0, Modusbasierte Ausführungssteuerung, APIC-Virtualisierung
Virtualisierungsbasierte Sicherheit – konfigurierte Dienste    Durch Hypervisor erzwungene Codeintegrität
Virtualisierungsbasierte Sicherheit – ausgeführte Dienste    Durch Hypervisor erzwungene Codeintegrität
Windows Defender-Anwendungssteuerungsrichtlinie    Erzwungen
Windows Defender-Anwendungssteuerungs-Richtlinie für den Benutzermodus    Aus
Unterstützung der Geräteverschlüsselung    Erweiterung zum Anzeigen erforderlich
Es wurde ein Hypervisor erkannt. Features, die für Hyper-V erforderlich sind, werden nicht angezeigt.

## create a lib . Problem with Windows 11...

I get this error unter WINDOW 11, with WINDOW10 it was o.k

> lib := cat(kernelopts(mapledir), "/lib/mylibn.mla"); savelib('MYEIGW', lib)

Error, cannot lock repository, C:\Program Files\Maple 13/lib/mylibn.mla

 >
 >
 (1)
 >

## How do I solve a differential equation in Maple?...

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,

## How can one put 2 print() statements on the same l...

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

Regards,
Matt A

## Matt's procedure for pairs of prime numbers...

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

## second question about a simple two column matrix...

Thank you very much, Maple Community !

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

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...

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...

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...

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

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

 > рестарт;
 >
 >
 > 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;
 (1)
 > G:=2*E; êîëè÷åñòâî ãðàíè÷íûõ óñëîâèé
 (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]):  Äó âèäà 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)
 (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)
 (4)
 >

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

 > GP:=i0-1+i1-1;
 (5)
 >

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

 > i:=1:
 > while i
 >

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

 > i:=number0+1:
 > while i
 >

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

 > i:=number0+number1+1:
 > while i
 >

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

 > i: = число0 + число1 + nops (v0) +1:
 > в то время как i
 >
 >
 >

новогодний

 > sistema: = {}:

новенькие

 > время: = {}:
 >
 > для я от 2 до G делаю
 > sistema: = sistema union {GRAN [i]}:
 > od:
 > для я от 2 до G делаю
 > peremenie: = peremenie union {C [i]}:
 > od:
 >
 > коэфф: = решить (система, изменение):
 > коэфф:
 >
 > система [1];
 (6)
 > ГРАН [1];
 (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...

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);
 (5)
 > series(Ek4,k=0,4);
 > series(sqrt(Ek4),k=0,4);

## prime 2 tuple with difference of 36...

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

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