Unanswered Questions

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

Here, I have a 3D map T=(T1,T2,T3) with

T1(x1,x2,x3)=x1*exp(6-2*x1-x2-0.3*x3),

T2(x1,x2,x3)=x2*exp(5-2*x1-2*x2-*x3),

T3(x1,x2,x3)=x3*exp(8-0.2*x1-x2-x3),

How can I use  IterativeMaps:-Attractor obtain the attractor for T.

Hello i need a tutor to help me with  Numerical analysis on Maple and the programming part. Im In UAE but we can do this online as well

Dear All,

The excel file consists of data (v2, Re, t) and I would like plot the variable "t" aginst the variables "v2" and "Re". Could anyone point me out?

Attached is the excel data.Collapsetime.xlsx 

Thank you.

Very kind wishes,

Wang Zhe

I am not unfamiliar with the Wolfram syntax but also not very good with it, and there is a particular element in a Mathematica code I have been given which I do not entirely understand how to efficiently write in Maple. The basic idea is to read in a list of expressions from an external file (LIST) and process the non zero elements and assign them to a function (COEF) which can be called later on. Here is the Mathematica exert:

k = 0;
i = 0;
a = b = \[Theta];
Do[k = k + 1; KK = LIST[[k]]; 
  If[KK =!= 0, i = i + 1; ff = Factor[KK]; 
   COEF[x,y, z, l_, m_, n_] = ff], {z, -2, 
   2}, {y, -2, 2}, {x, -2, 2}];

The LIST has the following form and only contains l, m and n and another factor E which is left undefined for now. It does not contain x, y or z. The LIST can contain any number of terms depending on the problem. Here is an example:

LIST={0, 0, 0, 0, 0, 0, 0, a^2 b m (-1 + n) n (a^2 + b^2 - 2 E), ... ,0,0, a^3 n(l+1+m) ... }

So the Do loop cycles through the LIST and extracts out the non zero terms. What I am unsure about is how it is looping over x,y and z when they do not appear in the LIST at all. I assume it is attaching a x,y,z combination to each COEF and they can be called like this:

COEF[0,1,1,0,2,3]

For the instance of when x=0, y=1, z=1, l=0, m=2 n=3. Is this correct? What would be the best way to replicate this in Maple?

- Yeti

i want to solve the system of equation ( 1 )  , (2)  ,  (3)   under the assumation that x , y have the CDF in (4)  ,  (5)
 

diff(L(lambda[1], lambda[2], alpha), lambda[1]) = n/lambda[1]+sum(x[i], i = 1 .. n)-(sum(2*x[i]*exp(lambda[1])/(exp(x__i*`λ__1`)-1+alpha), i = 1 .. n))

diff(L(lambda[1], lambda[2], alpha), lambda[1]) = n/lambda[1]+sum(x[i], i = 1 .. n)-(sum(2*x[i]*exp(lambda[1])/(exp(x__i*`λ__1`)-1+alpha), i = 1 .. n))

(1)

diff(L(lambda[1], lambda[2], alpha), lambda[2]) = m/lambda[2]+sum(y[j], j = 1 .. m)-(sum(2*y[j]*exp(lambda[2])/(exp(y__j*`λ__2`)-1+alpha), j = 1 .. m))

diff(L(lambda[1], lambda[2], alpha), lambda[2]) = m/lambda[2]+sum(y[j], j = 1 .. m)-(sum(2*y[j]*exp(lambda[2])/(exp(y__j*`λ__2`)-1+alpha), j = 1 .. m))

(2)

diff(L(lambda[1], lambda[2], alpha), alpha) = (n+m)/alpha-(sum(2/(exp(x[i]*`λ__1`)-1+alpha), i = 1 .. n))-(sum(2/(exp(y[j]*`λ__2`)-1+alpha), j = 1 .. m))

diff(L(lambda[1], lambda[2], alpha), alpha) = (n+m)/alpha-(sum(2/(exp(x[i]*`λ__1`)-1+alpha), i = 1 .. n))-(sum(2/(exp(y[j]*`λ__2`)-1+alpha), j = 1 .. m))

(3)

G(x, lambda[1], alpha) = 1-alpha/(exp(lambda[1]*x)-1+alpha)

G(x, lambda[1], alpha) = 1-alpha/(exp(lambda[1]*x)-1+alpha)

(4)

G(y, lambda[2], alpha) = 1-alpha/(exp(lambda[2]*x)-1+alpha)

G(y, lambda[2], alpha) = 1-alpha/(exp(lambda[2]*x)-1+alpha)

(5)

``

``


 

Download internet.mw

How do we write code for optimal problem using Pontryagin's maximum principle for simulation.

I want to make the blue output my procedure spits out a another color, and align it to the right, is this even possible? Or something like it?

I basically want to make a Maple procedure that does certain calculations and writes the explanation for each calculation. I do however want Maple to write these explanations as a text field like in a normal Maple worksheet, instead of the blue output in the middle. Is this possible? Or is there any alternative ideas you have that I could try? Would really appreciate any kind of help, thanks.

I have included a transfer function of a 3 degree of freedom system.  There are special loci, aside from s = 0, where the numerator will equate to a REAL scalar value, and the denominator will reduce to a simple product of the individual z roots.  Does anyone know if there is some physical significance, interpretation, or analog to this condition and special loci?untitled.mw

 

