Download pak.mw

Is there a command in maple that shows which files were read during start up? I am confused as to which maple init file is being read on my PC.  From http://www.maplesoft.com/support/help/Maple/view.aspx?path=worksheet/reference/initialization

it says

"Under Windows, the initialization file is called maple.ini.
If <Maple>\lib\maple.ini exists, it is loaded first (where <Maple> is your Maple installation directory). With a  network installation of Maple, the commands in this initialization file will be executed by all users on the network.
To execute a user's personal set of commands, only the first initialization file in one of the following paths will be loaded.
1) The binary directory of your current working directory (for example, "c:\Program_Files\Maple\bin.win\maple.ini")
2) The <Maple>\Users directory (for example, "c:\Program_Files\Maple\Users\maple.ini")
3) The user's personal profile directory (multiuser only) (for example, "c:\Documents_and_Settings\userid\maple.ini")
      
Maple reads and executes the network initialization file before the personal initialization file."

How do I find out, from inside Maple, which file(s) were read? Or make maple shows a trace of the loading process to see what files ini files it is reading?

Maple 18.2 on windows. btw, the above help page seems old. I am on windows 7, and I do not have "c:\Program_Files\Maple\" folder. And do not have "c:\Documents_and_Settings" folder. So the above help is not very useful. May be it was written during windows 95 times?

Hi all,

I have some "boolean variable" constraint equation like this:

a1*x1+a2*x2+...+an*xn>=b1*y1+b2*y2+...+bn*yn

where a1,a2,...,an and b1, b2, ..., bn are 1 or -1

These equations will be used in LPSolve or the other command to find a group of parameters which can fit them.

Now I used for-loop to deal with this kind of question, for example:

But there are more than 10 boolean variables in my case and It's very inefficient. On the other hand, using for-loop to determine the equation we solve in the command will lead to great confusion.

I think there should be some ways able to solve this kind of "boolean variables" question in Maple, such as, through assume command to define the type of "boolean variable".

But I have no idea how to do it.

I have the following construct:
for i from 1 to 10 do
cubeprod:=i^3;
if irem(cubeprod,3)=0  and if(modp(cubeprod,2)<>0 then
cubesum(cubeprod);// I need both if statements to be true in order to invoke cubesum(). I've noticed that using an and between both if statements is incorrect but :

for i from 1 to 10 do
cubeprod:=i^3;
if irem(cubeprod,3)=0 then 
 if modp(cubeprod,2)<>0 then cubesum(cubeprod) fi;
end if;
end do;

gives me an error as well. What is the right syntax to achieve this?

I have the following procedure to do the above. It works but it returns [9,10],[10,9],[12,1] for n=1729(for example). How do I modify this to 

a) to count 9,10 and 10,9 as the same and hence only show one of them

b) get 1,12 to show as a solution?

cubesum:=proc(n::nonnegint)
global listcub:=table();
local k:=0, x:=iroot(iquo(n,3),3),y:=x,x3:=x^3,y3:=y^3;
if 3*x3 <> n then x=x+1; x3:=x^3;y:=x;y3:=x3 end if;
while x3<=n do
y:=iroot(n-x3,3); y3:=y^3;
if(x3+y3 = n) then k:=k+1; listcub[k]:=[x,y]end if;
x:=x+1; x3:=x^3;
end do;
convert(listcub,list);
end proc:

Hi
I have my question makes any sense. I am from Denmark and not used to write math in english.

I have an characteristic matrix with an variable λ that takes on differen values.

How do I write λ in the matrix so Maple knows that when I call out a row with the variable λ in it and asssign

λ to a specific value, Maple changes λ to the specific value.

Example (I was thinking something like this):

A:=Matrix(2,2,[(1-λ_i),2,3,(4-λ_i)])

λ_1:=2

λ_2:=4

A[1],λ_1               (1-2) 2

A[1],λ_2               (1-4) 2

A[2],λ_1               3 (4-2)

A[2],λ_2               3 (4-4)

