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

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

How do I get Maple 13, 16 and 17 to start on 64 bi...

Asked by: maple_linux 20
March 17 2021
4  0

 

Hi,

I am using Maple 13, 16 and 17 on 64 bit ubuntu 14.04LTS without 

any problems.

 

I tried to start the above Maple versions after upgrading to ubuntu 18.04LTS. Unfortunately, Maple does not start and transmits an incorrect hostid error message. The error message for Maple 16 is pasted below.

 

The command ifconfig -a produces an ethernet address 00:21:70:a6:a4:4b that is identical to the hostid given in the license.dat file.

 

I believe that the license manager flexnet is 32 bit. I used the synaptic 

manager to install lsb-base and ia32-libs (synaptic package manager refused to install lsb-core and transmitted some error messages). However, this did not help and I obtained the same hostid error message.

 

I do not know what is causing the hostid error message. Maybe some script file/s in the Maple installation can be modified such that Maple starts correctly?

 

Any assistance in getting Maple 13, 16, 17 to start on ubuntu 18.04LTS 

would be much appreciated.

 

Thank you very much for your support.
CF

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 

> ...:~/maple> maple16
> Maple initialization error, Invalid host.
> The hostid of this system does not match the hostid
> specified in the license file.
> Feature: Maple16
> Hostid: 002170a6a44b
> License path:
> /home/apps/maple16/toolbox/MapleSim/license/MapleSim.dat:/home/apps/maple16/license/license.dat:
> FLEXnet Licensing error:-9,57
> For further information, refer to the FLEXnet Licensing documentation,
> available at
> "https://linkprotect.cudasvc.com/url?a=https%3a%2f%2fwww.acresso.com&c=E,1,5ci9cywR8uvrXaVJ4IkoXH0TFBBbJwSK_TulrfrUDs5AyVQ4TXsHOVm5KWKAIhp4OZb9C4OAnmgylY1JC6fpbdZK9SjQLVrW012Uzxq18iDYC56NG2eDyHS8RPA,&typo=1".
.......
>
> root@<mailto:root@>...:~# ifconfig -a
> enp0s25: flags=4099<UP,BROADCAST,MULTICAST> mtu 1500
> ether 00:21:70:a6:a4:4b txqueuelen 1000 (Ethernet)
> RX packets 0 bytes 0 (0.0 B)
> RX errors 0 dropped 0 overruns 0 frame 0
> TX packets 0 bytes 0 (0.0 B)
> TX errors 0 dropped 0 overruns 0 carrier 0 collisions 0
> device interrupt 22 memory 0xf6ae0000-f6b00000
>
> lo: flags=73<UP,LOOPBACK,RUNNING> mtu 65536
> inet 127.0.0.1 netmask 255.0.0.0
> inet6 ::1 prefixlen 128 scopeid 0x10<host>
> loop txqueuelen 1000 (Local Loopback)
> RX packets 454 bytes 35138 (35.1 KB)
> RX errors 0 dropped 0 overruns 0 frame 0
> TX packets 454 bytes 35138 (35.1 KB)
> TX errors 0 dropped 0 overruns 0 carrier 0 collisions 0
 

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 

Have you ever heard of Vector asterisks ...

Asked by: Matt C Anderson 165
plotting
March 15 2021
0  2

Look

manual_instructions.mw

 

manual_instructions.pdf

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