hello, i just try to plot the relation between my outputs (u, and phat) with i from 0 to 10 , but i have aproblem any suggestions? 
 

M := .4556;

.4556

(1)

K := 18;

18

(2)

c := .2865;

.2865

(3)

Nabla(t) := .1NULL

.1

(4)

Khat := 206.1055;

206.1055

(5)

NULL

N := 10NULL

10

(6)

``

NULL

NULL

a__1 := 4/.1^2*.4556+2/(.1)*.2865NULL

187.9700000

(7)

``

NULL

NULL

a__2 := 4/(.1)*.4556+.2865NULL

18.51050000

(8)

``

NULL

NULL

a__3 := .4556NULL

.4556

(9)

NULL

NULLNULL

fu := Array(0 .. 10):

p:=Array(0..10):
  p[0]:=0:
  for i from 0 to 4 do
      p[i+1]:=50*sin(3.14*(i+1)*(4)/0.6):
      phat[i+1]:= p[i+1]+((7)*u[i])+((8)*u__dot[i])+((9)*u__doubledot[i]):
      u[i+1]:= phat[i+1]/(5):
      u__dot[i+1]:=(20*(u[i+1]-u[i]))-u__dot[i]:
      u__doubledot[i+1]:= ((400*(u[i+1]-u[i]))-(40*u__dot[i])-(u__doubledot[i])):
end do;
for i from 5 to 9 do
      p[i+1]:=0.0:
      phat[i+1]:= p[i+1]+((7)*u[i])+((8)*u__dot[i])+((9)*u__doubledot[i]):
      u[i+1]:= phat[i+1]/(5):
      u__dot[i+1]:=(20*(u[i+1]-u[i]))-u__dot[i]:
      u__doubledot[i+1]:= ((400*(u[i+1]-u[i]))-(40*u__dot[i])-(u__doubledot[i])):
  end do;

 

 

24.98850513

 

24.98850513

 

.1212413309

 

2.424826618

 

48.49653236

 

43.28799198

 

133.0574982

 

.6455795610

 

8.061937982

 

64.24569494

 

49.99998414

 

349.8504158

 

1.697433673

 

12.97514426

 

34.01843056

 

43.32779001

 

618.0696023

 

2.998802081

 

13.05222390

 

-32.47683816

 

25.05744793

 

815.5490181

 

3.956949320

 

6.11072088

 

-106.3532218

 

0.

 

808.4457346

 

3.922485012

 

-6.80000704

 

-151.8613364

 

0.

 

542.2499525

 

2.630933927

 

-19.03101466

 

-92.7588160

 

0.

 

100.0021368

 

.4851987783

 

-23.88368831

 

-4.2946573

 

0.

 

-352.8528440

 

-1.712001106

 

-20.06030938

 

80.7622360

 

0.

 

-656.3359300

 

-3.184465868

 

-9.38898586

 

132.6642346

(10)

 

``

fd_table:=eval(seq[u(i),phat(i)],i=0..N);

seq[(table( [( 0 ) = 0, ( 1 ) = .1212413309, ( 2 ) = .6455795610, ( 3 ) = 1.697433673, ( 4 ) = 2.998802081, ( 5 ) = 3.956949320, ( 6 ) = 3.922485012, ( 7 ) = 2.630933927, ( 9 ) = -1.712001106, ( 8 ) = .4851987783, ( 10 ) = -3.184465868 ] ))('i'), (table( [( 1 ) = 24.98850513, ( 2 ) = 133.0574982, ( 3 ) = 349.8504158, ( 4 ) = 618.0696023, ( 5 ) = 815.5490181, ( 6 ) = 808.4457346, ( 7 ) = 542.2499525, ( 9 ) = -352.8528440, ( 8 ) = 100.0021368, ( 10 ) = -656.3359300 ] ))('i')]

(11)

``

plot([u(i+1), p(i+1)])

Error, (in plot) invalid input: assigned expects its 1st argument, n, to be of type assignable, but received table( [( 0 ) = 0, ( 1 ) = .1212413309, ( 2 ) = .6455795610, ( 3 ) = 1.697433673, ( 4 ) = 2.998802081, ( 5 ) = 3.956949320, ( 6 ) = 3.922485012, ( 7 ) = 2.630933927, ( 9 ) = -1.712001106, ( 8 ) = .4851987783, ( 10 ) = -3.184465868 ] )

 

``

``


 

Download hw_4_structural.mw

Hello everyone,

 

     I am having trouble trying to solve a system of differential equations. The modeling was made from the equilibrium equations of a pressure vessel. The set of equations is shown below:

     As you see it is a set of two second-order partial differential equations. So, we need four boundary conditions. This one is the first. It means that the left end of the pressure vessel is fixed.

This one is the second boundary condition. It means that the right end of the pressure vessel is free.

This one is the third boundary condition. It means that the inner surface of the pressure vessel is subject to an external load:

At last, we have the fourth boundary condition. It means that the outer surface of the pressure vessel is free.

     The first test I have been trying to do is the static case. In this case, the time terms (the right side of the two equations shown) is zero.

    The maple commands that I am using are the following:

 

