MaplePrimes Questions

This is most likely a simple question for the power users of this forum, but I do not manage to find a solution. I have date in an Excel file. The first column consists of dates (mm/dd/yyyy) and the second of time (hh:mm). I can easily concatenate both in Maple using cat("9/7/2023", "10:22") but how can I convert the obtained string into a date+time that Maple understands? 

Thank you in advance for your help.

 

 

I create a system of equations (with 9 linear equations and 7 variables).
I get 7 equations from the multiplication of a matrix M with a transposed vector of the 7 variables equals the transposed vector of the 7 variables. Other two simple equations are necessary because they are restrictions. Those two equations a re very simple and thave the 7 variable sin it.
To start I have been trying with fsolve but i haven't been able to solve it yet, as I also get the error: "Error, (in fsolve) number of equations, 9, does not match number of variables, 7"
Have you andy idea to solve this problem? Thanks in advance.

Below is the plot like I want it. The basic plot has been done with a simple plot command.

>plot(0, x = 0 .. 10, y = 0 .. 4, gridlines = true)

But the label of each axis was done manually. But I have tried to do it inside the plot command. A little help would be very appreciated.

So here is what I want the plot to look like:

Thank you in advance for the help.

Mario

Hi everyone,

Do you know if there are some known codes for plotting pseudospectra of squared (and finite) matrices with given spectrum ? 
Thanks in advance

Best,

Rachid

What will be the range of p and q to get the plot and to get the optimum solution?
If possible get a solution for particular value of p and q.
file attached: q1.mw

Hi all, any one help  me to find the values of constants by using given condition and then how to varify that the goiven condition varify the expression. I have found manually but want to varify through maple.

help.mw

I have a list:
mylist := [x^4 + (-4*m - 7)*x^3 + (m + 4)*x^2 + (3*m - 5)*x, x^4 + (-4*m - 7)*x^3 + (m + 5)*x^2 + (5*m - 7)*x, x^4 + (-4*m - 7)*x^3 + (m + 5)*x^2 + (7*m - 5)*x, x^4 + (-4*m - 7)*x^3 + (2*m + 5)*x^2 + (3*m - 5)*x, x^4 + (-4*m - 7)*x^3 + (3*m + 1)*x^2 + (7*m - 10)*x, x^4 + (-4*m - 5)*x^3 + (2*m + 1)*x^2 + (7*m - 9)*x]

I use

L := map~(normal, mylist);

and get. 

L := [x^4 - (4*m + 7)*x^3 + (m + 4)*x^2 + (3*m - 5)*x, x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (5*m - 7)*x, x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (7*m - 5)*x, x^4 - (4*m + 7)*x^3 + (2*m + 5)*x^2 + (3*m - 5)*x, x^4 - (4*m + 7)*x^3 + (3*m + 1)*x^2 + (7*m - 10)*x, x^4 - (4*m + 5)*x^3 + (2*m + 1)*x^2 + (7*m - 9)*x].

I use seq to list L[i] and diff(L[i])

[seq([L[i], diff(L[i], x)], i = 1 .. nops(L))];

and get

[[x^4 - (4*m + 7)*x^3 + (m + 4)*x^2 + (3*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 4)*x + 3*m - 5], [x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (5*m - 7)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 5)*x + 5*m - 7], [x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (7*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 5)*x + 7*m - 5], [x^4 - (4*m + 7)*x^3 + (2*m + 5)*x^2 + (3*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(2*m + 5)*x + 3*m - 5], [x^4 - (4*m + 7)*x^3 + (3*m + 1)*x^2 + (7*m - 10)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(3*m + 1)*x + 7*m - 10], [x^4 - (4*m + 5)*x^3 + (2*m + 1)*x^2 + (7*m - 9)*x, 4*x^3 - 3*(4*m + 5)*x^2 + 2*(2*m + 1)*x + 7*m - 9]]