Hi, I am an student and I am currently working on a system that sketches the relation of predator and prey of yellowstone's gray wolf and elk. I tried using the Lotka-Volterra model, but I wanted to add more parameters and add a carrying capacity for the system. Unfortunatley I cannot find a way to edit the Lotka model to my needs, and because I am new i do not know how to create my own model. This is the two equations I want to use: (D(x))(t) = alpha*x(t)-ax^2/k-b*x(t)*y(t)-gx(t), (D(y))(t) = -beta*y(t)+c*x(t)*y(t)-gy(t)

were k is carrying capacity.

Basically what I am asking is that if someone can help make the system workable on Maple and some steps of how to do it. 

What I tried is at

http://apfp.sourceforge.net/maple_interval_arith.pdf

I would like to be able to do comparisons of intervals.

At least take min or max - I thought I saw that on help pages but didn't get it to work.

Would also like to use with complex numbers - tried evalrC - but did not do what I expected.

Is there a way to tell Maple to export a figure to PDF with a proper bounding box, like all of its other graphics export formats?

To export a graphics produced by Maple's plot(), I right-click on the plot, select Export, and then one of the several choices of graphics formats.  All options, other than the PDF, work fine—they produce graphics files whose bounding boxes correspond to the extents of the image.  Exporting to PDF misbehaves—it produces the equivalent of a 8.5''x11'' paper and inserts the graphics somewhere near the upper left corner.  I am absolutely at a loss to see the utility of that.  What in the world is the use of such an export?

Is there a configuration setting that tells Maple to save to PDF with a proper bounding box?  I looked around but couldn't find one.

Below is my code:

> restart:
> q:=1.6*10^(-19):e0:=8.85*10^(-12):k:=1.38*10^(-23):t:=300:ni:=1.45*10^(16):eox:=3.9*e0:esi:=11.7*e0:
> tox:=2*10^(-9):tsi:=25*10^(-9):L:=1*10^(-6):cox:=(eox/tox):csi:=esi/tsi:a:=sqrt(q^2*ni/(2*k*t*esi)):V:=0:vgs:=2:u:=0.04:
> mm:=(cos((q*cox*(vgs-psi11)/(esi*k*t*a))*(exp(q*(V-psi11)/(2*k*t))))-((2*q*cox*(vgs-psi11)/(esi*k*t*a^2*(tsi)))*(exp(-q*(psi11-V)/(k*t))))):
> psi1:=solve(mm=0,psi11);
psi1 := 0.5434311697 - 0.003168514913 I
> mm:=(cos((q*cox*(vgs-psi22)/(esi*k*t*a))*(exp(q*(vds-psi22)/(2*k*t))))-((2*q*cox*(vgs-psi22)/(esi*k*t*a^2*(tsi)))*(exp(-q*(psi22-vds)/(k*t))))):
> qq2:=seq(solve(mm=0,psi22),vds=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2]):
> psi2:=(vector([qq2]));vds:=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2]:
psi2 := [0.5434311697 - 0.003168514913 I, 0.6398902333 - 0.003560933110 I,

0.7360983262 - 0.003949922699 I, 0.8586260376, 0.9544714953, 1.066786341,

1.145104060, 1.239782678, 1.333965608, 1.427671475, 1.521739902,

1.612717888 + 0.004670871971 I, 1.699606898 + 0.007295737134 I,

1.782802467 - 0.009645080463 I, 1.913366188, 1.969622161, 1.996631837,

1.999909512, 1.999998089, 1.999999960, 1.999999999]
> for i from 1 to 21 do
> Qgg[i] := Re( tsi*cox*L*(1/6*psi2[i]^3-1/2*psi2[i]^2-1/2*psi2[i]^2*vgs-1/2*psi2[i]^2*k*t/q+1/2*psi2[i]*vgs^2+psi2[i]*k*t*vgs/q+psi2[i]*vgs-8*psi2[i]*esi*k*t/(q*tsi*cox)-64*esi^2*k^2*t^2*ln(-8*esi*k*t-q*tsi*cox*vgs+q*tsi*cox*psi2[i])/(q^2*tsi^2*cox^2)-1/6*psi1^3+1/2*psi1^2+1/2*psi1^2*vgs+1/2*psi1^2*k*t/q-1/2*psi1*vgs^2-psi1*k*t*vgs/q-psi1*vgs+8*psi1*esi*k*t/(q*tsi*cox)+64*esi^2*k^2*t^2*ln(-8*esi*k*t-q*tsi*cox*vgs+q*tsi*cox*psi1)/(q^2*tsi^2*cox^2))/(k*t*(vg-psi2[i])/q+1/8*(vgs-psi2[i])^2+4*k^2*t^2*esi*ln(cox(vgs-psi2[i])+8*k*t*esi/(q*tsi))/(q^2*tsi*cox)-k*t*(vg-psi1)/q-1/8*(vgs-psi1)^2-4*k^2*t^2*esi*ln(cox(vgs-psi1)+8*k*t*esi/(q*tsi))/(q^2*tsi*cox)));
> end do:
> Qg:=simplify((Qgg));
/[ -15
Qg := TABLE\[1 = Float(infinity), 2 = -2.726896811 10 ,

-15 -15 -15
3 = -2.681258392 10 , 5 = -2.579829239 10 , 4 = -2.623431373 10 ,

-15 -15 -15
7 = -2.495397349 10 , 6 = -2.529657206 10 , 10 = -2.378078479 10 ,

-15 -15 -15
11 = -2.341729539 10 , 8 = -2.454901580 10 , 9 = -2.415751244 10 ,

-15 -15 -15
15 = -2.212746094 10 , 14 = -2.250797897 10 , 13 = -2.277803789 10 ,

-15 -15 -15
12 = -2.308132152 10 , 21 = -2.192572935 10 , 20 = -2.192572943 10 ,

-15 -15 -15
16 = -2.198907782 10 , 17 = -2.193213670 10 , 18 = -2.192589892 10 ,

-15]\
19 = -2.192573293 10 ]/

