%% example 4.3.1 q(40) = 2.7812; q(39) = 2.5982; q(38) = 2.4313; q(37) = 2.2791; q(36) = 2.1402; q(35) = 2.0136; q(34) = 1.8980; q(33) = 1.7927; q(32) = 1.6965; q(31) = 1.6089; q(30) = 1.5289; qMin = 30; q = q/1000;