restart; E := 200*10^9; nu := .33; G := E/(2*(1+nu)); RI := 0.254e-1; RO := 2*RI; p := proc (x) options operator, arrow; 50000000 end proc; sys := [E*(nu*(diff(v(x, r), x))/r+nu*(diff(diff(v(x, r), x), r))+(1-nu)*(diff(diff(u(x, r), x), x)))/(-2*nu^2-nu+1)+G*(diff(diff(u(x, r), r), r)+diff(diff(v(x, r), x), r)+(diff(u(x, r), r))/r+(diff(v(x, r), x))/r) = 0, E*((1-nu)*(diff(diff(v(x, r), r), r))+nu*(diff(diff(u(x, r), x), r))+(1-nu)*(diff(v(x, r), r))/r-(1-nu)*v(x, r)/r^2)/(-2*nu^2-nu+1)+G*(diff(diff(u(x, r), r), x)+diff(diff(v(x, r), x), x)) = 0]; BCs := {E*(nu*v(L, r)/r+nu*(D[2](v))(L, r)+(1-nu)*(D[1](u))(L, r))/(-2*nu^2-nu+1) = 0, E*(nu*v(x, RI)/RI+(1-nu)*(D[2](v))(x, RI)+nu*(D[1](u))(x, RI))/(-2*nu^2-nu+1) = -p(x), E*(nu*v(x, RO)/RO+(1-nu)*(D[2](v))(x, RO)+nu*(D[1](u))(x, RO))/(-2*nu^2-nu+1) = 0, u(0, r) = 0}

sol := pdsolve(sys, BCs, numeric)

 

I am getting the following error:

 

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions must depend upon exactly one of the independent variables: 0.1459531181e12*v(L, r)/r+0.1459531181e12*(D[2](v))(L, r)+0.2963290579e12*(D[1](u))(L, r) = 0

In this case, my boundary conditions do depend on more than one independent variable. How do I proceed?

 

Thank you in advance,

Pedro Guaraldi

 

 

Is there anyone who has seen maple 2017 provide some details about what new features are being introduced. Is there a platform where we can suggest what features we would like to be added or enhanced?

I'm trying to define some multilinear forms to study differential geometry. What I need is only symbolic. My intention is symplify computations involving multilinear forms.

For example, to create an symbolic inner product "g" I used the command "define" like in this post:

http://www.mapleprimes.com/questions/203480-Define-And-Use-Abstract-Linear-Operator

 

So I tipped:

define(g, orderless, multilinear);

 

My doubt is: how can I declare that g(x,y) is always scalar?

With it I would simplify things like g(g(z,w)*x,y) = g(z,w)*g(x,y)

 

In my case, specifically, I type:

v:=(X,Y,Z)->g(Y,Z)*X-g(X,Z)*Y;

r:=(X,Y,Z,W)->g(v(X,Y,Z),W)-g(Y,T)*g(v(X,T,Z),W)+g(X,T)*g(v(Y,T,Z),W);

expand(r(X,Y,Z,W));

and the result is:

g(W,g(Y,Z)*X)-g(W,g(X,Z)*Y)-g(T,Y)*g(W,g(T,Z)*X)+g(T,Y)*g(W,g(X,Z)*T)+g(T,X)*g(W,g(T,Z)*Y)-g(T,X)*g(W,g(Y,Z)*T)

But I would enjoy that it were:

g(Y,Z)*g(W,X)-g(X,Z)*g(W,Y)-g(T,Z)*g(T,Y)*g(W,X)+g(X,Z)*g(T,Y)*g(W,T)+g(T,Z)*g(T,X)*g(W,Y)-g(Y,Z)*g(T,X)*g(W,T)

 

Is there a way to declare that g(x,y) is always scalar?

Thanks.

 

 

I have a problem using dchange when my variable depend on two (or more variables) and I would like to apply the chain rule.

For example, when I use the command

I would expect something like 

But I get an error saying that the number of new variables and transformation equations must be the same.

Any idea how I could solve it? 

Thanls a lot for your help.

 

i need to solve for u[i+1] as i attached i wrote the equations but i cant get any answers for it, the delta t is 0.1 and i need to go for ten steps, thank you
 

M := .4556;

.4556

(1)

K := 18;

18

(2)

c := .2865;

.2865

(3)

`u__ double dot`[0] := 0;

0

(4)

u__[0] := 0;

0

(5)

P__[0] := 0;

0

(6)

Typesetting:-delayGradient(t) := .1;

.1

(7)

N := 10;

10

(8)

a__1 := 4/.1^2*.4556+2/(.1)*.2865;

187.9700000

(9)

a__2 := 4/(.1)*.4556+.2865;

18.51050000

(10)

a__3 := .4556;

.4556

(11)

khat := 18+187.9700000;

205.9700000``

(12)

`u__ dot`[0] := 0;

0

(13)

 

for i from 0 to 10 do phat[i+1] := p[i+1]+187.9700000*u[i]+18.51050000*u__dot[i]+.4556*`u__ double dot`[i] end do

p[1]+187.9700000*u[0]+18.51050000*u__dot[0]

 

p[2]+187.9700000*u[1]+18.51050000*u__dot[1]+.4556*`u__ double dot`[1]

 

p[3]+187.9700000*u[2]+18.51050000*u__dot[2]+.4556*`u__ double dot`[2]

 

p[4]+187.9700000*u[3]+18.51050000*u__dot[3]+.4556*`u__ double dot`[3]

 

p[5]+187.9700000*u[4]+18.51050000*u__dot[4]+.4556*`u__ double dot`[4]

 

