Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

hi

 

please help me :

 

 

restart; eq := diff(T(x, y), x) = a*(diff(T(x, y), `$`(y, 2)))/u(y); u := proc (y) options operator, arrow; (-1)*1.218493440*10^11*y^2+4.244913600*10^6*y+0.33e-1 end proc; eq; ICs := (D[1](T))(x, 0) = 1000, (D[1](T))(x, 0.25e-4) = 2000, T(0, y) = 0; T_sol := pdsolve({ICs, eq}, T(x, y)); T_sol

diff(T(x, y), x) = a*(diff(diff(T(x, y), y), y))/u(y)

 

proc (y) options operator, arrow; (-1)*1.218493440*100000000000*(y^2)+4.244913600*1000000*y+0.33e-1 end proc

 

diff(T(x, y), x) = a*(diff(diff(T(x, y), y), y))/(-0.1218493440e12*y^2+4244913.600*y+0.33e-1)

 

(D[1](T))(x, 0) = 1000, (D[1](T))(x, 0.25e-4) = 2000, T(0, y) = 0

 

Error, (in PDEtools:-Library:-NormalizeBoundaryConditions) unable to isolate the functions {T(0, y), (D[1](T))(x, 0), (D[1](T))(x, 0.25e-4)} in the given boundary conditions {T(0, y) = 0, (D[1](T))(x, 0) = 1000, (D[1](T))(x, 0.25e-4) = 2000}

 

T_sol

(1)

NULL

``

BC1 = diff(T(x, 0), y)=1000

BC2 = diff(T(x, 0.000025), y)=2000

IC = T(0,y)=0

where :

u(y)=-1.218493440*10^11*y^2+4.244913600*10^6*y+0.033

Download PDE_Sol.mw

Hi,

 

  I have a loop code, such as

 

****

i_max:=10;

for i from 1 to i_max do
  blah blah blah

end do;

****

 

  I would like to clean memory, something like restart suppose to do, after each cycle. Restart could only work at top level. How should I do to clean memory after each cycle?

 

Dear all,

It's very convenient to define a DE or PDE through Differential Operator D, for example,

((D[1, 1]+D[1, 2]+D[2, 2])(z))(x, y) = exp(x)*sin(y)

Is it possible to realize Inverse Operator Method of Operator D? How to solve the following equation if we rewrite the pde through inverse operator method?

(z)(x, y)=((D[1, 1]+D[1, 2]+D[2, 2])^(-1))exp(x)*sin(y)

 

Thanks a lot.

 

We are using your maple service and required your help in solving below problem '
1. Create a question type Maple Numeric with correct answer 27
2. If user gives 27d and 27e and on evaluate with your service getting response as true.

Please find the attached screens for more information. Also see request and response for numeric value 10 and if you see response even 10d is evaluating as true. This is wrong, please review and advice.

-------------------

mapleScript :: before substitution

expr1:="<math><cn>${#answer}</cn></math>"; expr2:="<math><cn>${#correctAnswer}</cn></math>"; expr3:="<math><cn>${#tolrence}</cn></math>";expr1 := StringTools[SubstituteAll](expr1, ",", "");expr1 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\1</ci>", expr1);expr1 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr1);expr1 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr1);expr1 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr1);expr2 := StringTools[SubstituteAll](expr2, ",", "");expr2 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr2);expr2 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr2);expr1 := MathML[ImportContent](expr1);expr2 := MathML[ImportContent](expr2);expr := combine(expr1-expr2, [radical,power,trig,ln,power,radical], symbolic);expr3:=MathML[ImportContent](expr3);result := evalb(abs(simplify(expr, symbolic))<=expr3);result := evalb(result = true);if not result then result := evalb(abs(simplify(expr, assume = real))<=expr3) end if;result := evalb(result = true);

mapleScript :: (user input) 10d

