Mathematical Physics - Volume II - Numerical Methods

Chapter 3. Comparison of finite element method and finite difference method

86

----------------------------------------------------------- 1 1 0 9 1 1 3 9 7 2 6 8 4 5 ----------------------------------------------------------- 1 -1. 0. 4. ----------------------------------------------------------- 0 ----------------------------------------------------------- 1 1 1 0 0. ----------------------------------------------------------- 4 1 1 0 1 0. 0. 1 1 0 2 0. -4. 1 1 0 3 0. -4. 1 1 0 4 0. 0. ----------------------------------------------------------- end ----------------------------------------------------------- After running the program, the “fem.out” database is formed:

Test primer br. 1: Jednacina Lap(u)=4. Resenje: u=x^2+y^2. node no x-coordinate y-coordinate

1 2 3 4 5 6 7 8 9

.0000

.0000 .0000 .0000

1.0000 2.0000

.0000

1.0000 1.0000 1.0000 2.0000 2.0000 2.0000

1.0000 2.0000

.0000

1.0000 2.0000

elem no ne mat

node numbers

1

9

1

1 3 9 7 2 6 8 4 5

mat no

k

b

f

1

-1.0000

.0000

4.0000

node no

essent. value

1

.0000

element no side

p

gamma