
w = [3.5 4 6 25 103 117 200 252 300 437 1340 2000 2700 5000 22500 30000 33000 50000 70000 100000 415000 450000 500000 3000000];
r = [787 660 588 670 347 300 420 352 300 269 251 205 187 120 100 85 81 70 72 70 45 38 40 48];

title('mammalian heart rate model')
xlabel('weight w^(-1/3)')
ylabel('pulse rate')
plot(w.^(-1/3),r,'o')