expr1:="<math><cn>10d</cn></math>"; expr2:="<math><cn>10</cn></math>"; expr3:="<math><cn>0</cn></math>";expr1 := StringTools[SubstituteAll](expr1, ",", "");expr1 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\1</ci>", expr1);expr1 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr1);expr1 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr1);expr1 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr1);expr2 := StringTools[SubstituteAll](expr2, ",", "");expr2 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr2);expr2 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr2);expr1 := MathML[ImportContent](expr1);expr2 := MathML[ImportContent](expr2);expr := combine(expr1-expr2, [radical,power,trig,ln,power,radical], symbolic);expr3:=MathML[ImportContent](expr3);result := evalb(abs(simplify(expr, symbolic))<=expr3);result := evalb(result = true);if not result then result := evalb(abs(simplify(expr, assume = real))<=expr3) end if;result := evalb(result = true);

evalResult ::
expr1 := "<math><cn>10d</cn></math>" expr2 := "<math><cn>10</cn></math>" expr3 := "<math><cn>0</cn></math>" expr1 := "<math><cn>10d</cn></math>" expr1 := "<math><cn>10d</cn></math>" expr1 := "<math><cn>10d</cn></math>" expr1 := "<math><cn>10d</cn></math>" expr1 := "<math><cn>10d</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr1 := 10. expr2 := 10 expr := 0. expr3 := 0 result := true result := true result := true

mapleScript :: (user input) 10
expr1:="<math><cn>10</cn></math>"; expr2:="<math><cn>10</cn></math>"; expr3:="<math><cn>0</cn></math>";expr1 := StringTools[SubstituteAll](expr1, ",", "");expr1 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\1</ci>", expr1);expr1 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr1);expr1 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr1);expr1 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr1);expr2 := StringTools[SubstituteAll](expr2, ",", "");expr2 := StringTools[RegSubs]("<apply>< *([a-zA-Z]+)[ /]*>" = "<apply><ci>\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<apply><ci>(plus|minus|divide|times|power)</ci>" = "<apply><\\1/>", expr2);expr2 := StringTools[RegSubs]("<apply><ci> *([a-zA-Z]+) *</ci>" = "<apply><ci>x_\\1</ci>", expr2);expr2 := StringTools[RegSubs]("<c[ni]> *([0-9]+\\.[0-9]+) *</c[ni]>" = "<apply><ci>convert</ci><cn>\\1</cn><ci>rational</ci><ci>exact</ci></apply>", expr2);expr1 := MathML[ImportContent](expr1);expr2 := MathML[ImportContent](expr2);expr := combine(expr1-expr2, [radical,power,trig,ln,power,radical], symbolic);expr3:=MathML[ImportContent](expr3);result := evalb(abs(simplify(expr, symbolic))<=expr3);result := evalb(result = true);if not result then result := evalb(abs(simplify(expr, assume = real))<=expr3) end if;result := evalb(result = true);

evalResult ::
expr1 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr3 := "<math><cn>0</cn></math>" expr1 := "<math><cn>10</cn></math>" expr1 := "<math><cn>10</cn></math>" expr1 := "<math><cn>10</cn></math>" expr1 := "<math><cn>10</cn></math>" expr1 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr2 := "<math><cn>10</cn></math>" expr1 := 10 expr2 := 10 expr := 0 expr3 := 0 result := true result := true result := true

When I try to export my data from a spreadsheet to Excel, instead of numbers I get a long text containing numbers. For example instead of having just "11625.7677908435" I get "Typesetting:-mprintslash([11625.7677908435], [HFloat(11625.767790843516)])". 

I want to get numerical solution of the Eqs.ode(see the folowlling ode and ibc)in Maple.However,when i run the following procedure,it prompts an error "Error, (in dsolve/numeric/bvp) cannot determine a suitable initial profile, please specify an approximate initial solution". How to solve the issue? Please help me.


restart:
n := 1.4; phi := 1; beta := .6931; psi := 1

