> restart;

> with( Student[Calculus1] );

[AntiderivativePlot, AntiderivativeTutor, ApproximateInt, ApproximateIntTutor, ArcLength, ArcLengthTutor, Asymptotes, Clear, CriticalPoints, CurveAnalysisTutor, DerivativePlot, DerivativeTutor, DiffTu...[AntiderivativePlot, AntiderivativeTutor, ApproximateInt, ApproximateIntTutor, ArcLength, ArcLengthTutor, Asymptotes, Clear, CriticalPoints, CurveAnalysisTutor, DerivativePlot, DerivativeTutor, DiffTu...[AntiderivativePlot, AntiderivativeTutor, ApproximateInt, ApproximateIntTutor, ArcLength, ArcLengthTutor, Asymptotes, Clear, CriticalPoints, CurveAnalysisTutor, DerivativePlot, DerivativeTutor, DiffTu...[AntiderivativePlot, AntiderivativeTutor, ApproximateInt, ApproximateIntTutor, ArcLength, ArcLengthTutor, Asymptotes, Clear, CriticalPoints, CurveAnalysisTutor, DerivativePlot, DerivativeTutor, DiffTu...[AntiderivativePlot, AntiderivativeTutor, ApproximateInt, ApproximateIntTutor, ArcLength, ArcLengthTutor, Asymptotes, Clear, CriticalPoints, CurveAnalysisTutor, DerivativePlot, DerivativeTutor, DiffTu...

>

> f := x/(x^4+1);

f := x/(x^4+1)

>

> A := 0;

A := 0

> B := 4;

B := 4

> N := 4;

N := 4

>

> ApproximateInt(f, A..B, partition = N, method = left, output = sum);

Sum(i/(i^4+1), i = 0 .. 3)

> ApproximateInt(f, A..B, partition = N, method = left, output = plot);

[Plot]

>

> ApproximateInt(f, A..B, partition = N, method = midpoint, output = sum);

Sum((i+1/2)/((i+1/2)^4+1), i = 0 .. 3)

> ApproximateInt(f, A..B, partition = N, method = midpoint, output = plot);

[Plot]

>

> ApproximateInt(f, A..B, partition = N, method = right, output = sum);

Sum(i/(i^4+1), i = 1 .. 4)

> ApproximateInt(f, A..B, partition = N, method = right, output = plot);

[Plot]

>

> ApproximateInt(f, A..B, partition = N, output = animation);

[Plot]

>

> ApproximateInt(f, A..B, partition = n, method = left, output = sum);

4*(Sum(4*i/(n*(256*i^4/n^4+1)), i = 0 .. n-1))/n

> limit( %, n=infinity );

1/2*arctan(16)

> evalf( % );

.7541887585

>

> q1 := ApproximateInt(f, a..b, partition = n, method = left, output = sum);

q1 := (b-a)*(Sum((a+i*(b-a)/n)/((a+i*(b-a)/n)^4+1), i = 0 .. n-1))/n

>

> #value( % );

> #simplify( %, n );

> #collect( %, n );

> q2 := limit( q1, n=infinity );

q2 := -1/2*arctan(a^2)+1/2*arctan(b^2)

>

> F := int( f, x );

F := 1/2*arctan(x^2)

> q3 := eval( F, x=b ) - eval( F, x=a );

q3 := -1/2*arctan(a^2)+1/2*arctan(b^2)

> q4 := int( f, x=a..b );

q4 := -1/2*arctan(a^2)+1/2*arctan(b^2)

>

>