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

## more proper divisors from 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

## Proper Divisors Procedure in Maple code

Maple 13

Hi again all

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

Hope this helps

Matt Anderson

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

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

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.

also see

Maple Worksheet - Error

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

Have a look.

Regards

Matt

## form of primes greater than four...

Hi all,

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

Is it 6*p +/- 1?

Maple Worksheet - Error

et/convert/                                  Maple_prime_expository_number_one_for_sharing.mw .

Regards,

Matt

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