> restart;
> n := 1.4; phi := 1; beta := .6931; psi := 1;

> s := proc (x) options operator, arrow; evalf(1+(phi*exp(beta*psi)*h(x))^n) end proc;

> Y := proc (x) options operator, arrow; evalf(f-(1/2-(1/2)/n)*ln(s(x))+2*ln(1-(1-s(x))^(-1+1/n))) end proc;


> ode := diff(h(x), `$`(x, 2))+(diff(Y(x), x))*(diff(h(x), x)+1) = 0;


> ibc := h(0) = 0, ((D(h))(10)+1)*s(10)^(-(1-1/n)*(1/2))*(1-(1-1/s(10))^(1-1/n))^2 = 0;

> p := dsolve({ibc, ode}, numeric);
Error, (in dsolve/numeric/bvp) cannot determine a suitable initial profile, please specify an approximate initial solution
>

Hi Maple friends.

expand( (a+b)^2 );

a^2+2*a*b+b^2

expand( (a+b)^3 );

a^3+3*a^2*b+3*a*b^2+b^3

expand( (a+4)^4 )

a^4+16*a^3+96*a^2+256*a+256 (???)

Pascal's triangle shows that 'expand( (a+4)^4 )' should have resulted in a^4+4*a^3*b+6*a^2*b^2+4*a*b^3+b^4

Where are the b variables in Maple's solution?

Thanks in advance.


``

restart:

l := 16; m := .75; EIy := 0.2e5; EIz := 0.400e7; GJ := 0.1e5; mj := .1; Pz := 5000;

16

 

.75

 

0.2e5

 

0.400e7

 

0.1e5

 

.1

 

5000

(1)

eq1 := m*(diff(diff(w(x, t), t), t))+diff(diff(EIy*(diff(diff(w(x, t), x), x))+(EIz-EIy)*(diff(diff(w(x, t), x), x))*(theta(x)+phi(x, t))^2, x), x)-Pz = 0;

.75*(diff(diff(w(x, t), t), t))+0.2e5*(diff(diff(diff(diff(w(x, t), x), x), x), x))+0.398e7*(diff(diff(diff(diff(w(x, t), x), x), x), x))*(theta(x)+phi(x, t))^2+0.1592e8*(diff(diff(diff(w(x, t), x), x), x))*(theta(x)+phi(x, t))*(diff(theta(x), x)+diff(phi(x, t), x))+0.796e7*(diff(diff(w(x, t), x), x))*(diff(theta(x), x)+diff(phi(x, t), x))^2+0.796e7*(diff(diff(w(x, t), x), x))*(theta(x)+phi(x, t))*(diff(diff(theta(x), x), x)+diff(diff(phi(x, t), x), x))-5000 = 0

(2)

w := unapply(q[1](t)*cosh(1.8751*x/l),x, t);

proc (x, t) options operator, arrow; q[1](t)*cosh(.1171937500*x) end proc

(3)

theta := proc (x) options operator, arrow; 0 end proc;

proc (x) options operator, arrow; 0 end proc

(4)

#theta(x);eq1;

phi := unapply(q[4](t)*sin(1.5708*x/l),x, t);

proc (x, t) options operator, arrow; q[4](t)*sin(0.9817500000e-1*x) end proc

(5)

eq2:=unapply(lhs(eq1),x,t);

proc (x, t) options operator, arrow; .75*(diff(diff(q[1](t), t), t))*cosh(.1171937500*x)+3.772661154*q[1](t)*cosh(.1171937500*x)-302.9569534*q[1](t)*cosh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)^2+2515.691179*q[1](t)*sinh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)*cos(0.9817500000e-1*x)+1053.716523*q[1](t)*cosh(.1171937500*x)*q[4](t)^2*cos(0.9817500000e-1*x)^2-5000 end proc

(6)

f[1] := int(expand(eq2(x,t)*cosh(1.8751*x/l)), x = 0 .. l);

