1 #include "f2c.h"
2 
3 double
r_dim(real * a,real * b)4 r_dim (real * a, real * b)
5 {
6   return (*a > *b ? *a - *b : 0);
7 }
8