p[6]+187.9700000*u[5]+18.51050000*u__dot[5]+.4556*`u__ double dot`[5]

 

p[7]+187.9700000*u[6]+18.51050000*u__dot[6]+.4556*`u__ double dot`[6]

 

p[8]+187.9700000*u[7]+18.51050000*u__dot[7]+.4556*`u__ double dot`[7]

 

p[9]+187.9700000*u[8]+18.51050000*u__dot[8]+.4556*`u__ double dot`[8]

 

p[10]+187.9700000*u[9]+18.51050000*u__dot[9]+.4556*`u__ double dot`[9]

 

p[11]+187.9700000*u[10]+18.51050000*u__dot[10]+.4556*`u__ double dot`[10]

(14)

for i from 0 to 10 do u[i+1] := (1/18)*phat[i+1] end do;

(1/18)*p[1]+10.44277778*u[0]+1.028361111*u__dot[0]

 

(1/18)*p[2]+.5801543211*p[1]+109.0516077*u[0]+10.73894656*u__dot[0]+1.028361111*u__dot[1]+0.2531111111e-1*`u__ double dot`[1]

 

(1/18)*p[3]+.5801543211*p[2]+6.058422650*p[1]+1138.801706*u[0]+112.1444325*u__dot[0]+10.73894656*u__dot[1]+.2643183086*`u__ double dot`[1]+1.028361111*u__dot[2]+0.2531111111e-1*`u__ double dot`[2]

 

(1/18)*p[4]+.5801543211*p[3]+6.058422650*p[2]+63.26676144*p[1]+11892.25315*u[0]+1171.099388*u__dot[0]+112.1444325*u__dot[1]+2.760217359*`u__ double dot`[1]+10.73894656*u__dot[2]+.2643183086*`u__ double dot`[2]+1.028361111*u__dot[3]+0.2531111111e-1*`u__ double dot`[3]

 

(1/18)*p[5]+.5801543211*p[4]+6.058422650*p[3]+63.26676144*p[2]+660.6807306*p[1]+124188.1569*u[0]+12229.53067*u__dot[0]+1171.099388*u__dot[1]+28.82433650*`u__ double dot`[1]+112.1444325*u__dot[2]+2.760217359*`u__ double dot`[2]+10.73894656*u__dot[3]+.2643183086*`u__ double dot`[3]+1.028361111*u__dot[4]+0.2531111111e-1*`u__ double dot`[4]

 

(1/18)*p[6]+.5801543211*p[5]+6.058422650*p[4]+63.26676144*p[3]+660.6807306*p[2]+6899.342050*p[1]+1296869.325*u[0]+127710.2711*u__dot[0]+12229.53067*u__dot[1]+301.0061407*`u__ double dot`[1]+1171.099388*u__dot[2]+28.82433650*`u__ double dot`[2]+112.1444325*u__dot[3]+2.760217359*`u__ double dot`[3]+10.73894656*u__dot[4]+.2643183086*`u__ double dot`[4]+1.028361111*u__dot[5]+0.2531111111e-1*`u__ double dot`[5]

 

(1/18)*p[7]+.5801543211*p[6]+6.058422650*p[5]+63.26676144*p[4]+660.6807306*p[3]+6899.342050*p[2]+72048.29583*p[1]+13542918.17*u[0]+1333649.981*u__dot[0]+127710.2711*u__dot[1]+3143.340237*`u__ double dot`[1]+12229.53067*u__dot[2]+301.0061407*`u__ double dot`[2]+1171.099388*u__dot[3]+28.82433650*`u__ double dot`[3]+112.1444325*u__dot[4]+2.760217359*`u__ double dot`[4]+10.73894656*u__dot[5]+.2643183086*`u__ double dot`[5]+1.028361111*u__dot[6]+0.2531111111e-1*`u__ double dot`[6]

 

(1/18)*p[8]+.5801543211*p[7]+6.058422650*p[6]+63.26676144*p[5]+660.6807306*p[4]+6899.342050*p[3]+72048.29583*p[2]+752384.3428*p[1]+141425684.9*u[0]+13927010.38*u__dot[0]+1333649.981*u__dot[1]+32825.20357*`u__ double dot`[1]+127710.2711*u__dot[2]+3143.340237*`u__ double dot`[2]+12229.53067*u__dot[3]+301.0061407*`u__ double dot`[3]+1171.099388*u__dot[4]+28.82433650*`u__ double dot`[4]+112.1444325*u__dot[5]+2.760217359*`u__ double dot`[5]+10.73894656*u__dot[6]+.2643183086*`u__ double dot`[6]+1.028361111*u__dot[7]+0.2531111111e-1*`u__ double dot`[7]

 

1476876999.*u[0]+3143.340237*`u__ double dot`[3]+301.0061407*`u__ double dot`[4]+28.82433650*`u__ double dot`[5]+2.760217359*`u__ double dot`[6]+.2643183086*`u__ double dot`[7]+0.2531111111e-1*`u__ double dot`[8]+342786.3064*`u__ double dot`[1]+32825.20357*`u__ double dot`[2]+13927010.38*u__dot[1]+1333649.981*u__dot[2]+127710.2711*u__dot[3]+12229.53067*u__dot[4]+1171.099388*u__dot[5]+112.1444325*u__dot[6]+10.73894656*u__dot[7]+1.028361111*u__dot[8]+145436674.5*u__dot[0]+(1/18)*p[9]+7856982.494*p[1]+752384.3428*p[2]+72048.29583*p[3]+6899.342050*p[4]+660.6807306*p[5]+63.26676144*p[6]+6.058422650*p[7]+.5801543211*p[8]

 

