Question: plotting problem

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

Please Wait...