Advanced Search

Harry Garst

259 Reputation

6 Badges

19 years, 24 days
Contact Harry Garst

MaplePrimes Activity


These are questions asked by Harry Garst

indexed versus symbol...

Harry Garst 259 Maple
indexed
February 24 2022
2 3

I don't know how to fix this. What's wrong? Is there an alternative way to include underscores in a cat command?

``

> 

restart; with(LinearAlgebra)
> 

komb := proc (n::integer, m::integer) local L, a, j; L := [seq(cat(a, "__", j, ",", j), j = 1 .. n)]; add(mul(k), `in`(k, combinat:-choose(L, m))); return % end proc

proc (n::integer, m::integer) local L, a, j, k; L := [seq(cat(a, "__", j, ",", j), j = 1 .. n)]; add(mul(k), `in`(k, combinat:-choose(L, m))); return % end proc

 (1)
> 

komb(5, 5)

`a__1,1`*`a__2,2`*`a__3,3`*`a__4,4`*`a__5,5`

 (2)
> 

whattype(`a__1,1`)

symbol

 (3)
> 

A := Matrix(5, 5, shape = diagonal, symbol = a)

Matrix(%id = 36893489895174306860)

 (4)
> 

abs(A)

a[1, 1]*a[2, 2]*a[3, 3]*a[4, 4]*a[5, 5]

 (5)
> 

whattype(a[1, 1])

indexed

 (6)
> 

simplify(komb(5, 5)-abs(A))

`a__1,1`*`a__2,2`*`a__3,3`*`a__4,4`*`a__5,5`-a[1, 1]*a[2, 2]*a[3, 3]*a[4, 4]*a[5, 5]

 (7)

``

Download indexed.mw

input for Hessian in package VectorCalcu...

Harry Garst 259 Maple 2020
vectorcalculus duplicate-question hessian
July 09 2021
1 1

Using the VectorCalculus package in Maple 2020.2, the Jacobian does accept my target function, whereas the Hessian does not. Both the Jacobian and the Hessian need an algebraic expression for its input.

What can I do to make it work?

kind regards, Harry

 

vectorcalculus_question.mw

 

Linear Algebra package incompatable with...

Harry Garst 259 Maple
July 08 2021
0 5

I am using Maple 2020.2 and it seems that the Linear Algebra package is incompatable with the Vector Calculus package:

Why not a friendly error message?
I want to use the Jacobian and the Hessian. Are there alternatives? MTM package can give a Jacobian, but no Hessian (but MTM is also tricky to use).

Harry

 

RootOf interpretation...

Harry Garst 259 Maple
solve rootof
March 11 2021
1 1

 

I was hoping to find a nice symbolic solution for the root of this equation:

a,c,d,k are parameters.
Unfortunately, I got an RootOf expression, which does not make sense to me. Assumptions do no help me further.

Next week, I have to present my results, but I don't know how to interpret is (the help file does not bring any relief). Is there something to say about the root? Suppose, I wanted to use side restriction on a problem? E.g., find a solution conditional on the fact that the derivative is zero. 


Clearly, in the plot it seems simple.

In one case there is a solution Maple provides: a  (a parameter)

This seems strange, because numerically there is one solution (and clearly not a).

> 

restart
> 

version(); interface(version)

 User Interface: 1502365
         Kernel: 1502365
        Library: 1502365

  

1502365

  

`Standard Worksheet Interface, Maple 2020.2, Windows 10, November 11 2020 Build ID 1502365`

 (1)
> 

f := proc (t) options operator, arrow; c/(1+exp(k*(t-a)))-c*(t-a)*k*exp(k*(t-a))/(1+exp(k*(t-a)))^2+d/(1+exp(-k*(t-a)))+d*(t-a)*k*exp(-k*(t-a))/(1+exp(-k*(t-a)))^2 end proc

proc (t) options operator, arrow; c/(1+exp(k*(t-a)))-c*(t-a)*k*exp(k*(t-a))/(1+exp(k*(t-a)))^2+d/(1+exp(-k*(t-a)))+d*(t-a)*k*exp(-k*(t-a))/(1+exp(-k*(t-a)))^2 end proc

 (2)
> 

assume(k > 0)
> 

additionally(a > 0)
> 

additionally(t > 0, t < 21)
> 

about(a, k, t)

Originally a, renamed a~:
  is assumed to be: RealRange(Open(0),infinity)

Originally k, renamed k~:
  is assumed to be: RealRange(Open(0),infinity)

Originally t, renamed t~:
  is assumed to be: RealRange(Open(0),Open(21))
 

  
> 

interface(showassumed = 0)

0

 (3)
> 

solve(f(t) = 0, t)

Warning, solve may be ignoring assumptions on the input variables.

  

(k*a+RootOf(-(exp(_Z))^4*d+(exp(_Z))^3*_Z*c-(exp(_Z))^3*_Z*d-(exp(_Z))^3*c-3*(exp(_Z))^3*d+2*c*_Z*(exp(_Z))^2-2*(exp(_Z))^2*_Z*d-3*c*(exp(_Z))^2-3*d*(exp(_Z))^2+c*_Z*exp(_Z)-d*_Z*exp(_Z)-3*c*exp(_Z)-d*exp(_Z)-c))/k

 (4)
> 

[allvalues(%)]

[(k*a+I*Pi+(2*I)*Pi*_Z1)/k, (k*a+RootOf(-d*(exp(_Z))^2+c*_Z*exp(_Z)-d*_Z*exp(_Z)-c*exp(_Z)-d*exp(_Z)-c))/k]

 (5)
> 

Re([op(%)][1])

a

 (6)
> 

Re([op(`%%`)][2])

(k*a+Re(RootOf(-d*(exp(_Z))^2+c*_Z*exp(_Z)-d*_Z*exp(_Z)-c*exp(_Z)-d*exp(_Z)-c)))/k

 (7)
> 

a := 'a'; b := 'b'; c := 'c'; d := 'd'; k := 'k'; t := 't'

a

  

b

  

c

  

d

  

t

 (8)

Now, numerically:

> 

b := 40; k := .3; a := 10; c := 2.8; d := -1.5

40

  

.3

  

10

  

2.8

  

-1.5

 (9)
> 

RealDomain:-solve(f(t) = 0, t)

11.02361521

 (10)
> 

plot(f(t), t = 0 .. 20, color = "DarkBlue", thickness = 2)

 
> 

``


 

Download Please_NO_RootOf.mw

Maple 2020 frozen...

Harry Garst 259 Maple 2020
December 05 2020
0 0

The last months Maple 2020 suffers from frozen interface behavior:

I simply wanted to save a file, but that resulted in a frozen window:

After 10 minutes the file explorer window appeared, but also this screen is unresponsive.

I reinstalled Maple last month, but that didn't help. Remarkably, the problem comes and goes. It seems to me that some java problem is making my life difficult.
Any ideas?
kind regards,

Harry

 

ps I ran the process monitor of Sysinternals, but I have no idea if this has any relevance to the problems:

11460 RegOpenKey HKCR\Local Settings\Software\Microsoft\Windows\CurrentVersion\AppModel\PackageRepository\Extensions\windows.protocol\Maple.cwmaple.2020 NAME NOT FOUND Desired Access: Read

11460 ReqQueryValue HKCR\Maple.cwmaple.2020\URL Protocol NAME NOT FOUND Length: 12

11460 RegOpenKey HKCU\Software\Classes\Maple.cwmaple.2020  NAME NOT FOUND Desired Access: Maximum Allowed

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