Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

how to find infinity norm for function of fraction exponenctial

I am wanting to plot a phase portrait of du/dt against u.

 

I have attempted to draw this on maple however am having great difficulty. Can anyone recommend what I may be doing wrong or how I should go about doing this?

Load the Tetrads package, and choose Schwarzshild spacetime:

with(Physics):
with(Tetrads):
Setup(metric = schwarzschild,quiet):

Consider then the following code/output:

expr := D_[~mu](e_[mu,a]);
seq(simplify(SumOverRepeatedIndices(expr             )),a = 1..4);
seq(simplify(SumOverRepeatedIndices(D_[~mu](e_[mu,a]))),a = 1..4);

In my mind, the last two code lines are entirely identical, but their corresponding outputs are certainly not. I don't get it. The above can be compared with the following code/output in which the vierbein indices mu and a have been interchanged:

expr := D_[~mu](e_[a,mu]);
seq(simplify(SumOverRepeatedIndices(expr             )),a = 1..4);
seq(simplify(SumOverRepeatedIndices(D_[~mu](e_[a,mu]))),a = 1..4);

Here, the outputs are identical. The mystery becomes total when realizing that according to Maple, e_[a,mu] and e_[mu,a] are considered completely identical [a notational choice I do not find reassuring].

 

I've been doing some soul searching and started to think that the maple interface and i are kindred spirts, conjoined spirit twins if you will, and i can confirm this if this is the manner in which the totient function is calculated:

 

For example how to input the mixed fraction 3 1/7? What are keys for super- and subscript?

I got decimal places problem,

not known correct or not

the value are different but difference is constant 

i do not know how many places needed to get exact result

i do not believe the difference is constant

because the matrix are different

but even if using 36 decimal places still constant, 

i notice increasing decimal places , the constant difference is changed

is it possible  to output fraction when calculate eigenvector?

Partial_for_R0.mw
Anybody knows what happen to this coding?

restart

unprotect(Pi)

k := .5;

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(1)

R := k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))

65688.31239

(2)

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), k)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ε`)

Error, invalid input: diff received 0.2e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), tau)

Error, invalid input: diff received .95, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), phi)

Error, invalid input: diff received 0.25e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), chi)

Error, invalid input: diff received 0.598e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), p)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ϑ`)

Error, invalid input: diff received 0.8e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), mu)

Error, invalid input: diff received 0.1e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), alpha)

Error, invalid input: diff received 0.57e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), rho)

Error, invalid input: diff received 0.5e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), beta)

Error, invalid input: diff received 0.115e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), eta)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), q)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

NULL

NULL


 

Download Partial_for_R0.mw
 

restart

unprotect(Pi)

k := .5;

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(1)

R := k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))

65688.31239

(2)

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), k)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ε`)

Error, invalid input: diff received 0.2e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), tau)

Error, invalid input: diff received .95, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), phi)

Error, invalid input: diff received 0.25e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), chi)

Error, invalid input: diff received 0.598e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), p)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ϑ`)

Error, invalid input: diff received 0.8e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), mu)

Error, invalid input: diff received 0.1e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), alpha)

Error, invalid input: diff received 0.57e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), rho)

Error, invalid input: diff received 0.5e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), beta)

Error, invalid input: diff received 0.115e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), eta)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), q)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

NULL

NULL


 

Download Partial_for_R0.mw

 

 

Basic_Reproduction_Number.mw
Hello, anybody knows how to compute sensitivity analysis for basic reproduction number using Maple coding? By taking the partial derivatives of all the parameters involved in the basic reproduction number using normalised forward sensitivity index formula. Thanks.

restart

unprotect(Pi)

R[0] := k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

(1)

``


 

Download Basic_Reproduction_Number.mw

 

I am having difficulties with a recursive function call with named parameters.

The proc is defind as follows:

Subs:=proc(eqn::seq(equation),elemt::Element,{num::boolean:=false},$) option remember;
...
qs[i]=Subs(eqn,elemt[qs[i]]); # This is where I need to add something like 'num'=num
...
end proc;

My difficulty is with the "num" flag. In the code this flag governs whether to try to evalf() certain results or not. That and the recursion per se all works. What does not work is when I add the num option to the internal call. I have been trying 'num'=num, ''num''=num, "num"=num and `num`=num.  I either get an error that something like true=true is not a valid argument here, or that "num" is not valid.

Any hint?

Thanks,

Mac Dude

Let G be a permutation group like S3. I have a loop in which x runs over the elements of this group.

When I ask the cycle type of x or the number of orbits of x, Maple gives an incorrect result in some cases. This is because when x=(12), for example, Maple considers x to be an element of the permutation group S2 instead of S3, i.e. it fails to understand that (12) should really be (12)(3) in this case.

Is there a way I can instruct Maple to consider the permutation (12) as a member of S3 instead of S2? 

I am posting the code as requested. Something like this:

with(GroupTheory);

for x in Elements(S3) do print(x, numelems(Orbits(PermutationGroup(x)))) end do;

