% solves the boyce problem % requires the m-file retire.m % uses starting guess rho=0.1 rho = fzero(@retire , 0.1)