samen

25 Reputation

5 Badges

3 years, 68 days

MaplePrimes Activity


These are replies submitted by samen


 

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-3*k

(1)

for i from 0 while i <= N do u[i, 0] := h*i+1 end do:

for j from 0 while j <= N+1 do u[0, j] := .1; u[N+1, j] := .5 end do:

printlevel := 2:

NULL

NULL

sys := ([seq])(seq(eq[i, j], j = 0 .. N), i = 1 .. N):

nops(sys);

vars:=indets(sys) minus {k}:

nn := Matrix(N+1, N+1,(i, j)-> u[i-1, j-1]):

##

p:=proc(kk) local u_res,A;

  u_res:=solve(eval(sys,k=kk),vars);

  A:=eval(nn,u_res);

  plots:-matrixplot(A)

end proc;

## Testing p for k=0.001:

p(0.001);

## Animating the plot for k=0.0001..0.001:

plots:-animate(p,[k],k=0.0001..0.001);

 

90

 

proc (kk) local u_res, A; u_res := solve(eval(sys, k = kk), vars); A := eval(nn, u_res); plots:-matrixplot(A) end proc

 

 

 

 

``


 

Download Crank_scheme_4.mw

@Preben Alsholm 

Thank you...
 

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-3*k

(1)

for i from 0 while i <= N do u[i, 0] := h*i+1 end do:

for j from 0 while j <= N+1 do u[0, j] := .1; u[N+1, j] := .5 end do:

printlevel := 2:

``

``

sys := ([seq])(seq(eq[i, j], j = 0 .. N), i = 1 .. N):

nops(sys);

vars:=indets(sys) minus {k}:

nn := Matrix(N+1, N+1,(i, j)-> u[i-1, j-1]):

##

p:=proc(kk) local u_res,A;

  u_res:=solve(eval(sys,k=kk),vars);

  A:=eval(nn,u_res);

  plots:-matrixplot(A)

end proc;

## Testing p for k=10000:

p(10000);

## Animating the plot for k=0..10000:

plots:-animate(p,[k],k=0..10000);

 

90

 

proc (kk) local u_res, A; u_res := solve(eval(sys, k = kk), vars); A := eval(nn, u_res); plots:-matrixplot(A) end proc

 

 

 

 

NULL


 

Download Crank_scheme.mw

@Preben Alsholm 

i want to find the change in results with chnage in 'k'. system consider lemda=0 bus lemda has some value there for no change in values...any command for this ???

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-6

(1)

for i from 0 while i <= N+1 do u[i, 0] := h*i+1 end do;

1

 

11/10

 

6/5

 

13/10

 

7/5

 

3/2

 

8/5

 

17/10

 

9/5

 

19/10

 

2

(2)

for j from 0 while j <= N+1 do u[0, j] := .1; u[N+1, j] := .5 end do;

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

 

.1

 

.5

(3)

printlevel := 2:

sys := ([seq])(seq(eq[i, j], j = 0 .. N), i = 1 .. N):

fsolve(sys);

{u[1, 1] = 1.1000, u[1, 2] = 1.1000, u[1, 3] = 1.1000, u[1, 4] = 1.1000, u[1, 5] = 1.1000, u[1, 6] = 1.1000, u[1, 7] = 1.1000, u[1, 8] = 1.1000, u[1, 9] = 1.1000, u[1, 10] = 1.1000, u[2, 1] = 1.2000, u[2, 2] = 1.2000, u[2, 3] = 1.2000, u[2, 4] = 1.2000, u[2, 5] = 1.2000, u[2, 6] = 1.2000, u[2, 7] = 1.2000, u[2, 8] = 1.2000, u[2, 9] = 1.2000, u[2, 10] = 1.2000, u[3, 1] = 1.3000, u[3, 2] = 1.3000, u[3, 3] = 1.3000, u[3, 4] = 1.3000, u[3, 5] = 1.3000, u[3, 6] = 1.3000, u[3, 7] = 1.3000, u[3, 8] = 1.3000, u[3, 9] = 1.3000, u[3, 10] = 1.3000, u[4, 1] = 1.4000, u[4, 2] = 1.4000, u[4, 3] = 1.4000, u[4, 4] = 1.4000, u[4, 5] = 1.4000, u[4, 6] = 1.4000, u[4, 7] = 1.4000, u[4, 8] = 1.4000, u[4, 9] = 1.4000, u[4, 10] = 1.4000, u[5, 1] = 1.5000, u[5, 2] = 1.5000, u[5, 3] = 1.5000, u[5, 4] = 1.5000, u[5, 5] = 1.5000, u[5, 6] = 1.5000, u[5, 7] = 1.5000, u[5, 8] = 1.5000, u[5, 9] = 1.5000, u[5, 10] = 1.5000, u[6, 1] = 1.6000, u[6, 2] = 1.6000, u[6, 3] = 1.6000, u[6, 4] = 1.6000, u[6, 5] = 1.6000, u[6, 6] = 1.6000, u[6, 7] = 1.6000, u[6, 8] = 1.6000, u[6, 9] = 1.6000, u[6, 10] = 1.6000, u[7, 1] = 1.7000, u[7, 2] = 1.7000, u[7, 3] = 1.7000, u[7, 4] = 1.7000, u[7, 5] = 1.7000, u[7, 6] = 1.7000, u[7, 7] = 1.7000, u[7, 8] = 1.7000, u[7, 9] = 1.7000, u[7, 10] = 1.7000, u[8, 1] = 1.8000, u[8, 2] = 1.8000, u[8, 3] = 1.8000, u[8, 4] = 1.8000, u[8, 5] = 1.8000, u[8, 6] = 1.8000, u[8, 7] = 1.8000, u[8, 8] = 1.8000, u[8, 9] = 1.8000, u[8, 10] = 1.8000, u[9, 1] = 1.9000, u[9, 2] = 1.9000, u[9, 3] = 1.9000, u[9, 4] = 1.9000, u[9, 5] = 1.9000, u[9, 6] = 1.9000, u[9, 7] = 1.9000, u[9, 8] = 1.9000, u[9, 9] = 1.9000, u[9, 10] = 1.9000}

(4)

nn := Matrix(N+1, N+1, proc (i, j) options operator, arrow; u[i-1, j-1] end proc); evalm(nn)

