LDLIBS=-lm -lgsl -lgslcblas calculate_limits: calculate_limits.c