## Plotting pseudospectra...

Hi everyone,

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

Best,

Rachid

## What will be the range of p and q to get the plot ...

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

## How to find the value of constants?...

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

## How can I insert L'(2) =0 and L''(2) > 0 in this l...

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}]]

## How can I extract all the nonzero powers from poly...

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?

## Need to Solve System of ODE Using RK-4...

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};

## How do I write Homotopy Perturbation method for a ...

HPM.mw

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

## 1 dependent and 2 independent variables....

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

## Error using SingularValues in a procedure...

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();

## Point plot for system of nonlinear equations...

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?

 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 (1)
 >
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)

 >

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!

## How do I stimulate the nonlinear models in maple?...

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

## Newton iteration not converging...

Good day all,

please I am new Maple user, i used the same code to solve boundary layer flow. Itried to use the same code in solving nanofluid but the solution is not coverging at some point. But the author i want to compared my results with had solution at those points.

Is there any way i can get solution of this model at this with these values?. Please i need your help. Thank you in advance

secod_code.mw

 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >

 >
 (5)
 >
 >
 > for k to 10 do       sol_All := dsolve                  ( eval                    ( {bcs1, ode1, ode2,ode3},                      [xi= L[k],R=1, CT=0.2, Gr=0.5, Gc=0.5, lambda1=0.2, lambda2=0.2, Pr=0.72, Nb=0.5,  Le=1, dc=0.2, de=0.2, Nt=0.5, Zeta=0.5, varsigma=0.5, M=1.5, bet=1.5]                    ),                    [f(eta), theta(eta),phi(eta)],                    numeric,                    output = listprocedure                  );       Y_sol || k := rhs(sol_All[5]);       YP_sol || k := -rhs(sol_All[6]); feta || k := rhs(sol_All[4]);       fpeta || k := rhs(sol_All[3]); phieta || k := rhs(sol_All[8]);       phipeta || k := rhs(sol_All[7])   end do:
 >
 >
 >
 >
 > plot( [ seq((Y_sol||j)(eta), j = 1..6)],          eta = 0 .. 6,          labels = [eta, theta(eta)],          axes = boxed       ); plot( [ seq((YP_sol||j)(eta), j = 1..6)],          eta = 0 .. 8,          labels = [eta, thetaprime(eta)],          axes = boxed       );  plot( [ seq((feta||j)(eta), j = 1..6)],          eta = 0 .. 8,          labels = [eta, f(eta)],          axes = boxed       );   plot( [ seq((fpeta||j)(eta), j = 1..6)],          eta = 0 .. 5,          labels = [eta, fprime(eta)],          axes = boxed       ); plot( [ seq((phipeta||j)(eta), j = 1..6)],          eta = 0 .. 5,          labels = [eta, phiprime(eta)],          axes = boxed       ); plot( [ seq((phieta||j)(eta), j = 1..6)],          eta = 0 .. 5,          labels = [eta, phi(eta)],          axes = boxed       );

## How to protect/unprotect parameters...

Sometimes the canvas in the 2d-workspace of a subsystem (or component) is grayed out and the parameters in the parameter pane cannot be changed.

Example from the Model Gallery where the mass cannot be changed

Another example of a CAD geometry where I would like to change color and transparency:  seeshell.msim

How to unprotect the parameters?

I have observed that copy and pasting (Edit2: in a new document) and deleting the original component and copy/paste back into the original document and reconnection of the pasted  component works. Are there shorter ways to unprotect?

And: How to do the reverse action (i.e. protecting a component)?

(I could not find a "when" and "why" for the protection in the documentation. Maybe I overlooked something)

Edit:

Other reasons why to unprotect are suppressions of components and highlight of multibody ports and axes.

Adding or reformating automatically generated ports (Edit3).

## Analytical Solution of Nonlinear ODE ...

Hello Everyone;

Can anybody help me to solve the analytical solution of the following ODE?
AAS := C*diff(y(t), t) + (-B0*y(t)^3 - B1*y(t)^2 - B2*y(t) - B3);y[0] := 1

C,B0,B1,B2 are constants.

## A resistant integro-differential ode....

How did Maple arrive at the following, implicit solution,

MSimplicitesolution := (-9*y(x)*(-y(x)*x^2)^(1/3)*6^(2/3)*x^2 + 16*_C1*x^(8/3) - 24*x^2 + 6)/(16*x^(8/3)) = 0,

to the following 'ode',

ODE:= (2/3)*(int[(y'(x)*(x^2)/((x^2) -1)]) =int(-sqrt [2*y(x)])?

'odeadvisor' suggested y=G(x,y'(x), but I could not see how this could be implimented with this equation.

(Unfortunately, I am unable to download the Maple worksheet onto this sheet at this time.  Any assitance

would be appreciated.

 First 88 89 90 91 92 93 94 Last Page 90 of 2344
﻿