array( 1 .. 10, 1 .. 10, [( 1, 2 ) = (.1), ( 9, 8 ) = (u[8, 7]), ( 4, 2 ) = (u[3, 1]), ( 1, 10 ) = (.1), ( 2, 7 ) = (u[1, 6]), ( 5, 2 ) = (u[4, 1]), ( 1, 3 ) = (.1), ( 6, 3 ) = (u[5, 2]), ( 9, 5 ) = (u[8, 4]), ( 2, 9 ) = (u[1, 8]), ( 5, 5 ) = (u[4, 4]), ( 7, 10 ) = (u[6, 9]), ( 1, 5 ) = (.1), ( 2, 2 ) = (u[1, 1]), ( 1, 1 ) = (.1), ( 6, 6 ) = (u[5, 5]), ( 4, 8 ) = (u[3, 7]), ( 1, 6 ) = (.1), ( 2, 3 ) = (u[1, 2]), ( 8, 8 ) = (u[7, 7]), ( 10, 4 ) = (u[9, 3]), ( 9, 1 ) = (9/5), ( 7, 7 ) = (u[6, 6]), ( 1, 7 ) = (.1), ( 2, 4 ) = (u[1, 3]), ( 10, 5 ) = (u[9, 4]), ( 3, 4 ) = (u[2, 3]), ( 6, 9 ) = (u[5, 8]), ( 4, 5 ) = (u[3, 4]), ( 1, 8 ) = (.1), ( 2, 5 ) = (u[1, 4]), ( 4, 1 ) = (13/10), ( 1, 9 ) = (.1), ( 7, 4 ) = (u[6, 3]), ( 2, 6 ) = (u[1, 5]), ( 5, 1 ) = (7/5), ( 3, 7 ) = (u[2, 6]), ( 4, 6 ) = (u[3, 5]), ( 8, 10 ) = (u[7, 9]), ( 2, 8 ) = (u[1, 7]), ( 8, 1 ) = (17/10), ( 3, 2 ) = (u[2, 1]), ( 6, 7 ) = (u[5, 6]), ( 1, 4 ) = (.1), ( 5, 8 ) = (u[4, 7]), ( 8, 5 ) = (u[7, 4]), ( 10, 2 ) = (u[9, 1]), ( 7, 6 ) = (u[6, 5]), ( 6, 2 ) = (u[5, 1]), ( 9, 2 ) = (u[8, 1]), ( 5, 9 ) = (u[4, 8]), ( 8, 4 ) = (u[7, 3]), ( 3, 8 ) = (u[2, 7]), ( 9, 7 ) = (u[8, 6]), ( 5, 3 ) = (u[4, 2]), ( 10, 9 ) = (u[9, 8]), ( 10, 7 ) = (u[9, 6]), ( 3, 3 ) = (u[2, 2]), ( 6, 8 ) = (u[5, 7]), ( 4, 10 ) = (u[3, 9]), ( 7, 1 ) = (8/5), ( 5, 7 ) = (u[4, 6]), ( 8, 6 ) = (u[7, 5]), ( 7, 5 ) = (u[6, 4]), ( 3, 6 ) = (u[2, 5]), ( 4, 7 ) = (u[3, 6]), ( 8, 9 ) = (u[7, 8]), ( 7, 2 ) = (u[6, 1]), ( 2, 1 ) = (11/10), ( 3, 9 ) = (u[2, 8]), ( 9, 6 ) = (u[8, 5]), ( 5, 4 ) = (u[4, 3]), ( 6, 1 ) = (3/2), ( 9, 3 ) = (u[8, 2]), ( 5, 10 ) = (u[4, 9]), ( 8, 3 ) = (u[7, 2]), ( 6, 4 ) = (u[5, 3]), ( 9, 4 ) = (u[8, 3]), ( 5, 6 ) = (u[4, 5]), ( 8, 2 ) = (u[7, 1]), ( 7, 9 ) = (u[6, 8]), ( 3, 10 ) = (u[2, 9]), ( 9, 9 ) = (u[8, 8]), ( 10, 6 ) = (u[9, 5]), ( 4, 3 ) = (u[3, 2]), ( 2, 10 ) = (u[1, 9]), ( 10, 10 ) = (u[9, 9]), ( 3, 5 ) = (u[2, 4]), ( 6, 10 ) = (u[5, 9]), ( 9, 10 ) = (u[8, 9]), ( 4, 4 ) = (u[3, 3]), ( 3, 1 ) = (6/5), ( 7, 3 ) = (u[6, 2]), ( 6, 5 ) = (u[5, 4]), ( 4, 9 ) = (u[3, 8]), ( 10, 1 ) = (19/10), ( 8, 7 ) = (u[7, 6]), ( 10, 8 ) = (u[9, 7]), ( 7, 8 ) = (u[6, 7]), ( 10, 3 ) = (u[9, 2])  ] )

(5)

 