i have to plot Qg with respect to vds. the value of Qg are stored in a table which have mixed entries. also please suggest how can i set the x-axis of the plot in decimal notation?

Regards

Hello,

please explain how to write a code to calculate and output the actual area using integration for y=X^3 over range (0,2) using left-hand rule and 200 subdivisions?

Thank you 

I currently have a function quadsum(n) that determines the [x,y] solutions of the above equation for an integer n. :

quadsum:= proc(n::nonnegint)
local
k:= 0, mylist:= table(),
x:= isqrt(iquo(n,2)), y:= x, x2:= x^2, y2:= y^2;
if 2*x2 <> n then x:= x+1; x2:= x2+2*x-1; y:= x; y2:= x2; end if;
while x2 <= n do
y:= isqrt(n-x2); y2:= y^2;
if x2+y2 = n then k:= k+1; mylist[k]:= [x,y] end if;
x:= x+1; x2:= x2+2*x-1;
end do;
convert(mylist, list)
end proc:

How would I alter this so that I get [x,y] for n= (5^a).(13^b).(17^c)(29^d) for non-negative integers a,b,c,d?

Hi everyone,

Consider the following worksheet (Maple 13):

Euler.mw

 (or as a Google Drive link)

https://drive.google.com/file/d/0Bzr3EyK8arkOUkZPcEE5b0NySVU/view?usp=sharing

Choose execute entire worksheet on both Maple 13 and Maple 18.

Maple 13 is blindingly fast at 40 seconds, entire wroksheet upon completion

Maple 18 crawls and hangs forever.

Can anyone recognise what the problem with Maple 18 is?

Running on a two core 1.7GHz AMD Athlon.

Many thanks

--

Yiannis

Hi

I can not my software,I get this error"the specified procedure could not be found"

what should I do?

Given a  vector v (parallel to the principle axis) and a center point c,

and generating vectors u(beta) orthogonal to v for 0 < beta < 2pi,

and real positive parameters a and b 

I want to plot the hyperboloid  c + a*v*sec(alpha) + b*u(beta)*tan(alpha)

for -pi/2 < alpha < pi/2.

I am able to generate something of a plot, but can't control the size and appearance.