int(.75*(diff(diff(q[1](t), t), t))*cosh(.1171937500*x)^2+3.772661154*q[1](t)*cosh(.1171937500*x)^2-302.9569534*q[1](t)*cosh(.1171937500*x)^2*q[4](t)^2*sin(0.9817500000e-1*x)^2+2515.691179*cosh(.1171937500*x)*q[1](t)*sinh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)*cos(0.9817500000e-1*x)+1053.716523*q[1](t)*cosh(.1171937500*x)^2*q[4](t)^2*cos(0.9817500000e-1*x)^2-5000*cosh(.1171937500*x), x = 0 .. 16)

(7)

 

 

``


Download 1.mw

Hello,

I would like to open a modelica model with MapleSIM. When I do "import modelica", I received the error message : "import failed : no library found".

I believed the modelica library was already installed on maplesim but may be not. Do I need to install modelica library ? and if yes, how ?

More generally speaking, how can I do to open a modelica model on maplesim?

Thank you for your help.

Hi, I would like to know if there is a command on Maple 18 that can show all the decompositions of a given symbolic tensor.

 

 

 

Hallo,

how can I  construct a Permutationgroup with given GroupOrder and Degree.

for instance: GroupOrder 12 and Degree 5, I found by accident

G:=<(1,2,3),(1,5,2)>

Best Regards

 

Kurt Ewald

Hi:

how i can plot the function in terms of different values of n in a figure:

f(x)=1+2*(1+x)^n

n=1,2,3,4

x-axis->x

y-axis->f(x)

 

Good Afternoon Everyone,

Would anybody be kind to help me with the interpolation of two variables in a table to obtain a value for f18.1 (Stiffness) and f18.2 (Damping) given SR_ratio and Lamda_ratio.

 

Regards

Moses


 

NULL

restart

with(CurveFitting)

with(ExcelTools):

``

Slenderness Ratio (Variable 1):

SR__ratio := 84.01

``

Ratio of modulus of elasticity of Pile vs shear modulus of soil (Variable 2):

`&lambda;__ratio` := 229.58

``

Slenderness ratio of pile from stiffness of pile table:

SR__stiff := [10.8696, 21.7391, 32.6087, 43.4783, 46.7391, 54.3478, 65.2174, 76.0870, 86.9565, 100.0]

``

Ratio of Young's Modulus of Pile to Shear modulus of soil:

EG__ratio := [250.0, 500.0, 1000.0, 2500.0, 5000.0, 10000.0]

``

Vertical stiffness coefficient for floating pile as per Novak (1983) f18.1.250:

`f__18.1.250` := [0.332e-1, 0.509e-1, 0.571e-1, 0.582e-1, 0.582e-1, 0.582e-1, 0.582e-1, 0.582e-1, 0.582e-1, 0.582e-1]

``

Vertical stiffness coefficient for floating pile as per Novak (1983) f18.1.500:

`f__18.1.500` := [0.187e-1, 0.301e-1, 0.364e-1, 0.405e-1, 0.416e-1, 0.416e-1, 0.416e-1, 0.416e-1, 0.416e-1, 0.416e-1]

``

Vertical stiffness coefficient for floating pile as per Novak (1983) f18.1.1000:

`f__18.1.1000` := [0.104e-1, 0.166e-1, 0.218e-1, 0.26e-1, 0.27e-1, 0.281e-1, 0.291e-1, 0.301e-1, 0.301e-1, 0.301e-1]

``

Vertical stiffness coefficient for floating pile as per Novak (1983) f18.1.2500:

`f__18.1.2500` := [0.52e-2, 0.83e-2, 0.104e-1, 0.125e-1, 0.135e-1, 0.145e-1, 0.166e-1, 0.177e-1, 0.187e-1, 0.197e-1]

``

Vertical stiffness coefficient for floating pile as per Novak (1983) f18.1.10000:

`f__18.1.10000` := [0.21e-2, 0.31e-2, 0.42e-2, 0.42e-2, 0.52e-2, 0.52e-2, 0.62e-2, 0.62e-2, 0.73e-2, 0.83e-2]

``

Slenderness ratio of pile from damping of pile table:

SR__damp := [10.8696, 16.3043, 21.7391, 27.1739, 32.6087, 38.0435, 43.4783, 48.913, 54.3478, 59.7826, 65.2174, 70.6522, 76.0870, 81.5217, 86.9565, 92.3913, 100.0]

``

Vertical damping coefficient for floating pile as per Novak (1983) f18.2.250:

`f__18.2.250` := [.1032, .1137, .1126, .1095, .1053, .1021, 0.989e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.979e-1, 0.989e-1, 0.989e-1]

``

Vertical damping coefficient for floating pile as per Novak (1983) f18.2.500:

`f__18.2.500` := [0.558e-1, 0.695e-1, 0.811e-1, 0.832e-1, 0.811e-1, 0.789e-1, 0.758e-1, 0.737e-1, 0.726e-1, 0.716e-1, 0.705e-1, 0.695e-1, 0.695e-1, 0.695e-1, 0.695e-1, 0.705e-1, 0.705e-1]

``

 Vertical damping coefficient for floating pile as per Novak (1983) f18.2.1000:

`f__18.2.1000` := [0.295e-1, 0.421e-1, 0.495e-1, 0.537e-1, 0.568e-1, 0.589e-1, 0.579e-1, 0.568e-1, 0.558e-1, 0.537e-1, 0.526e-1, 0.516e-1, 0.516e-1, 0.505e-1, 0.505e-1, 0.505e-1, 0.495e-1]

``

 Vertical damping coefficient for floating pile as per Novak (1983) f18.2.2500:

`f__18.2.2500` := [0.126e-1, 0.179e-1, 0.232e-1, 0.263e-1, 0.305e-1, 0.326e-1, 0.347e-1, 0.368e-1, 0.379e-1, 0.379e-1, 0.379e-1, 0.379e-1, 0.379e-1, 0.368e-1, 0.358e-1, 0.358e-1, 0.337e-1]

``

   Vertical damping coefficient for floating pile as per Novak (1983) f18.2.10000:

`f__18.2.10000` := [0.32e-2, 0.53e-2, 0.74e-2, 0.84e-2, 0.105e-1, 0.116e-1, 0.137e-1, 0.147e-1, 0.147e-1, 0.168e-1, 0.168e-1, 0.179e-1, 0.189e-1, 0.189e-1, 0.200e-1, 0.211e-1, 0.211e-1]

``

 

NULL

``

``

`f__18.1` := ArrayInterpolation(SR__ratio, `&lambda;__ratio`, SR__stiff, EG__ratio, `f__18.1.250`)

`f__18.2__ _` := ArrayInterpolation(SR__ratio, `&lambda;__ratio`, SR__damp, EG__ratio, `f__18.2.250`)

NULL

NULL


Download Dynamic_Table.mwDynamic_Table.mw

Good Afternoon Everyone,

Would anybody be kind to help me with the interpolation of two variables in a table to obtain a value for f18.1 (Stiffness) and f18.2 (Damping) given SR_ratio and Lamda_ratio.

 

Regards

Moses

Hi:

how i can plot the function f(x) in term of 1/x in maple?for example i have the function f(x)=x^2 and x=0..1,now i plot f(x) in term of 1/x.

It is possible to add groups of questions to an assignment in Maple T.A. But how do you see the grades divided into these groups? As an example, suppose we have two groups of questions say 5 questions in the group "algebra" and 5 questions in the group "geometry". The class grades show the grades for all 10 questions all together, but I would like to see the grades for two groups individually.

 

 

First 256 257 258 259 260 261 262 Last Page 258 of 362