0.1542269831e11*u[0]+32825.20357*`u__ double dot`[3]+3143.340237*`u__ double dot`[4]+301.0061407*`u__ double dot`[5]+28.82433650*`u__ double dot`[6]+2.760217359*`u__ double dot`[7]+.2643183086*`u__ double dot`[8]+0.2531111111e-1*`u__ double dot`[9]+3579641.223*`u__ double dot`[1]+342786.3064*`u__ double dot`[2]+145436674.5*u__dot[1]+13927010.38*u__dot[2]+1333649.981*u__dot[3]+127710.2711*u__dot[4]+12229.53067*u__dot[5]+1171.099388*u__dot[6]+112.1444325*u__dot[7]+10.73894656*u__dot[8]+1.028361111*u__dot[9]+1518762873.*u__dot[0]+.5801543211*p[9]+(1/18)*p[10]+82048722.17*p[1]+7856982.494*p[2]+752384.3428*p[3]+72048.29583*p[4]+6899.342050*p[5]+660.6807306*p[6]+63.26676144*p[7]+6.058422650*p[8]

 

0.1610558112e12*u[0]+342786.3064*`u__ double dot`[3]+32825.20357*`u__ double dot`[4]+3143.340237*`u__ double dot`[5]+301.0061407*`u__ double dot`[6]+28.82433650*`u__ double dot`[7]+2.760217359*`u__ double dot`[8]+.2643183086*`u__ double dot`[9]+0.2531111111e-1*`u__ double dot`[10]+37381397.82*`u__ double dot`[1]+3579641.223*`u__ double dot`[2]+1518762873.*u__dot[1]+145436674.5*u__dot[2]+13927010.38*u__dot[3]+1333649.981*u__dot[4]+127710.2711*u__dot[5]+12229.53067*u__dot[6]+1171.099388*u__dot[7]+112.1444325*u__dot[8]+10.73894656*u__dot[9]+1.028361111*u__dot[10]+0.1586010318e11*u__dot[0]+6.058422650*p[9]+.5801543211*p[10]+(1/18)*p[11]+856816572.8*p[1]+82048722.17*p[2]+7856982.494*p[3]+752384.3428*p[4]+72048.29583*p[5]+6899.342050*p[6]+660.6807306*p[7]+63.26676144*p[8]

(15)

for i from 0 to 10 do u__dot[i+1] := 2*u[i+1]/(.1)-u[i] end do;

1.111111111*p[1]+207.8555556*u[0]+20.56722222*u__dot[0]

 

1.111111111*p[2]+34.40000000*p[1]+6445.600778*u[0]+636.7611999*u__dot[0]+.5062222222*`u__ double dot`[1]

 

1.111111111*p[3]+34.40000000*p[2]+1065.601376*p[1]+199664.3295*u[0]+19724.81428*u__dot[0]+15.67264000*`u__ double dot`[1]+.5062222222*`u__ double dot`[2]

 

1.111111111*p[4]+34.40000000*p[3]+1065.601376*p[2]+33008.92305*p[1]+6184962.451*u[0]+611011.6704*u__dot[0]+485.4879870*`u__ double dot`[1]+15.67264000*`u__ double dot`[2]+.5062222222*`u__ double dot`[3]

 

191590358.6*u[0]+15.67264000*`u__ double dot`[3]+.5062222222*`u__ double dot`[4]+15038.86535*`u__ double dot`[1]+485.4879870*`u__ double dot`[2]+18927187.66*u__dot[0]+1022510.881*p[1]+33008.92305*p[2]+1065.601376*p[3]+34.40000000*p[4]+1.111111111*p[5]

 

5934856645.*u[0]+485.4879870*`u__ double dot`[3]+15.67264000*`u__ double dot`[4]+.5062222222*`u__ double dot`[5]+465855.9575*`u__ double dot`[1]+15038.86535*`u__ double dot`[2]+586303748.8*u__dot[0]+31674117.34*p[1]+1022510.881*p[2]+33008.92305*p[3]+1065.601376*p[4]+34.40000000*p[5]+1.111111111*p[6]

 

0.1838428805e12*u[0]+15038.86535*`u__ double dot`[3]+485.4879870*`u__ double dot`[4]+15.67264000*`u__ double dot`[5]+.5062222222*`u__ double dot`[6]+14430727.85*`u__ double dot`[1]+465855.9575*`u__ double dot`[2]+0.1816181527e11*u__dot[0]+981162868.1*p[1]+31674117.34*p[2]+1022510.881*p[3]+33008.92305*p[4]+1065.601376*p[5]+34.40000000*p[6]+1.111111111*p[7]

 

0.5694864547e13*u[0]+465855.9575*`u__ double dot`[3]+15038.86535*`u__ double dot`[4]+485.4879870*`u__ double dot`[5]+15.67264000*`u__ double dot`[6]+.5062222222*`u__ double dot`[7]+447017802.5*`u__ double dot`[1]+14430727.85*`u__ double dot`[2]+0.5625949595e12*u__dot[0]+0.3039328810e11*p[1]+981162868.1*p[2]+31674117.34*p[3]+1022510.881*p[4]+33008.92305*p[5]+1065.601376*p[6]+34.40000000*p[7]+1.111111111*p[8]

 