How can I insert L'(2), L''(2) and solve the systems L'(2) = 0 and L''(2) > 0 to get the solutions m?
like this
[seq([L[i], diff(L[i], x), solve([L'(2) = 0,L''(2)>0],m) ], i = 1 .. nops(L))]

I also tried
[seq([L[i], diff(L[i], x), eval(diff(L[i], x), x = 2), solve([eval(diff(L[i], x), x = 2) = 0], m)], i = 1 .. nops(L))]

to obtain 
[[x^4 - (4*m + 7)*x^3 + (m + 4)*x^2 + (3*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 4)*x + 3*m - 5, -41 - 41*m, {m = -1}], [x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (5*m - 7)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 5)*x + 5*m - 7, -39 - 39*m, {m = -1}], [x^4 - (4*m + 7)*x^3 + (m + 5)*x^2 + (7*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(m + 5)*x + 7*m - 5, -37 - 37*m, {m = -1}], [x^4 - (4*m + 7)*x^3 + (2*m + 5)*x^2 + (3*m - 5)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(2*m + 5)*x + 3*m - 5, -37 - 37*m, {m = -1}], [x^4 - (4*m + 7)*x^3 + (3*m + 1)*x^2 + (7*m - 10)*x, 4*x^3 - 3*(4*m + 7)*x^2 + 2*(3*m + 1)*x + 7*m - 10, -58 - 29*m, {m = -2}], [x^4 - (4*m + 5)*x^3 + (2*m + 1)*x^2 + (7*m - 9)*x, 4*x^3 - 3*(4*m + 5)*x^2 + 2*(2*m + 1)*x + 7*m - 9, -33 - 33*m, {m = -1}]]

I need to extract the list of all powers present in multivariable polynomial. For example, for x+y+x^2y^2+y^d I want to get the list [(1,0),(0,1),(2,2),(0,d)]. How to perform this?

Hello Everyone;

I need help to solve the following system using rk-4 method

restart;
NULL;
NULL;
C := 1.0;
gK := 36.0;
gNa := 120.0;
gL := 0.3;
VK := -77.0;
VNa := 50.0;
VL := -54.0;
III := 20;
alpha_n := 0.01*(v(t) + 55.0)/(1 - exp(-1.0/10.0*v(t) - 11.0/2.0));
beta_n := 0.125*exp((-v(t))/80.0 + (-1)*13.0/16.0);
alpha_m := 0.1*(v(t) + 40.0)/(1 - exp(-1.0/10.0*v(t) - 4.0));
beta_m := 4.0*exp(-1.0/18.0*v(t) - 65.0/18.0);
alpha_h := 0.07*exp((-1)*(v(t) + 65.0)/20.0);
beta_h := 1/(1.0 + exp((-v(t) + 35.0)/10.0));

dsys1 := {diff(h(t), t) = alpha_h*(1 - h(t)) - beta_h*h(t), diff(m(t), t) = alpha_m*(1 - m(t)) - beta_m*m(t), diff(n(t), t) = alpha_n*(1 - n(t)) - beta_n(t), diff(v(t), t) = III - gK*n(t)^4*(v(t) - VK) - gNa*m(t)^3*h(t)*(v(t) - VNa) - gL*(v(t) - VL), h(0) = 0.9996937394, m(0) = 0.02890553447, n(0) = 0.2445865495, v(0) = -70};

HPM.mw

  • Please help me
  • How do I write Homotopy Perturbation method for a partial differential equation for this question in Maple?

Download HPM.mw

 

modifed_practice.mw

Impact of Shape-Dependent Hybrid Nanofluid on Transient Efficiency of a Fully W
et Porous Longitudinal Fin

dear sir please help me to solve the graph i given reference pdf also. i have implimented the code but getting error in ploting 

Thank you

Using the SingularValues command with any output type specified causes this error

Error, (in sanity) invalid input: LinearAlgebra:-SingularValues expects value for keyword parameter output to be of type {list(identical(U,S,Vt,list)), identical(U,S,Vt,list)}, but received [U, S, Vt]

Below is a minimal example to trigger the error. Is there something wrong about how I am calling the function?

restart:
with(LinearAlgebra):

sanity := proc()
  local A,U,S,Vt;
  A :=RandomMatrix(3,10);
  U,S,Vt := SingularValues(A, output=['U','S','Vt']);
end proc:

sanity();

Dear Maple experts,

I have a system of several nonlinear equaitons. My code can solve it for a given parameter. But when I want to plot it, it takes too much time with no results. So, I decided to plot it for several given points. I get the answer for the points individually, but I don't know how to apply this to 'plot' command. Would you please help?

restartNULL

with(plots); with(RealDomain)

NULL

c := 1; cr := 0.3e-1*c; u := 5; sExp := 0.6e-1*c; s := .65*c; v := 3*c

NULL

NULL

FirmModelPartial1 := proc (alpha, beta, delta) local L1s, qs, ps, prs, hs, `&kappa;s`, `&lambda;__1s`, `&lambda;__2s`, `&lambda;__3s`, q, p, pr, h, kappa, `&lambda;__1`, `&lambda;__2`, `&lambda;__3`, FirmpfSiS, RecpfSiS, UnsoldSiS, EnvironSiS, p0, OldSoldPrim, xi, prof1, prof2, L1, L2, E1, E2, E3, E4, E5, E6, E7, E8, E9, E10, E11, E12; prof1 := (ps-c)*qs+((1/2)*beta^2*`&kappa;s`^2/(u*(1-alpha))-(1/2)*(qs+beta*`&kappa;s`)^2/u)*(ps-s)+hs*prs*(beta*`&kappa;s`-(1/2)*beta^2*`&kappa;s`^2/(u*(1-alpha))); L1s := prof1+`&lambda;__1s`*(1-sExp/prs-hs)+`&lambda;__2s`*(qs-`&kappa;s`)+`&lambda;__3s`*(qs-alpha*beta*`&kappa;s`/(1-alpha)); E1 := qs*(diff(L1s, qs)) = 0; E2 := hs*(diff(L1s, hs)) = 0; E3 := `&lambda;__1s`*(1-sExp/prs-hs) = 0; E4 := `&lambda;__2s`*(qs-`&kappa;s`) = 0; E5 := `&lambda;__3s`*(qs-alpha*beta*`&kappa;s`/(1-alpha)) = 0; E6 := qs = alpha*u*(v-ps)/(v-s); E7 := prs = ps-delta*v; E8 := `&kappa;s` = (beta*prs*(1-hs)+sExp*(1-beta))/(beta^2*(prs*(1-hs)-sExp)/(u*(1-alpha))+2*cr); p, q, pr, h, kappa, `&lambda;__1`, `&lambda;__2`, `&lambda;__3` := (eval([ps, qs, prs, hs, `&kappa;s`, `&lambda;__1s`, `&lambda;__2s`, `&lambda;__3s`], solve({0 <= qs-alpha*beta*`&kappa;s`/(1-alpha), 0 <= qs-`&kappa;s`, 0 <= 1-sExp/prs-hs, 0 <= `&lambda;__1s`, 0 <= `&lambda;__2s`, 0 <= `&lambda;__3s`, diff(L1s, qs) <= 0, diff(L1s, hs) <= 0, c < ps, (1/2)*beta*`&kappa;s`/(u*(1-alpha)) < 1, sExp+delta*v < ps, E1, E2, E3, E4, E5, E6, E7, E8}, [ps, qs, prs, hs, `&kappa;s`, `&lambda;__1s`, `&lambda;__2s`, `&lambda;__3s`])[1]))[]; xi := kappa/q; FirmpfSiS := max(0, eval(prof1, [ps = p, qs = q, prs = p-delta*v, hs = h, `&kappa;s` = kappa])); RecpfSiS := ((1-h)*pr-sExp)*(beta*kappa-(1/2)*beta^2*kappa^2/(u*(1-alpha)))+(sExp-cr)*kappa; UnsoldSiS := (1/2)*(q+beta*kappa)^2/u-(1/2)*beta^2*kappa^2/(u*(1-alpha)); EnvironSiS := q+UnsoldSiS; return p, q, FirmpfSiS, RecpfSiS, EnvironSiS, h, UnsoldSiS, h, xi end proc

NULL

NULL

FirmModelPartial1(.1, .2, .1)[3]

Error, (in FirmModelPartial1) invalid subscript selector

 

NULL

FirmModelPartial1(.2, .2, .1)[3]

Error, (in FirmModelPartial1) invalid subscript selector

 

FirmModelPartial1(.3, .2, .1)[3]

.2200894743

(1)

NULL

FirmModelPartial1(.4, .2, .1)[3]

.2892440018

(2)

FirmModelPartial1(.5, .2, .1)[3]

.3801240291

(3)

FirmModelPartial1(.6, .2, .1)[3]

.4921660779

(4)

FirmModelPartial1(.7, .2, .1)[3]

.6227094789

(5)

FirmModelPartial1(.8, .2, .1)[3]

.7668402612

(6)

FirmModelPartial1(.9, .2, .1)[3]

.9161921261

(7)
 

plot('FirmModelPartial1(alpha, .2, .2)[3]', alpha = .1 .. .90, adaptive = false, sample = [.1, .2, .3, .4, .5, .6, .7, .8, .9])

 

Download Reza_Tchno_(v3).mw

Besdies, it seems for some points there is no feasible solution. So, there is no need to have plot for them.

I attach the code. Thank you so much!

dS/dt=qnv-alpha_1*S*H-mu_sh*S

dH/dt=lamda_h-alpha_1*S*H-mu_h*H

dI/dt=alpha_1*S*H-beta_1*I-mu_i*I

dS/dt=beta_1*I-delta_1*S-mu_s*S

DATA:

q=0.75; n=15;v=20; alpha_You1=0.001; S=2000; H=3000; mu_sh=0.000000000012; lambda=3000;mu_h=0.94; beta_1=0.125; I=0; mi_i=0.95; delta_1=0.0975; S=0; mu_s=0.029

First 172 173 174 175 176 177 178 Last Page 174 of 2428