#include float octcps(cps) float cps; { return(log(cps/1.021975)/.69314718); /* .69etc is log of 2., 1.02etc is offset for middle C */ }