0.1764086927e15*u[0]+14430727.85*`u__ double dot`[3]+465855.9575*`u__ double dot`[4]+15038.86535*`u__ double dot`[5]+485.4879870*`u__ double dot`[6]+15.67264000*`u__ double dot`[7]+.5062222222*`u__ double dot`[8]+0.1384718205e11*`u__ double dot`[1]+447017802.5*`u__ double dot`[2]+0.1742739279e14*u__dot[0]+1.111111111*p[9]+0.9414868743e12*p[1]+0.3039328810e11*p[2]+981162868.3*p[3]+31674117.34*p[4]+1022510.881*p[5]+33008.92305*p[6]+1065.601376*p[7]+34.40000000*p[8]

 

0.5464577183e16*u[0]+447017802.5*`u__ double dot`[3]+14430727.85*`u__ double dot`[4]+465855.9575*`u__ double dot`[5]+15038.86535*`u__ double dot`[6]+485.4879870*`u__ double dot`[7]+15.67264000*`u__ double dot`[8]+.5062222222*`u__ double dot`[9]+0.4289414197e12*`u__ double dot`[1]+0.1384718205e11*`u__ double dot`[2]+0.5398448996e15*u__dot[0]+34.40000000*p[9]+1.111111111*p[10]+0.2916425269e14*p[1]+0.9414868743e12*p[2]+0.3039328810e11*p[3]+981162868.3*p[4]+31674117.34*p[5]+1022510.881*p[6]+33008.92305*p[7]+1065.601376*p[8]

 

0.1692751266e18*u[0]+0.1384718205e11*`u__ double dot`[3]+447017802.5*`u__ double dot`[4]+14430727.85*`u__ double dot`[5]+465855.9575*`u__ double dot`[6]+15038.86535*`u__ double dot`[7]+485.4879870*`u__ double dot`[8]+15.67264000*`u__ double dot`[9]+.5062222222*`u__ double dot`[10]+0.1328723353e14*`u__ double dot`[1]+0.4289414197e12*`u__ double dot`[2]+0.1672266869e17*u__dot[0]+1065.601376*p[9]+34.40000000*p[10]+1.111111111*p[11]+0.9034152876e15*p[1]+0.2916425269e14*p[2]+0.9414868743e12*p[3]+0.3039328810e11*p[4]+981162868.3*p[5]+31674117.34*p[6]+1022510.881*p[7]+33008.92305*p[8]

(16)

 

``

for i from 0 to 10 do `u__ double dot`[i+1] := 4*(u[i+1]-u[i])/.1^2-4*`u__ dot`[i+1]/(.1)-`u__ double dot`[i] end do;

22.22222222*p[1]+3777.111112*u[0]+411.3444444*u__dot[0]-40.00000000*`u__ dot`[1]

 

22.22222222*p[2]+869.6543212*p[1]+159407.8005*u[0]+16097.73632*u__dot[0]-364.9777776*`u__ dot`[1]-40.00000000*`u__ dot`[2]

 

22.22222222*p[3]+869.6543212*p[2]+35344.36401*p[1]+6472746.341*u[0]+654241.8501*u__dot[0]-15483.60256*`u__ dot`[1]-364.9777776*`u__ dot`[2]-40.00000000*`u__ dot`[3]

 

-628632.0873*`u__ dot`[1]-15483.60256*`u__ dot`[2]-364.9777776*`u__ dot`[3]-40.00000000*`u__ dot`[4]+262941585.9*u[0]+26576866.06*u__dot[0]+1435772.456*p[1]+35344.36401*p[2]+869.6543212*p[3]+22.22222222*p[4]

 

-25536885.15*`u__ dot`[1]-628632.0873*`u__ dot`[2]-15483.60256*`u__ dot`[3]-364.9777776*`u__ dot`[4]-40.00000000*`u__ dot`[5]+0.1068138112e11*u[0]+1079622608.*u__dot[0]+58324875.48*p[1]+1435772.456*p[2]+35344.36401*p[3]+869.6543212*p[4]+22.22222222*p[5]

 

-1037375646.*`u__ dot`[1]-25536885.15*`u__ dot`[2]-628632.0873*`u__ dot`[3]-15483.60256*`u__ dot`[4]-364.9777776*`u__ dot`[5]-40.00000000*`u__ dot`[6]+0.4339059231e12*u[0]+0.4385712279e11*u__dot[0]+2369310541.*p[1]+58324875.48*p[2]+1435772.456*p[3]+35344.36401*p[4]+869.6543212*p[5]+22.22222222*p[6]

 

-0.4214093965e11*`u__ dot`[1]-1037375646.*`u__ dot`[2]-25536885.15*`u__ dot`[3]-628632.0873*`u__ dot`[4]-15483.60256*`u__ dot`[5]-364.9777776*`u__ dot`[6]-40.00000000*`u__ dot`[7]+0.1762640503e14*u[0]+0.1781592203e13*u__dot[0]+0.9624765415e11*p[1]+2369310541.*p[2]+58324875.48*p[3]+1435772.456*p[4]+35344.36401*p[5]+869.6543212*p[6]+22.22222222*p[7]

 

