*Hi I am working on the following problem (See below the line) from the text A Introduction to the Mathematics of Biology (Ch 2 pgs.21-23). I am working in Maple 18 and the code for this problem is from Maple 11 I think. If you look at the last line of code, I ran into a snag. Can anyone help????? We are in a study Gourp and alll are stumped.*

*_______________________________________________________________________________________________________*

**Problem # 2 - Find a fit for the cumulative US AIDS data as a polynomial function. Alsoo find an exponential fit for the data. (assume the Data I have inputed is correct if you don't have the book).**

AIDS := [97, 206, 406, 700, 1289, 1654, 2576, 3392, 4922, 6343, 8359, 9968, 12990, 14397, 16604, 17124, 19585, 19707, 21392, 20846, 23690, 24610, 26228, 22768];

print(`output redirected...`); # input placeholder

[97, 206, 406, 700, 1289, 1654, 2576, 3392, 4922, 6343, 8359,

9968, 12990, 14397, 16604, 17124, 19585, 19707, 21392, 20846,

23690, 24610, 26228, 22768]

CAC := [seq(sum(AIDS[j]/(1000.0), j = 1 .. i), i = 1 .. 24)];

print(`output redirected...`); # input placeholder

[0.09700000000, 0.3030000000, 0.7090000000, 1.409000000,

2.698000000, 4.352000000, 6.928000000, 10.32000000,

15.24200000, 21.58500000, 29.94400000, 39.91200000,

52.90200000, 67.29900000, 83.90300000, 101.0270000,

120.6120000, 140.3190000, 161.7110000, 182.5570000,

206.2470000, 230.8570000, 257.0850000, 279.8530000]

Time := [seq(1981+(i-1)*(1/2), i = 1 .. 24)];

ln(CAC)=k*lnt+A

LnCAC := map(ln, CAC);

print(`output redirected...`); # input placeholder

[-2.333044300, -1.194022473, -0.3438997525, 0.3428802329,

0.9925107578, 1.470635510, 1.935571171, 2.334083760,

2.724054775, 3.071998629, 3.399328971, 3.686677031,

3.968441145, 4.209145378, 4.429661370, 4.615387808,

4.792578782, 4.943918402, 5.085810791, 5.207062453,

5.329074480, 5.441798471, 5.549406770, 5.634264465]

LnTime := map(ln, [seq((i+1)/(2*(1/10)), i = 1 .. 24)]);

print(`output redirected...`); # input placeholder

[ln(10), ln(15), ln(20), 2 ln(5), ln(30), ln(35), ln(40), ln(45),

ln(50), ln(55), ln(60), ln(65), ln(70), ln(75), ln(80), ln(85),

ln(90), ln(95), 2 ln(10), ln(105), ln(110), ln(115), ln(120),

3 ln(5)]

with(stats);

fit[leastsquare[[x, y], y = k*x+lnA]]([LnTime, LnCAC]);

print(`output redirected...`); # input placeholder

y = 3.293411005 x - 10.12289000

k := op(1, op(1, rhs(%))); LnA := op(2, rhs(`%%`)); A := exp(LnA);

**print(`output redirected...`); # input placeholder**

**Error, invalid input: rhs received exp(LnA), which is not valid for its 1st argument, expr**

**Error, invalid input: rhs received exp(LnA), which is not valid for its 1st argument, expr**

** exp(LnA)**