// I would like this to be a long double,
// but the c library only provides
// plain double trig functions etc.
typedef double daton;