-0.1711876309e13*`u__ dot`[1]-0.4214093965e11*`u__ dot`[2]-1037375645.*`u__ dot`[3]-25536885.15*`u__ dot`[4]-628632.0873*`u__ dot`[5]-15483.60256*`u__ dot`[6]-364.9777776*`u__ dot`[7]-40.00000000*`u__ dot`[8]+0.7160311434e15*u[0]+0.7237298245e14*u__dot[0]+0.3909834009e13*p[1]+0.9624765415e11*p[2]+2369310541.*p[3]+58324875.48*p[4]+1435772.456*p[5]+35344.36401*p[6]+869.6543212*p[7]+22.22222222*p[8]

 

-0.6954093867e14*`u__ dot`[1]-0.1711876309e13*`u__ dot`[2]-0.4214093962e11*`u__ dot`[3]-1037375646.*`u__ dot`[4]-25536885.15*`u__ dot`[5]-628632.0873*`u__ dot`[6]-15483.60256*`u__ dot`[7]-364.9777776*`u__ dot`[8]-40.00000000*`u__ dot`[9]+0.2908707689e17*u[0]+0.2939981766e16*u__dot[0]+22.22222222*p[9]+0.1588277878e15*p[1]+0.3909834009e13*p[2]+0.9624765415e11*p[3]+2369310540.*p[4]+58324875.44*p[5]+1435772.456*p[6]+35344.36401*p[7]+869.6543212*p[8]

 

-0.2824936664e16*`u__ dot`[1]-0.6954093867e14*`u__ dot`[2]-0.1711876309e13*`u__ dot`[3]-0.4214093965e11*`u__ dot`[4]-1037375646.*`u__ dot`[5]-25536885.15*`u__ dot`[6]-628632.0873*`u__ dot`[7]-15483.60256*`u__ dot`[8]-364.9777776*`u__ dot`[9]-40.00000000*`u__ dot`[10]+0.1181593915e19*u[0]+0.1194298271e18*u__dot[0]+869.6543212*p[9]+22.22222222*p[10]+0.6452004379e16*p[1]+0.1588277878e15*p[2]+0.3909834009e13*p[3]+0.9624765415e11*p[4]+2369310540.*p[5]+58324875.44*p[6]+1435772.456*p[7]+35344.36401*p[8]

 

-0.1147563911e18*`u__ dot`[1]-0.2824936664e16*`u__ dot`[2]-0.6954093867e14*`u__ dot`[3]-0.1711876310e13*`u__ dot`[4]-0.4214093965e11*`u__ dot`[5]-1037375646.*`u__ dot`[6]-25536885.15*`u__ dot`[7]-628632.0873*`u__ dot`[8]-15483.60256*`u__ dot`[9]-364.9777776*`u__ dot`[10]-40.00000000*`u__ dot`[11]+0.4799946674e20*u[0]+0.4851555123e19*u__dot[0]+35344.36401*p[9]+869.6543212*p[10]+22.22222222*p[11]+0.2620974648e18*p[1]+0.6452004379e16*p[2]+0.1588277877e15*p[3]+0.3909834009e13*p[4]+0.9624765415e11*p[5]+2369310540.*p[6]+58324875.44*p[7]+1435772.456*p[8]

(17)

slon := fsolve({0, p[11]+187.9700000*u[10]+18.51050000*u__dot[10]+.4556*`u__ double dot`[10], 0.1692751266e18*u[0]+0.1384718205e11*`u__ double dot`[3]+447017802.5*`u__ double dot`[4]+14430727.85*`u__ double dot`[5]+465855.9575*`u__ double dot`[6]+15038.86535*`u__ double dot`[7]+485.4879870*`u__ double dot`[8]+15.67264000*`u__ double dot`[9]+.5062222222*`u__ double dot`[10]+0.1328723353e14*`u__ double dot`[1]+0.4289414197e12*`u__ double dot`[2]+0.1672266869e17*u__dot[0]+1065.601376*p[9]+34.40000000*p[10]+1.111111111*p[11]+0.9034152876e15*p[1]+0.2916425269e14*p[2]+0.9414868743e12*p[3]+0.3039328810e11*p[4]+981162868.3*p[5]+31674117.34*p[6]+1022510.881*p[7]+33008.92305*p[8], -0.1147563911e18*`u__ dot`[1]-0.2824936664e16*`u__ dot`[2]-0.6954093867e14*`u__ dot`[3]-0.1711876310e13*`u__ dot`[4]-0.4214093965e11*`u__ dot`[5]-1037375646.*`u__ dot`[6]-25536885.15*`u__ dot`[7]-628632.0873*`u__ dot`[8]-15483.60256*`u__ dot`[9]-364.9777776*`u__ dot`[10]-40.00000000*`u__ dot`[11]+0.4799946674e20*u[0]+0.4851555123e19*u__dot[0]+35344.36401*p[9]+869.6543212*p[10]+22.22222222*p[11]+0.2620974648e18*p[1]+0.6452004379e16*p[2]+0.1588277877e15*p[3]+0.3909834009e13*p[4]+0.9624765415e11*p[5]+2369310540.*p[6]+58324875.44*p[7]+1435772.456*p[8], 0.1610558112e12*u[0]+342786.3064*`u__ double dot`[3]+32825.20357*`u__ double dot`[4]+3143.340237*`u__ double dot`[5]+301.0061407*`u__ double dot`[6]+28.82433650*`u__ double dot`[7]+2.760217359*`u__ double dot`[8]+.2643183086*`u__ double dot`[9]+0.2531111111e-1*`u__ double dot`[10]+37381397.82*`u__ double dot`[1]+3579641.223*`u__ double dot`[2]+1518762873.*u__dot[1]+145436674.5*u__dot[2]+13927010.38*u__dot[3]+1333649.981*u__dot[4]+127710.2711*u__dot[5]+12229.53067*u__dot[6]+1171.099388*u__dot[7]+112.1444325*u__dot[8]+10.73894656*u__dot[9]+1.028361111*u__dot[10]+0.1586010318e11*u__dot[0]+6.058422650*p[9]+.5801543211*p[10]+(1/18)*p[11]+856816572.8*p[1]+82048722.17*p[2]+7856982.494*p[3]+752384.3428*p[4]+72048.29583*p[5]+6899.342050*p[6]+660.6807306*p[7]+63.26676144*p[8]});

