Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi,

I am the administrator of Maple in my school, and all the students use Maple in part of their exams. Is it possible  to block the access to ChatGPT thru eg. the firewall or otherwise during exams. 

The reason for this question is that the students must have access to some internet sources during exams, but definately not CharGPT.

Kind regards 

Per Kirkegaard

Hello, I try to solve the equations of the odometric model with the Maple 2024 but I have not the answers as with the hands, can you help me to verify it ?

dsolve(diff(phi(t), t) = tan(10*t)/5)

dsolve(diff(x(t), t) = V*cos(ln(1 + tan(10*t)^2)/100))

dsolve(diff(y(t), t) = V*sin(ln(1 + tan(10*t)^2)/100))

Best regards, Edern Ollivier.

Hello Everyone;

I need to find the bifurcation point and further bifarcation diagram for the given model. But facing error. Can anybody help to do this? Can you refer some library for bifurcation analysis of ODE's? Code is attched. Thanks in Advance. 

123.mw

 

 

 

 

restart

C_m := 1.0; g_K := 36.0; I_inj := 0; g_L := .3; E_Na := 50.0; E_K := -77.0; E_L := -54.4

alpha_m := (.1*(V-25.0))/(1-exp(-(V-25.0)*(1/10))); beta_m := 4*exp(-V/(18.0)); alpha_h := 0.7e-1*exp(-V/(20.0)); beta_h := 1/(1+exp(-(V-30)*(1/10))); alpha_n := (0.1e-1*(V-10.0))/(1-exp(-(V-10.0)/(10.0))); beta_n := .125*exp(-V/(80.0)); I_Na := g_Na*m^3*h*(V-E_Na); I_K := g_K*n^4*(V-E_K); I_L := g_L*(V-E_L)

.125*exp(-0.1250000000e-1*V)

(1.1)

eq1 := (I_inj-I_Na-I_K-I_L)/C_m; m := alpha_m/(alpha_m+beta_m); n := alpha_n/(alpha_n+beta_n); h := alpha_h/(alpha_h+beta_h)

-16.32000000-1.000000000*g_Na*m^3*h*(V-50.0)-36.00000000*n^4*(V+77.0)-.3000000000*V

 

.1*(V-25.0)/((1-exp(-(1/10)*V+2.500000000))*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V)))

 

0.1e-1*(V-10.0)/((1-exp(-.1000000000*V+1.000000000))*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V)))

 

0.7e-1*exp(-0.5000000000e-1*V)/(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3)))

(1.2)

bif_eq1 := eq1 = 0;

-16.32000000-0.7000000000e-4*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)*(V-50.0)/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))-0.3600000000e-6*(V-10.0)^4*(V+77.0)/((1-exp(-.1000000000*V+1.000000000))^4*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V))^4)-.3000000000*V = 0

bif_eq2 := diff( eq1, V) = 0;

-0.2100000000e-3*g_Na*(V-25.0)^2*exp(-0.5000000000e-1*V)*(V-50.0)/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))+0.2100000000e-4*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)*(V-50.0)*exp(-(1/10)*V+2.500000000)/((1-exp(-(1/10)*V+2.500000000))^4*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))+0.2100000000e-3*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)*(V-50.0)*(.1/(1-exp(-(1/10)*V+2.500000000))-0.1000000000e-1*(V-25.0)*exp(-(1/10)*V+2.500000000)/(1-exp(-(1/10)*V+2.500000000))^2-.2222222222*exp(-0.5555555556e-1*V))/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^4*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))+0.3500000000e-5*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)*(V-50.0)/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))+0.7000000000e-4*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)*(V-50.0)*(-0.3500000000e-2*exp(-0.5000000000e-1*V)+(1/10)*exp(-(1/10)*V+3)/(1+exp(-(1/10)*V+3))^2)/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3)))^2)-0.7000000000e-4*g_Na*(V-25.0)^3*exp(-0.5000000000e-1*V)/((1-exp(-(1/10)*V+2.500000000))^3*(.1*(V-25.0)/(1-exp(-(1/10)*V+2.500000000))+4*exp(-0.5555555556e-1*V))^3*(0.7e-1*exp(-0.5000000000e-1*V)+1/(1+exp(-(1/10)*V+3))))-0.1440000000e-5*(V-10.0)^3*(V+77.0)/((1-exp(-.1000000000*V+1.000000000))^4*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V))^4)+0.1440000000e-6*(V-10.0)^4*(V+77.0)*exp(-.1000000000*V+1.000000000)/((1-exp(-.1000000000*V+1.000000000))^5*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V))^4)+0.1440000000e-5*(V-10.0)^4*(V+77.0)*(0.1e-1/(1-exp(-.1000000000*V+1.000000000))-0.1000000000e-2*(V-10.0)*exp(-.1000000000*V+1.000000000)/(1-exp(-.1000000000*V+1.000000000))^2-0.1562500000e-2*exp(-0.1250000000e-1*V))/((1-exp(-.1000000000*V+1.000000000))^4*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V))^5)-0.3600000000e-6*(V-10.0)^4/((1-exp(-.1000000000*V+1.000000000))^4*(0.1e-1*(V-10.0)/(1-exp(-.1000000000*V+1.000000000))+.125*exp(-0.1250000000e-1*V))^4)-.3000000000 = 0

 

 

 

bif_sol := solve({ bif_eq1,bif_eq2}, {V, g_Na});

Warning, solutions may have been lost

 

 

as the solutions, which are then expressed as the points mu, y via

   