Plot(

``

``


 

Download plot_crnk.mw

@tomleslie 


 

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-6

(1)

for i from 0 while i <= N do u[i, 0] := h*i+1 end do:

for j from 0 while j < N+1 do u[0, j] := .1; u[N+1, j] := .5 end do:

printlevel := 2:

sys := ([seq])(seq(eq[i, j], j = 0 .. 9), i = 1 .. 9):

fsolve(sys);

{eq[1, 0] = 0., eq[1, 1] = 0., eq[1, 2] = 0., eq[1, 3] = 0., eq[1, 4] = 0., eq[1, 5] = 0., eq[1, 6] = 0., eq[1, 7] = 0., eq[1, 8] = 0., eq[1, 9] = 0., eq[2, 0] = 0., eq[2, 1] = 0., eq[2, 2] = 0., eq[2, 3] = 0., eq[2, 4] = 0., eq[2, 5] = 0., eq[2, 6] = 0., eq[2, 7] = 0., eq[2, 8] = 0., eq[2, 9] = 0., eq[3, 0] = 0., eq[3, 1] = 0., eq[3, 2] = 0., eq[3, 3] = 0., eq[3, 4] = 0., eq[3, 5] = 0., eq[3, 6] = 0., eq[3, 7] = 0., eq[3, 8] = 0., eq[3, 9] = 0., eq[4, 0] = 0., eq[4, 1] = 0., eq[4, 2] = 0., eq[4, 3] = 0., eq[4, 4] = 0., eq[4, 5] = 0., eq[4, 6] = 0., eq[4, 7] = 0., eq[4, 8] = 0., eq[4, 9] = 0., eq[5, 0] = 0., eq[5, 1] = 0., eq[5, 2] = 0., eq[5, 3] = 0., eq[5, 4] = 0., eq[5, 5] = 0., eq[5, 6] = 0., eq[5, 7] = 0., eq[5, 8] = 0., eq[5, 9] = 0., eq[6, 0] = 0., eq[6, 1] = 0., eq[6, 2] = 0., eq[6, 3] = 0., eq[6, 4] = 0., eq[6, 5] = 0., eq[6, 6] = 0., eq[6, 7] = 0., eq[6, 8] = 0., eq[6, 9] = 0., eq[7, 0] = 0., eq[7, 1] = 0., eq[7, 2] = 0., eq[7, 3] = 0., eq[7, 4] = 0., eq[7, 5] = 0., eq[7, 6] = 0., eq[7, 7] = 0., eq[7, 8] = 0., eq[7, 9] = 0., eq[8, 0] = 0., eq[8, 1] = 0., eq[8, 2] = 0., eq[8, 3] = 0., eq[8, 4] = 0., eq[8, 5] = 0., eq[8, 6] = 0., eq[8, 7] = 0., eq[8, 8] = 0., eq[8, 9] = 0., eq[9, 0] = 0., eq[9, 1] = 0., eq[9, 2] = 0., eq[9, 3] = 0., eq[9, 4] = 0., eq[9, 5] = 0., eq[9, 6] = 0., eq[9, 7] = 0., eq[9, 8] = 0., eq[9, 9] = 0.}

(2)

``


 

Download Crank_scheme_2.mw

@mehdi jafari 


 

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-6

(1)

for i from 0 while i <= N do u[i, 0] := h*i+1 end do:

for j from 0 while j < N+1 do u[0, j] := .1; u[N+1, j] := .5 end do:

printlevel := 2:

2.2000 = -0.11976e-7+2.*u[1, 1]-0.11976e-6*u[2, 1]

 

0.11976e-7+2.*u[1, 1]+0.11976e-6*u[2, 1] = -0.11976e-7+2.*u[1, 2]-0.11976e-6*u[2, 2]

 

0.11976e-7+2.*u[1, 2]+0.11976e-6*u[2, 2] = -0.11976e-7+2.*u[1, 3]-0.11976e-6*u[2, 3]

 

0.11976e-7+2.*u[1, 3]+0.11976e-6*u[2, 3] = -0.11976e-7+2.*u[1, 4]-0.11976e-6*u[2, 4]

 

0.11976e-7+2.*u[1, 4]+0.11976e-6*u[2, 4] = -0.11976e-7+2.*u[1, 5]-0.11976e-6*u[2, 5]

 

0.11976e-7+2.*u[1, 5]+0.11976e-6*u[2, 5] = -0.11976e-7+2.*u[1, 6]-0.11976e-6*u[2, 6]

 

0.11976e-7+2.*u[1, 6]+0.11976e-6*u[2, 6] = -0.11976e-7+2.*u[1, 7]-0.11976e-6*u[2, 7]

 

0.11976e-7+2.*u[1, 7]+0.11976e-6*u[2, 7] = -0.11976e-7+2.*u[1, 8]-0.11976e-6*u[2, 8]

 

0.11976e-7+2.*u[1, 8]+0.11976e-6*u[2, 8] = -0.11976e-7+2.*u[1, 9]-0.11976e-6*u[2, 9]

 

0.11976e-7+2.*u[1, 9]+0.11976e-6*u[2, 9] = -0.11976e-6*u[0, 10]+2.*u[1, 10]-0.11976e-6*u[2, 10]

 

2.4000 = -0.11976e-6*u[1, 1]+2.*u[2, 1]-0.11976e-6*u[3, 1]

 

0.11976e-6*u[1, 1]+2.*u[2, 1]+0.11976e-6*u[3, 1] = -0.11976e-6*u[1, 2]+2.*u[2, 2]-0.11976e-6*u[3, 2]

 

0.11976e-6*u[1, 2]+2.*u[2, 2]+0.11976e-6*u[3, 2] = -0.11976e-6*u[1, 3]+2.*u[2, 3]-0.11976e-6*u[3, 3]

 

0.11976e-6*u[1, 3]+2.*u[2, 3]+0.11976e-6*u[3, 3] = -0.11976e-6*u[1, 4]+2.*u[2, 4]-0.11976e-6*u[3, 4]

 

0.11976e-6*u[1, 4]+2.*u[2, 4]+0.11976e-6*u[3, 4] = -0.11976e-6*u[1, 5]+2.*u[2, 5]-0.11976e-6*u[3, 5]

 

0.11976e-6*u[1, 5]+2.*u[2, 5]+0.11976e-6*u[3, 5] = -0.11976e-6*u[1, 6]+2.*u[2, 6]-0.11976e-6*u[3, 6]

 

0.11976e-6*u[1, 6]+2.*u[2, 6]+0.11976e-6*u[3, 6] = -0.11976e-6*u[1, 7]+2.*u[2, 7]-0.11976e-6*u[3, 7]

 

0.11976e-6*u[1, 7]+2.*u[2, 7]+0.11976e-6*u[3, 7] = -0.11976e-6*u[1, 8]+2.*u[2, 8]-0.11976e-6*u[3, 8]

 

0.11976e-6*u[1, 8]+2.*u[2, 8]+0.11976e-6*u[3, 8] = -0.11976e-6*u[1, 9]+2.*u[2, 9]-0.11976e-6*u[3, 9]

 

0.11976e-6*u[1, 9]+2.*u[2, 9]+0.11976e-6*u[3, 9] = -0.11976e-6*u[1, 10]+2.*u[2, 10]-0.11976e-6*u[3, 10]

 

2.6000 = -0.11976e-6*u[2, 1]+2.*u[3, 1]-0.11976e-6*u[4, 1]

 

0.11976e-6*u[2, 1]+2.*u[3, 1]+0.11976e-6*u[4, 1] = -0.11976e-6*u[2, 2]+2.*u[3, 2]-0.11976e-6*u[4, 2]

 

0.11976e-6*u[2, 2]+2.*u[3, 2]+0.11976e-6*u[4, 2] = -0.11976e-6*u[2, 3]+2.*u[3, 3]-0.11976e-6*u[4, 3]

 

0.11976e-6*u[2, 3]+2.*u[3, 3]+0.11976e-6*u[4, 3] = -0.11976e-6*u[2, 4]+2.*u[3, 4]-0.11976e-6*u[4, 4]

 

0.11976e-6*u[2, 4]+2.*u[3, 4]+0.11976e-6*u[4, 4] = -0.11976e-6*u[2, 5]+2.*u[3, 5]-0.11976e-6*u[4, 5]

 

0.11976e-6*u[2, 5]+2.*u[3, 5]+0.11976e-6*u[4, 5] = -0.11976e-6*u[2, 6]+2.*u[3, 6]-0.11976e-6*u[4, 6]

 

0.11976e-6*u[2, 6]+2.*u[3, 6]+0.11976e-6*u[4, 6] = -0.11976e-6*u[2, 7]+2.*u[3, 7]-0.11976e-6*u[4, 7]

 

0.11976e-6*u[2, 7]+2.*u[3, 7]+0.11976e-6*u[4, 7] = -0.11976e-6*u[2, 8]+2.*u[3, 8]-0.11976e-6*u[4, 8]

 

0.11976e-6*u[2, 8]+2.*u[3, 8]+0.11976e-6*u[4, 8] = -0.11976e-6*u[2, 9]+2.*u[3, 9]-0.11976e-6*u[4, 9]

 

0.11976e-6*u[2, 9]+2.*u[3, 9]+0.11976e-6*u[4, 9] = -0.11976e-6*u[2, 10]+2.*u[3, 10]-0.11976e-6*u[4, 10]

 

2.8000 = -0.11976e-6*u[3, 1]+2.*u[4, 1]-0.11976e-6*u[5, 1]

 

0.11976e-6*u[3, 1]+2.*u[4, 1]+0.11976e-6*u[5, 1] = -0.11976e-6*u[3, 2]+2.*u[4, 2]-0.11976e-6*u[5, 2]

 

0.11976e-6*u[3, 2]+2.*u[4, 2]+0.11976e-6*u[5, 2] = -0.11976e-6*u[3, 3]+2.*u[4, 3]-0.11976e-6*u[5, 3]

 

0.11976e-6*u[3, 3]+2.*u[4, 3]+0.11976e-6*u[5, 3] = -0.11976e-6*u[3, 4]+2.*u[4, 4]-0.11976e-6*u[5, 4]

 

0.11976e-6*u[3, 4]+2.*u[4, 4]+0.11976e-6*u[5, 4] = -0.11976e-6*u[3, 5]+2.*u[4, 5]-0.11976e-6*u[5, 5]

 

0.11976e-6*u[3, 5]+2.*u[4, 5]+0.11976e-6*u[5, 5] = -0.11976e-6*u[3, 6]+2.*u[4, 6]-0.11976e-6*u[5, 6]

 

0.11976e-6*u[3, 6]+2.*u[4, 6]+0.11976e-6*u[5, 6] = -0.11976e-6*u[3, 7]+2.*u[4, 7]-0.11976e-6*u[5, 7]

 

0.11976e-6*u[3, 7]+2.*u[4, 7]+0.11976e-6*u[5, 7] = -0.11976e-6*u[3, 8]+2.*u[4, 8]-0.11976e-6*u[5, 8]

 

0.11976e-6*u[3, 8]+2.*u[4, 8]+0.11976e-6*u[5, 8] = -0.11976e-6*u[3, 9]+2.*u[4, 9]-0.11976e-6*u[5, 9]

 

0.11976e-6*u[3, 9]+2.*u[4, 9]+0.11976e-6*u[5, 9] = -0.11976e-6*u[3, 10]+2.*u[4, 10]-0.11976e-6*u[5, 10]

 

3.0000 = -0.11976e-6*u[4, 1]+2.*u[5, 1]-0.11976e-6*u[6, 1]

 

0.11976e-6*u[4, 1]+2.*u[5, 1]+0.11976e-6*u[6, 1] = -0.11976e-6*u[4, 2]+2.*u[5, 2]-0.11976e-6*u[6, 2]

 

0.11976e-6*u[4, 2]+2.*u[5, 2]+0.11976e-6*u[6, 2] = -0.11976e-6*u[4, 3]+2.*u[5, 3]-0.11976e-6*u[6, 3]

 

0.11976e-6*u[4, 3]+2.*u[5, 3]+0.11976e-6*u[6, 3] = -0.11976e-6*u[4, 4]+2.*u[5, 4]-0.11976e-6*u[6, 4]

 

0.11976e-6*u[4, 4]+2.*u[5, 4]+0.11976e-6*u[6, 4] = -0.11976e-6*u[4, 5]+2.*u[5, 5]-0.11976e-6*u[6, 5]

 

0.11976e-6*u[4, 5]+2.*u[5, 5]+0.11976e-6*u[6, 5] = -0.11976e-6*u[4, 6]+2.*u[5, 6]-0.11976e-6*u[6, 6]

 

0.11976e-6*u[4, 6]+2.*u[5, 6]+0.11976e-6*u[6, 6] = -0.11976e-6*u[4, 7]+2.*u[5, 7]-0.11976e-6*u[6, 7]

 

0.11976e-6*u[4, 7]+2.*u[5, 7]+0.11976e-6*u[6, 7] = -0.11976e-6*u[4, 8]+2.*u[5, 8]-0.11976e-6*u[6, 8]

 

0.11976e-6*u[4, 8]+2.*u[5, 8]+0.11976e-6*u[6, 8] = -0.11976e-6*u[4, 9]+2.*u[5, 9]-0.11976e-6*u[6, 9]

 

0.11976e-6*u[4, 9]+2.*u[5, 9]+0.11976e-6*u[6, 9] = -0.11976e-6*u[4, 10]+2.*u[5, 10]-0.11976e-6*u[6, 10]

 

3.2000 = -0.11976e-6*u[5, 1]+2.*u[6, 1]-0.11976e-6*u[7, 1]

 

0.11976e-6*u[5, 1]+2.*u[6, 1]+0.11976e-6*u[7, 1] = -0.11976e-6*u[5, 2]+2.*u[6, 2]-0.11976e-6*u[7, 2]

 

0.11976e-6*u[5, 2]+2.*u[6, 2]+0.11976e-6*u[7, 2] = -0.11976e-6*u[5, 3]+2.*u[6, 3]-0.11976e-6*u[7, 3]

 

0.11976e-6*u[5, 3]+2.*u[6, 3]+0.11976e-6*u[7, 3] = -0.11976e-6*u[5, 4]+2.*u[6, 4]-0.11976e-6*u[7, 4]

 

0.11976e-6*u[5, 4]+2.*u[6, 4]+0.11976e-6*u[7, 4] = -0.11976e-6*u[5, 5]+2.*u[6, 5]-0.11976e-6*u[7, 5]

 

0.11976e-6*u[5, 5]+2.*u[6, 5]+0.11976e-6*u[7, 5] = -0.11976e-6*u[5, 6]+2.*u[6, 6]-0.11976e-6*u[7, 6]

 

0.11976e-6*u[5, 6]+2.*u[6, 6]+0.11976e-6*u[7, 6] = -0.11976e-6*u[5, 7]+2.*u[6, 7]-0.11976e-6*u[7, 7]

 

0.11976e-6*u[5, 7]+2.*u[6, 7]+0.11976e-6*u[7, 7] = -0.11976e-6*u[5, 8]+2.*u[6, 8]-0.11976e-6*u[7, 8]

 

0.11976e-6*u[5, 8]+2.*u[6, 8]+0.11976e-6*u[7, 8] = -0.11976e-6*u[5, 9]+2.*u[6, 9]-0.11976e-6*u[7, 9]

 

0.11976e-6*u[5, 9]+2.*u[6, 9]+0.11976e-6*u[7, 9] = -0.11976e-6*u[5, 10]+2.*u[6, 10]-0.11976e-6*u[7, 10]

 

3.4000 = -0.11976e-6*u[6, 1]+2.*u[7, 1]-0.11976e-6*u[8, 1]

 

0.11976e-6*u[6, 1]+2.*u[7, 1]+0.11976e-6*u[8, 1] = -0.11976e-6*u[6, 2]+2.*u[7, 2]-0.11976e-6*u[8, 2]

 

0.11976e-6*u[6, 2]+2.*u[7, 2]+0.11976e-6*u[8, 2] = -0.11976e-6*u[6, 3]+2.*u[7, 3]-0.11976e-6*u[8, 3]

 

0.11976e-6*u[6, 3]+2.*u[7, 3]+0.11976e-6*u[8, 3] = -0.11976e-6*u[6, 4]+2.*u[7, 4]-0.11976e-6*u[8, 4]

 

0.11976e-6*u[6, 4]+2.*u[7, 4]+0.11976e-6*u[8, 4] = -0.11976e-6*u[6, 5]+2.*u[7, 5]-0.11976e-6*u[8, 5]

 

0.11976e-6*u[6, 5]+2.*u[7, 5]+0.11976e-6*u[8, 5] = -0.11976e-6*u[6, 6]+2.*u[7, 6]-0.11976e-6*u[8, 6]

 

0.11976e-6*u[6, 6]+2.*u[7, 6]+0.11976e-6*u[8, 6] = -0.11976e-6*u[6, 7]+2.*u[7, 7]-0.11976e-6*u[8, 7]

 

0.11976e-6*u[6, 7]+2.*u[7, 7]+0.11976e-6*u[8, 7] = -0.11976e-6*u[6, 8]+2.*u[7, 8]-0.11976e-6*u[8, 8]

 

0.11976e-6*u[6, 8]+2.*u[7, 8]+0.11976e-6*u[8, 8] = -0.11976e-6*u[6, 9]+2.*u[7, 9]-0.11976e-6*u[8, 9]

 

0.11976e-6*u[6, 9]+2.*u[7, 9]+0.11976e-6*u[8, 9] = -0.11976e-6*u[6, 10]+2.*u[7, 10]-0.11976e-6*u[8, 10]

 

3.6000 = -0.11976e-6*u[7, 1]+2.*u[8, 1]-0.11976e-6*u[9, 1]

 

0.11976e-6*u[7, 1]+2.*u[8, 1]+0.11976e-6*u[9, 1] = -0.11976e-6*u[7, 2]+2.*u[8, 2]-0.11976e-6*u[9, 2]

 

0.11976e-6*u[7, 2]+2.*u[8, 2]+0.11976e-6*u[9, 2] = -0.11976e-6*u[7, 3]+2.*u[8, 3]-0.11976e-6*u[9, 3]

 

0.11976e-6*u[7, 3]+2.*u[8, 3]+0.11976e-6*u[9, 3] = -0.11976e-6*u[7, 4]+2.*u[8, 4]-0.11976e-6*u[9, 4]

 

0.11976e-6*u[7, 4]+2.*u[8, 4]+0.11976e-6*u[9, 4] = -0.11976e-6*u[7, 5]+2.*u[8, 5]-0.11976e-6*u[9, 5]

 

0.11976e-6*u[7, 5]+2.*u[8, 5]+0.11976e-6*u[9, 5] = -0.11976e-6*u[7, 6]+2.*u[8, 6]-0.11976e-6*u[9, 6]

 

0.11976e-6*u[7, 6]+2.*u[8, 6]+0.11976e-6*u[9, 6] = -0.11976e-6*u[7, 7]+2.*u[8, 7]-0.11976e-6*u[9, 7]

 

0.11976e-6*u[7, 7]+2.*u[8, 7]+0.11976e-6*u[9, 7] = -0.11976e-6*u[7, 8]+2.*u[8, 8]-0.11976e-6*u[9, 8]

 

0.11976e-6*u[7, 8]+2.*u[8, 8]+0.11976e-6*u[9, 8] = -0.11976e-6*u[7, 9]+2.*u[8, 9]-0.11976e-6*u[9, 9]

 

0.11976e-6*u[7, 9]+2.*u[8, 9]+0.11976e-6*u[9, 9] = -0.11976e-6*u[7, 10]+2.*u[8, 10]-0.11976e-6*u[9, 10]

 

3.8000 = -0.11976e-6*u[8, 1]+2.*u[9, 1]-0.59880e-7

 

0.11976e-6*u[8, 1]+2.*u[9, 1]+0.59880e-7 = -0.11976e-6*u[8, 2]+2.*u[9, 2]-0.59880e-7

 

0.11976e-6*u[8, 2]+2.*u[9, 2]+0.59880e-7 = -0.11976e-6*u[8, 3]+2.*u[9, 3]-0.59880e-7

 

0.11976e-6*u[8, 3]+2.*u[9, 3]+0.59880e-7 = -0.11976e-6*u[8, 4]+2.*u[9, 4]-0.59880e-7

 

0.11976e-6*u[8, 4]+2.*u[9, 4]+0.59880e-7 = -0.11976e-6*u[8, 5]+2.*u[9, 5]-0.59880e-7

 

0.11976e-6*u[8, 5]+2.*u[9, 5]+0.59880e-7 = -0.11976e-6*u[8, 6]+2.*u[9, 6]-0.59880e-7

 

0.11976e-6*u[8, 6]+2.*u[9, 6]+0.59880e-7 = -0.11976e-6*u[8, 7]+2.*u[9, 7]-0.59880e-7

 

0.11976e-6*u[8, 7]+2.*u[9, 7]+0.59880e-7 = -0.11976e-6*u[8, 8]+2.*u[9, 8]-0.59880e-7

 

0.11976e-6*u[8, 8]+2.*u[9, 8]+0.59880e-7 = -0.11976e-6*u[8, 9]+2.*u[9, 9]-0.59880e-7

 

0.11976e-6*u[8, 9]+2.*u[9, 9]+0.59880e-7 = -0.11976e-6*u[8, 10]+2.*u[9, 10]-0.11976e-6*u[10, 10]

(2)

``


 

Download Crank_scheme.mw
 

restart; Digits := 5; with(plots); with(LinearAlgebra)

a := 0:

0.11976e-6

(1)

for i from 0 while i <= N do u[i, 0] := h*i+1 end do:

for j from 0 while j < N+1 do u[0, j] := .1; u[N+1, j] := .5 end do:

printlevel := 2:

2.2000 = -0.11976e-7+2.*u[1, 1]-0.11976e-6*u[2, 1]

 

0.11976e-7+2.*u[1, 1]+0.11976e-6*u[2, 1] = -0.11976e-7+2.*u[1, 2]-0.11976e-6*u[2, 2]

 

0.11976e-7+2.*u[1, 2]+0.11976e-6*u[2, 2] = -0.11976e-7+2.*u[1, 3]-0.11976e-6*u[2, 3]

 

0.11976e-7+2.*u[1, 3]+0.11976e-6*u[2, 3] = -0.11976e-7+2.*u[1, 4]-0.11976e-6*u[2, 4]

 

0.11976e-7+2.*u[1, 4]+0.11976e-6*u[2, 4] = -0.11976e-7+2.*u[1, 5]-0.11976e-6*u[2, 5]

 

0.11976e-7+2.*u[1, 5]+0.11976e-6*u[2, 5] = -0.11976e-7+2.*u[1, 6]-0.11976e-6*u[2, 6]

 

0.11976e-7+2.*u[1, 6]+0.11976e-6*u[2, 6] = -0.11976e-7+2.*u[1, 7]-0.11976e-6*u[2, 7]

 

0.11976e-7+2.*u[1, 7]+0.11976e-6*u[2, 7] = -0.11976e-7+2.*u[1, 8]-0.11976e-6*u[2, 8]

 

0.11976e-7+2.*u[1, 8]+0.11976e-6*u[2, 8] = -0.11976e-7+2.*u[1, 9]-0.11976e-6*u[2, 9]

 

0.11976e-7+2.*u[1, 9]+0.11976e-6*u[2, 9] = -0.11976e-6*u[0, 10]+2.*u[1, 10]-0.11976e-6*u[2, 10]

 

2.4000 = -0.11976e-6*u[1, 1]+2.*u[2, 1]-0.11976e-6*u[3, 1]

 

0.11976e-6*u[1, 1]+2.*u[2, 1]+0.11976e-6*u[3, 1] = -0.11976e-6*u[1, 2]+2.*u[2, 2]-0.11976e-6*u[3, 2]

 

0.11976e-6*u[1, 2]+2.*u[2, 2]+0.11976e-6*u[3, 2] = -0.11976e-6*u[1, 3]+2.*u[2, 3]-0.11976e-6*u[3, 3]

 

0.11976e-6*u[1, 3]+2.*u[2, 3]+0.11976e-6*u[3, 3] = -0.11976e-6*u[1, 4]+2.*u[2, 4]-0.11976e-6*u[3, 4]

 

0.11976e-6*u[1, 4]+2.*u[2, 4]+0.11976e-6*u[3, 4] = -0.11976e-6*u[1, 5]+2.*u[2, 5]-0.11976e-6*u[3, 5]

 

0.11976e-6*u[1, 5]+2.*u[2, 5]+0.11976e-6*u[3, 5] = -0.11976e-6*u[1, 6]+2.*u[2, 6]-0.11976e-6*u[3, 6]

 

0.11976e-6*u[1, 6]+2.*u[2, 6]+0.11976e-6*u[3, 6] = -0.11976e-6*u[1, 7]+2.*u[2, 7]-0.11976e-6*u[3, 7]

 

0.11976e-6*u[1, 7]+2.*u[2, 7]+0.11976e-6*u[3, 7] = -0.11976e-6*u[1, 8]+2.*u[2, 8]-0.11976e-6*u[3, 8]

 

0.11976e-6*u[1, 8]+2.*u[2, 8]+0.11976e-6*u[3, 8] = -0.11976e-6*u[1, 9]+2.*u[2, 9]-0.11976e-6*u[3, 9]

 

0.11976e-6*u[1, 9]+2.*u[2, 9]+0.11976e-6*u[3, 9] = -0.11976e-6*u[1, 10]+2.*u[2, 10]-0.11976e-6*u[3, 10]

 

2.6000 = -0.11976e-6*u[2, 1]+2.*u[3, 1]-0.11976e-6*u[4, 1]

 

0.11976e-6*u[2, 1]+2.*u[3, 1]+0.11976e-6*u[4, 1] = -0.11976e-6*u[2, 2]+2.*u[3, 2]-0.11976e-6*u[4, 2]

 

0.11976e-6*u[2, 2]+2.*u[3, 2]+0.11976e-6*u[4, 2] = -0.11976e-6*u[2, 3]+2.*u[3, 3]-0.11976e-6*u[4, 3]

 

0.11976e-6*u[2, 3]+2.*u[3, 3]+0.11976e-6*u[4, 3] = -0.11976e-6*u[2, 4]+2.*u[3, 4]-0.11976e-6*u[4, 4]

 

0.11976e-6*u[2, 4]+2.*u[3, 4]+0.11976e-6*u[4, 4] = -0.11976e-6*u[2, 5]+2.*u[3, 5]-0.11976e-6*u[4, 5]

 

0.11976e-6*u[2, 5]+2.*u[3, 5]+0.11976e-6*u[4, 5] = -0.11976e-6*u[2, 6]+2.*u[3, 6]-0.11976e-6*u[4, 6]

 

0.11976e-6*u[2, 6]+2.*u[3, 6]+0.11976e-6*u[4, 6] = -0.11976e-6*u[2, 7]+2.*u[3, 7]-0.11976e-6*u[4, 7]

 

0.11976e-6*u[2, 7]+2.*u[3, 7]+0.11976e-6*u[4, 7] = -0.11976e-6*u[2, 8]+2.*u[3, 8]-0.11976e-6*u[4, 8]

 

0.11976e-6*u[2, 8]+2.*u[3, 8]+0.11976e-6*u[4, 8] = -0.11976e-6*u[2, 9]+2.*u[3, 9]-0.11976e-6*u[4, 9]

 

0.11976e-6*u[2, 9]+2.*u[3, 9]+0.11976e-6*u[4, 9] = -0.11976e-6*u[2, 10]+2.*u[3, 10]-0.11976e-6*u[4, 10]

 

2.8000 = -0.11976e-6*u[3, 1]+2.*u[4, 1]-0.11976e-6*u[5, 1]

 

0.11976e-6*u[3, 1]+2.*u[4, 1]+0.11976e-6*u[5, 1] = -0.11976e-6*u[3, 2]+2.*u[4, 2]-0.11976e-6*u[5, 2]

 

0.11976e-6*u[3, 2]+2.*u[4, 2]+0.11976e-6*u[5, 2] = -0.11976e-6*u[3, 3]+2.*u[4, 3]-0.11976e-6*u[5, 3]

 

0.11976e-6*u[3, 3]+2.*u[4, 3]+0.11976e-6*u[5, 3] = -0.11976e-6*u[3, 4]+2.*u[4, 4]-0.11976e-6*u[5, 4]

 

0.11976e-6*u[3, 4]+2.*u[4, 4]+0.11976e-6*u[5, 4] = -0.11976e-6*u[3, 5]+2.*u[4, 5]-0.11976e-6*u[5, 5]

 

0.11976e-6*u[3, 5]+2.*u[4, 5]+0.11976e-6*u[5, 5] = -0.11976e-6*u[3, 6]+2.*u[4, 6]-0.11976e-6*u[5, 6]

 

0.11976e-6*u[3, 6]+2.*u[4, 6]+0.11976e-6*u[5, 6] = -0.11976e-6*u[3, 7]+2.*u[4, 7]-0.11976e-6*u[5, 7]

 

0.11976e-6*u[3, 7]+2.*u[4, 7]+0.11976e-6*u[5, 7] = -0.11976e-6*u[3, 8]+2.*u[4, 8]-0.11976e-6*u[5, 8]

 

0.11976e-6*u[3, 8]+2.*u[4, 8]+0.11976e-6*u[5, 8] = -0.11976e-6*u[3, 9]+2.*u[4, 9]-0.11976e-6*u[5, 9]

 

0.11976e-6*u[3, 9]+2.*u[4, 9]+0.11976e-6*u[5, 9] = -0.11976e-6*u[3, 10]+2.*u[4, 10]-0.11976e-6*u[5, 10]

 

3.0000 = -0.11976e-6*u[4, 1]+2.*u[5, 1]-0.11976e-6*u[6, 1]

 

0.11976e-6*u[4, 1]+2.*u[5, 1]+0.11976e-6*u[6, 1] = -0.11976e-6*u[4, 2]+2.*u[5, 2]-0.11976e-6*u[6, 2]

 

0.11976e-6*u[4, 2]+2.*u[5, 2]+0.11976e-6*u[6, 2] = -0.11976e-6*u[4, 3]+2.*u[5, 3]-0.11976e-6*u[6, 3]

 

0.11976e-6*u[4, 3]+2.*u[5, 3]+0.11976e-6*u[6, 3] = -0.11976e-6*u[4, 4]+2.*u[5, 4]-0.11976e-6*u[6, 4]

 

0.11976e-6*u[4, 4]+2.*u[5, 4]+0.11976e-6*u[6, 4] = -0.11976e-6*u[4, 5]+2.*u[5, 5]-0.11976e-6*u[6, 5]

 

0.11976e-6*u[4, 5]+2.*u[5, 5]+0.11976e-6*u[6, 5] = -0.11976e-6*u[4, 6]+2.*u[5, 6]-0.11976e-6*u[6, 6]

 

0.11976e-6*u[4, 6]+2.*u[5, 6]+0.11976e-6*u[6, 6] = -0.11976e-6*u[4, 7]+2.*u[5, 7]-0.11976e-6*u[6, 7]

 

0.11976e-6*u[4, 7]+2.*u[5, 7]+0.11976e-6*u[6, 7] = -0.11976e-6*u[4, 8]+2.*u[5, 8]-0.11976e-6*u[6, 8]

 

0.11976e-6*u[4, 8]+2.*u[5, 8]+0.11976e-6*u[6, 8] = -0.11976e-6*u[4, 9]+2.*u[5, 9]-0.11976e-6*u[6, 9]

 

0.11976e-6*u[4, 9]+2.*u[5, 9]+0.11976e-6*u[6, 9] = -0.11976e-6*u[4, 10]+2.*u[5, 10]-0.11976e-6*u[6, 10]

 

3.2000 = -0.11976e-6*u[5, 1]+2.*u[6, 1]-0.11976e-6*u[7, 1]

 

0.11976e-6*u[5, 1]+2.*u[6, 1]+0.11976e-6*u[7, 1] = -0.11976e-6*u[5, 2]+2.*u[6, 2]-0.11976e-6*u[7, 2]

 

0.11976e-6*u[5, 2]+2.*u[6, 2]+0.11976e-6*u[7, 2] = -0.11976e-6*u[5, 3]+2.*u[6, 3]-0.11976e-6*u[7, 3]

 

0.11976e-6*u[5, 3]+2.*u[6, 3]+0.11976e-6*u[7, 3] = -0.11976e-6*u[5, 4]+2.*u[6, 4]-0.11976e-6*u[7, 4]

 

0.11976e-6*u[5, 4]+2.*u[6, 4]+0.11976e-6*u[7, 4] = -0.11976e-6*u[5, 5]+2.*u[6, 5]-0.11976e-6*u[7, 5]

 

0.11976e-6*u[5, 5]+2.*u[6, 5]+0.11976e-6*u[7, 5] = -0.11976e-6*u[5, 6]+2.*u[6, 6]-0.11976e-6*u[7, 6]

 

0.11976e-6*u[5, 6]+2.*u[6, 6]+0.11976e-6*u[7, 6] = -0.11976e-6*u[5, 7]+2.*u[6, 7]-0.11976e-6*u[7, 7]

 

0.11976e-6*u[5, 7]+2.*u[6, 7]+0.11976e-6*u[7, 7] = -0.11976e-6*u[5, 8]+2.*u[6, 8]-0.11976e-6*u[7, 8]

 

0.11976e-6*u[5, 8]+2.*u[6, 8]+0.11976e-6*u[7, 8] = -0.11976e-6*u[5, 9]+2.*u[6, 9]-0.11976e-6*u[7, 9]

 

0.11976e-6*u[5, 9]+2.*u[6, 9]+0.11976e-6*u[7, 9] = -0.11976e-6*u[5, 10]+2.*u[6, 10]-0.11976e-6*u[7, 10]

 

3.4000 = -0.11976e-6*u[6, 1]+2.*u[7, 1]-0.11976e-6*u[8, 1]

 

0.11976e-6*u[6, 1]+2.*u[7, 1]+0.11976e-6*u[8, 1] = -0.11976e-6*u[6, 2]+2.*u[7, 2]-0.11976e-6*u[8, 2]

 

0.11976e-6*u[6, 2]+2.*u[7, 2]+0.11976e-6*u[8, 2] = -0.11976e-6*u[6, 3]+2.*u[7, 3]-0.11976e-6*u[8, 3]

 

0.11976e-6*u[6, 3]+2.*u[7, 3]+0.11976e-6*u[8, 3] = -0.11976e-6*u[6, 4]+2.*u[7, 4]-0.11976e-6*u[8, 4]

 

0.11976e-6*u[6, 4]+2.*u[7, 4]+0.11976e-6*u[8, 4] = -0.11976e-6*u[6, 5]+2.*u[7, 5]-0.11976e-6*u[8, 5]

 

0.11976e-6*u[6, 5]+2.*u[7, 5]+0.11976e-6*u[8, 5] = -0.11976e-6*u[6, 6]+2.*u[7, 6]-0.11976e-6*u[8, 6]

 

0.11976e-6*u[6, 6]+2.*u[7, 6]+0.11976e-6*u[8, 6] = -0.11976e-6*u[6, 7]+2.*u[7, 7]-0.11976e-6*u[8, 7]

 

0.11976e-6*u[6, 7]+2.*u[7, 7]+0.11976e-6*u[8, 7] = -0.11976e-6*u[6, 8]+2.*u[7, 8]-0.11976e-6*u[8, 8]

 

0.11976e-6*u[6, 8]+2.*u[7, 8]+0.11976e-6*u[8, 8] = -0.11976e-6*u[6, 9]+2.*u[7, 9]-0.11976e-6*u[8, 9]

 

0.11976e-6*u[6, 9]+2.*u[7, 9]+0.11976e-6*u[8, 9] = -0.11976e-6*u[6, 10]+2.*u[7, 10]-0.11976e-6*u[8, 10]

 

3.6000 = -0.11976e-6*u[7, 1]+2.*u[8, 1]-0.11976e-6*u[9, 1]

 

0.11976e-6*u[7, 1]+2.*u[8, 1]+0.11976e-6*u[9, 1] = -0.11976e-6*u[7, 2]+2.*u[8, 2]-0.11976e-6*u[9, 2]

 

0.11976e-6*u[7, 2]+2.*u[8, 2]+0.11976e-6*u[9, 2] = -0.11976e-6*u[7, 3]+2.*u[8, 3]-0.11976e-6*u[9, 3]

 

0.11976e-6*u[7, 3]+2.*u[8, 3]+0.11976e-6*u[9, 3] = -0.11976e-6*u[7, 4]+2.*u[8, 4]-0.11976e-6*u[9, 4]

 

0.11976e-6*u[7, 4]+2.*u[8, 4]+0.11976e-6*u[9, 4] = -0.11976e-6*u[7, 5]+2.*u[8, 5]-0.11976e-6*u[9, 5]

 

0.11976e-6*u[7, 5]+2.*u[8, 5]+0.11976e-6*u[9, 5] = -0.11976e-6*u[7, 6]+2.*u[8, 6]-0.11976e-6*u[9, 6]

 

0.11976e-6*u[7, 6]+2.*u[8, 6]+0.11976e-6*u[9, 6] = -0.11976e-6*u[7, 7]+2.*u[8, 7]-0.11976e-6*u[9, 7]

 

0.11976e-6*u[7, 7]+2.*u[8, 7]+0.11976e-6*u[9, 7] = -0.11976e-6*u[7, 8]+2.*u[8, 8]-0.11976e-6*u[9, 8]

 

0.11976e-6*u[7, 8]+2.*u[8, 8]+0.11976e-6*u[9, 8] = -0.11976e-6*u[7, 9]+2.*u[8, 9]-0.11976e-6*u[9, 9]

 

0.11976e-6*u[7, 9]+2.*u[8, 9]+0.11976e-6*u[9, 9] = -0.11976e-6*u[7, 10]+2.*u[8, 10]-0.11976e-6*u[9, 10]

 

3.8000 = -0.11976e-6*u[8, 1]+2.*u[9, 1]-0.59880e-7

 

0.11976e-6*u[8, 1]+2.*u[9, 1]+0.59880e-7 = -0.11976e-6*u[8, 2]+2.*u[9, 2]-0.59880e-7

 

0.11976e-6*u[8, 2]+2.*u[9, 2]+0.59880e-7 = -0.11976e-6*u[8, 3]+2.*u[9, 3]-0.59880e-7

 

0.11976e-6*u[8, 3]+2.*u[9, 3]+0.59880e-7 = -0.11976e-6*u[8, 4]+2.*u[9, 4]-0.59880e-7

 

0.11976e-6*u[8, 4]+2.*u[9, 4]+0.59880e-7 = -0.11976e-6*u[8, 5]+2.*u[9, 5]-0.59880e-7

 

0.11976e-6*u[8, 5]+2.*u[9, 5]+0.59880e-7 = -0.11976e-6*u[8, 6]+2.*u[9, 6]-0.59880e-7

 

0.11976e-6*u[8, 6]+2.*u[9, 6]+0.59880e-7 = -0.11976e-6*u[8, 7]+2.*u[9, 7]-0.59880e-7

 

0.11976e-6*u[8, 7]+2.*u[9, 7]+0.59880e-7 = -0.11976e-6*u[8, 8]+2.*u[9, 8]-0.59880e-7

 

0.11976e-6*u[8, 8]+2.*u[9, 8]+0.59880e-7 = -0.11976e-6*u[8, 9]+2.*u[9, 9]-0.59880e-7

 

0.11976e-6*u[8, 9]+2.*u[9, 9]+0.59880e-7 = -0.11976e-6*u[8, 10]+2.*u[9, 10]-0.11976e-6*u[10, 10]

(2)

``


 

Download Crank_scheme.mw

 

1 2 Page 2 of 2