{p[1] = 0.2999999998e-1*`u__ dot`[2]+0.9374999995e-13*`u__ dot`[3]-0.1499999999e-4*`u__ dot`[4]+0.1312499999e-15*`u__ dot`[5]+0.4999999999e-8*`u__ dot`[6]-0.6999999938e-9*`u__ dot`[7]+0.9899999882e-9*`u__ dot`[8]-0.1827374978e-8*`u__ dot`[9]+0.3743087460e-8*`u__ dot`[10]+0.7499999993e-16*`u__ dot`[11]+0.4597499946e-10*p[9]-0.1051124990e-9*p[10]-0.2053929019e-9*p[11]-0.2499999999e-1*p[2]-0.2499999998e-3*p[3]-0.2499999999e-4*p[4]-0.2499999999e-6*p[5]-0.5000000000e-8*p[6]-0.9999999992e-9*p[7]-0.3749999970e-10*p[8], p[2] = p[2], p[3] = p[3], p[4] = p[4], p[5] = p[5], p[6] = p[6], p[7] = p[7], p[8] = p[8], p[9] = p[9], p[10] = p[10], p[11] = p[11], u[0] = 0.281894999e-4*`u__ dot`[2]+0.7830419785e-6*`u__ dot`[3]+0.318079236e-8*`u__ dot`[4]+0.9022187373e-9*`u__ dot`[5]-0.1628286346e-10*`u__ dot`[6]+0.1100239238e-9*`u__ dot`[7]-0.2048575058e-9*`u__ dot`[8]+0.3826992048e-9*`u__ dot`[9]-0.7123420880e-9*`u__ dot`[10]+0.1380032829e-8*`u__ dot`[11]-0.9758339382e-11*p[9]+0.1808837852e-10*p[10]-0.3663799046e-10*p[11]-0.989828263e-4*p[2]-0.884545513e-6*p[3]-0.8052117590e-7*p[4]-0.1650827489e-8*p[5]-0.2205906206e-10*p[6]-0.3199234662e-11*p[7]+0.5195642082e-11*p[8], `u__ dot`[1] = -0.4999999998e-1*`u__ dot`[2]+0.2499999997e-3*`u__ dot`[3]+0.4999999998e-5*`u__ dot`[4]+0.7499999995e-6*`u__ dot`[5]-0.1499999999e-7*`u__ dot`[6]+0.7949999997e-7*`u__ dot`[7]-0.1479199999e-6*`u__ dot`[8]+0.2763269999e-6*`u__ dot`[9]-0.5141608198e-6*`u__ dot`[10]+0.9999999996e-6*`u__ dot`[11]-0.7045749997e-8*p[9]+0.1305056749e-7*p[10]-0.2665930549e-7*p[11]+0.1499999999e-2*p[3]-0.2499999999e-4*p[4]-0.9999999996e-9*p[7]+0.3769999998e-8*p[8], `u__ dot`[2] = `u__ dot`[2], `u__ dot`[3] = `u__ dot`[3], `u__ dot`[4] = `u__ dot`[4], `u__ dot`[5] = `u__ dot`[5], `u__ dot`[6] = `u__ dot`[6], `u__ dot`[7] = `u__ dot`[7], `u__ dot`[8] = `u__ dot`[8], `u__ dot`[9] = `u__ dot`[9], `u__ dot`[10] = `u__ dot`[10], `u__ dot`[11] = `u__ dot`[11], u__dot[0] = -0.2499999998e-2*`u__ dot`[2]+0.1249999999e-4*`u__ dot`[3]+0.1249999999e-5*`u__ dot`[4]+0.1749999998e-7*`u__ dot`[5]-0.2499999998e-9*`u__ dot`[6]+0.8349999992e-9*`u__ dot`[7]-0.1525399998e-8*`u__ dot`[8]+0.2848549998e-8*`u__ dot`[9]-0.5316285694e-8*`u__ dot`[10]+0.9999999990e-8*`u__ dot`[11]-0.7260249992e-10*p[9]+0.1354106748e-9*p[10]-0.2570080548e-9*p[11]+0.9999999994e-3*p[2]+0.2499999997e-4*p[3]+0.7499999994e-6*p[4]+0.9999999994e-8*p[5]+0.4999999995e-10*p[7]+0.3949999996e-10*p[8]}

(18)

``


 

Download hw_4_structural.mw

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