[Back to ODE Powertool Table of Contents]

 

 

I have solution to system of equations that results (I assign it to) in:

I2sol := -I*omega0*k*L*vin*1/(L^2*k^2*omega0^2 + R^2)

I then try to solve it for k by doing

solve(I2sol = I2, k)

but that doesn't work.  What is the "right" Maple way to rearrange I2 such that the expression is the solution(s) for k?

Hey! I need help ASAP, because my maple file has been corrupted and i dont know what to do. Do you guys know how to recover a file? i can save it again as_mw. but should i change it to xml? or how? i have the link to my maple file attached, so if someone can help me, it could be helpful! Because i have an upcoming exam. Thanks youu

If a type is not known, an error is thrown

type([],foo)
Error, type `foo` does not exist

Since no error is thrown, these types are known

type({},'{}');# why that output?
                             false

type([],'[]')
                              true

I would have expceted {} and [] to be listed as subtypes of set and list since their counterparts (nonemptylist and nonemptyset) exist. Technically the types {} and [] are not needed since negating

not(type([],nonemptylist));
not(type({},nonemptyset));
                              true

                              true

works.  However, the types exist, hence my question

I keep getting an error: invalid subscript selector on specifically this part of my code:

MuligIndgange := subsop(RemoveList\[i\]=NULL, MuligIndgange).
Whenever I comment out the MuligIndgange := part it works.

I can't for the life of me see why it's happening, and whenever I use this code outside my proc it works aswell.

Any help is appreciated, and if you need more context I'll be happy to give it.

for i to numelems(MuligIndgange) do
     if M1[MuligIndgange[i][1], MuligIndgange[i][2]] =/= 0 then
         RemoveList := [op(RemoveList), i];
      end if;
end do;
for i from numelems(RemoveList) by -1 to 1 do
    MuligIndgange := subsop(RemoveList[i] = NULL, MuligIndgange);
end do;

I have seen this quite a bit in blocks of code. The `>` symbol seems to appear erratically. I don't know how to specifically reproduce this. Does it mean something? I would post the worksheet but it will not run without the package.

Hello

I have the following procedure that uses the Lie Derivatives of a vector field to build a set of equations.

LieDerList:=proc(h::symbol,f::list,vars::list)
description "This function returns the system of equations based on the Lie derivative.":
local i,n:=numelems(vars),L:=Array(0..n):
L[0]:=h:
for i from 1 to n do
    L[i]:=inner(f,map((a,b) -> diff(b,a),vars,L[i-1])):
end do:
return(zip((w,v)->simplify(w-v),[seq(L)],[seq](cat(h,i),i=0..n))):
end proc:

Below it is an example on how to call the procedure.

I used CodeTools:-ThreadSafetyCheck to check all the procedures used within LieDerList and LieDerList itself, but nothing wrong came out. However when I try to run 

LieEq4:=Threads:-Map(w->LieDerList(x,w,[x,y,z]),models4):

where models4 is a list of 1765 elements, maple returns "Error, (in factors) attempting to assign to `LinearAlgebra:-Modular:-Create` which is protected". If I change Threads to Grid, there is no problem at all.  

What am I overlooking? Is there a method to ensure the procedure is thread-safe?

Many thanks.   

PS.  I found one problem - inner, which is related to LinearAlgebra package, is not thread-safe.  

For a lesson I'm preparing, I want to illustrate some probability concepts using Maple.
In particular, I need to use the fact that the Expectation operator(let say the Mean operator) is a linear operator with respect to random variables. 
However, I don't want to particularize my demonstration by using this or that statistical distribution but simply the notion of random variable.
I have therefore created a specific Distribution named  MinimalAbstractDistribution in which only the mean and variance are defined.

When Statistics:-Mean is applied to the expression (A*p+q) where p and q are names and A is a random variable with distribution MinimalAbstractDistribution, linearity is effectively used.
But not when it is applied to A/p or A-q.

Why that?
Is there a way of defining a statistical distribution so that Mean behaves as expected?

(You will easily understand that no workaround of the form

Mean(A+q);
eval(%, q=-q);

# or
add(Mean~([op(A-q)]))

can be accepted in a lesson).

Thanks for your attention
 

restart

with(Statistics):

 

MinimalAbstractDistribution := proc(i)
  Distribution(Mean=mu[i], Variance=sigma[i]^2)
end proc:

A := RandomVariable(MinimalAbstractDistribution(1))

_R

(1)

[Mean, Variance](A);

[mu[1], sigma[1]^2]

(2)

Mean(p*A+q);
Mean(p*(A+q))

p*mu[1]+q

 

p*(mu[1]+q)

(3)

# But

Mean(A-q);
Mean(A/p)

FAIL

 

FAIL

(4)
 

 

Download Mean_operator.mw

Can I draw such figures using Maple?

If possible please guide how I can do it.

How to evaluate the right eigenvector of a given matrix in maple?

Simple question (I'm sure) here.  I have one worksheet where the variable does not show up in the ouput of an assignment like this:

whereas I was expecting the output to include the variable being assigned like this in most other worksheets:

I am having some difficulting locating the cause (and therefore the setting) to have the behavior set back to including the variable name in the ouput.  Thanks in advance for your help.

det_roots.docx

How to find the roots of the equation ?

like if the roots of the equations are m ,n,o,p then the value of these roots in terms of λ0  λ1  λ2  λ3  µ1  µ2 and  µ3

Calculate the values of the following determinants

First 24 25 26 27 28 29 30 Last Page 26 of 2229