added general lib
[physik/computational_physics.git] / general.c
1 /* general functions */
2
3 double fak(int l) {
4         double result=1;
5         while(l--) result*=l;
6         return result;
7 }
8
9 double fak2(int l) {
10         double result=1;
11         while(l>0) {
12                 result*=l;
13                 l-=2;
14         }
15         return result;
16 }
17
18 double absolute_value(double l) {
19         return l<0?-l:l;
20 }