(I don't know how to format this, sorry)

This code is simply asking how many orbits exist in the action of each element of S3. The result will be 1 for x=(12), however, because Maple will assume that (12) is acting on the set {1,2}, even though I am considering (12) as an element of S3. I want it to consider (12)=(12)(3) and let it act on the set {1,2,3}, thereby producing 2 orbits, not 1.

Ok embarassing I posted this one and now it outputs this java error file log and yep i dont know that much about java at all. 

 

Sometimes it works, sometimes it closes maple automatically and spits out the .txt error log i have attached if thats any help... MAPLE_EXAMPLE_16.mw  i honestly have never ever seen maple behave this way in more than a decade of playing around in it with stuff like this.

EDIT: Here is a version with smaller plot components that seems to be working.... MAPLE_EXAMPLE_17.mw

It says i am not allowed to upload a log file but this is what it looks like:

#
# A fatal error has been detected by the Java Runtime Environment:
#
# java.lang.OutOfMemoryError: requested 1024000 bytes for GrET in C:\BUILD_AREA\jdk6_18\hotspot\src\share\vm\utilities\growableArray.cpp. Out of swap space?
#
 
V  [jvm.dll+0x15df8a]
V  [jvm.dll+0x1e1e14]
V  [jvm.dll+0x1a1aad]
V  [jvm.dll+0xc834f]
V  [jvm.dll+0xca01c]
V  [jvm.dll+0xca370]
V  [jvm.dll+0xce42a]
V  [jvm.dll+0x1d8592]
V  [jvm.dll+0xc9398]
V  [jvm.dll
Java Threads: ( => current thread )
  0x522e7c00 JavaThread "Timer-28" [_thread_blocked, id=23236, stack(0x5f2f0000,0x5f6f0000)]
  0x522e7400 JavaThread "Timer-22" [_thread_blocked, id=24248, stack(0x6b440000,0x6b840000)]
  0x522e2c00 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=11140, stack(0x6a2d0000,0x6a6d0000)]
  0x522e7000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=6792, stack(0x64d30000,0x65130000)]
  0x522e6400 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=20280, stack(0x6a6d0000,0x6aad0000)]
  0x522e4400 JavaThread "Timer-18" daemon [_thread_blocked, id=23888, stack(0x63f20000,0x64320000)]
  0x522e6800 JavaThread "Timer-17" [_thread_blocked, id=17320, stack(0x69ed0000,0x6a2d0000)]
  0x522e5000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=13480, stack(0x69ad0000,0x69ed0000)]
  0x522e4800 JavaThread "Timer-16" daemon [_thread_blocked, id=18184, stack(0x68ed0000,0x692d0000)]
  0x522e3000 JavaThread "Timer-15" [_thread_blocked, id=23296, stack(0x692d0000,0x696d0000)]
  0x522e3800 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=16388, stack(0x66820000,0x66c20000)]
  0x522e1400 JavaThread "Timer-14" daemon [_thread_blocked, id=24500, stack(0x61c50000,0x62050000)]
  0x522e0c00 JavaThread "Timer-13" [_thread_blocked, id=15848, stack(0x659b0000,0x65db0000)]
  0x522e0000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=7260, stack(0x655b0000,0x659b0000)]
  0x466cb000 JavaThread "Timer-12" daemon [_thread_blocked, id=1112, stack(0x64320000,0x64720000)]
  0x4699b800 JavaThread "Timer-11" [_thread_blocked, id=24540, stack(0x64930000,0x64d30000)]
  0x4699b400 JavaThread "Timer-10" daemon [_thread_blocked, id=21144, stack(0x60bf0000,0x60ff0000)]
  0x4699d400 JavaTh
 

Hi. I have an excel data file.

IDEALLY, I want to import the results of the formula in cells C3:C7.

But it doesn't like it. Reports 0's

When I convert the range to result to values in cells D3:D7 it imports no problem. But I don't want to do this step.

any ideas?

Book1.xls

import.mw

 

 

 

Lett_L_rot_AL_MapleQuestion.mws

If the attached program is run, at the end are two animation outputs.  One is of two letters A; the other of the single second A being animated.   The second letter A animation works fine, but in the former it "blips" on Continuous animation, and on 'single cycle' the A totally disappears.  I'm using the same number of frames, 48, for each. 
   Prior to that is the letter L which is meant to give the illusion of it rotating, and the letters A getting out of the way.

   Any help on ridding the "blip" most appreciated. 

 

Hello, 

I study mainly subjects that fall under umbrella of number theory, but i have specified a little further in the worksheet. This is really a request for assistance, because in as much as i have met so many brilliant people online via social media etc,  I would always love to meet more, and especially ones who are more experienced in this field. 

 

Basically i am too cheap and old to think about going to a good university, so I am trying to get free advice from the people who have probably completed doctorates in the relevant field. Got to be honest I say.

 

Anyway my contact email is at the top of the attached worksheet.

 

First thing that stood out to me about the distributions produced in this worksheet is how sparse the number of points is for N=17 relative to all the other values of N.

EXAMPLE_FOR_MAPLE3.mw

 

Edit: Another example worksheet added.

MAPLE_EXAMPLE_13.mw

Hello, 
I have this code to solve partial equation, and I need to writdata in file .dat or.txt using table and loops or if there is another way 
Any help please,
Thank you 

 

``

restart

``

eq := diff(Tf(x, y), x) = diff(Tf(x, y), y, y):

Cl := {1 = Tf(0, y), (D[2](Tf))(x, 0) = 1, (D[2](Tf))(x, 1) = 1};

{1 = Tf(0, y), (D[2](Tf))(x, 0) = 1, (D[2](Tf))(x, 1) = 1}

(1)

``

pds := pdsolve(eq, Cl, numeric);

_m221385792

(2)

``


 

Download mapleprime.mw

First 841 842 843 844 845 846 847 Last Page 843 of 2216