16.1	b) 	plot x,3-x
		integral 3-2x dx from 0 to 1.5
	c) 	plot[{x^2/2+2,4},{x,0,2}]
		integral 4-(x^2/2+2) dx from 0 to 2
	f) 	plot[{x^2,x^3},{x,0,1}]
		integral x^2-x^3 dx from 0 to 1
	g) 	plot[{3/x^2,0},{x,2,3}]
		integral 3/x^2 dx from 2 to 3
	h) 	plot[{sqrt(x-1),3-x,0},{x,1,3}]
		integral sqrt(x-1) dx from 1 to 2 + integral 3-x dx from 2 to 3
	i) 	plot[{x^5/2,0},{x,-1,2},{y,-1,16}]
		-integral x^5/2 dx from -1 to 0 + integral x^5/2 dx from 0 to 2
	j) 	plot[{x^2,-(-x)^(1/3),x^(1/3)},{x,-2,2}]
		integral x^2-x^{1/3} dx from -2 to 0 + integral x^{1/3}-x^2 dx from 0 to 1+integral x^2-x^{1/3} dx from 1 to 2
	k) 	plot[{sqrt(x^4(1-x^3)),-sqrt(x^4(1-x^3))},{x,-1,2}]
		2*integral x^2*sqrt(1-x^3) dx from 0 to 1

16.2	polar plot r=2sin theta where theta from pi/4 to pi/2
	1/2*integral 4sin^2x dx from pi/4 to pi/2

16.5	PolarPlot[{r=2sin (3*theta),r=2|sin (3*theta)|},{theta,0,pi}]
	1/2*integral 4sin^2(3x) dx from 0 to pi	

16.6	PolarPlot[{r=1+cos theta,r=cos theta},{theta,0,2pi}]
	PolarPlot[{r=1+cos theta,r=cos theta},{theta,0,pi}]
	1/2*(integral (1+cos x)^2 dx from 0 to 2pi) - 1/2*(integral cos^2 x dx from 0 to pi)
	
16.7	a) 	plot x^(1/3)/3+1, x=0..2
		integral sqrt(1+x/4 )dx from 0 to 2
	b)	plot x^(2/3), y=0..4
		integral sqrt(1+9x/4 )dx from 0 to 4

16.9	ParametricPlot[{x=cos (theta),y=sin (theta)},{theta,-9/4pi,-5/4pi}]
	integral r dx from -9/4pi to -5/4pi
	
16.10	ParametricPlot[{x=cos^3 (theta),y=sin^3 (theta)},{theta,0pi,2pi}]
	4* integral 3a*sin x * cos x dx from 0 to pi/2
	
16.12	integral pi*(1-x^2)^2 dx from -1 to 1
	
16.14	integral 1/2 * 3x/5 * 4x/5 dx from 0 to 5

16.15	integral (2*sqrt(r^2-x^2))^2 dx from -r to r (bisilinder ehk Steinmetzi keha)
	
16.17	c) 	plot[{x^2+1,x+1},{x,0,1}]
		pi * integral (x+1)^2-(x^2+1)^2 dx from 0 to 1
	e)	plot[{cos x},{x,0,pi}]
		pi * integral (cos x)^2 dx from 0 to pi

16.19	c)	plot[{3sqrt(x),3},{x,0,1}]
		pi * integral x^4/81 dx from 0 to 3
	e)	plot[{sqrt(4-x^2)},{x,0,pi}]
		pi * integral 4-x^2 dx from 0 to 2
	
16.21	plot[2x^2,x+1},{x,-0.5,1}]
	pi * integral (x+1)^2-4x^4 dx from -0.5 to 1
	
16.25	integral_0^20 sqrt(400-20t) dt
	
16.26	Plot[Piecewise[{{2000t-1000t^2, 0.1>t>0.0}, {190, 0.3>t>0.1}, {190-900(t-0.3), 0.4>t>0.30}}], {t, 0, 0.4}]
	integral 2000t-1000t^2 dt from 0 to 0.1 + integral 190 dt from 0.1 to 0.3 + integral 190-900t dt from 0 to 0.1
		
16.27	plot[{t^3-6t^2+8t},{t,0,6}]
	integral |t^3-6t^2+8t| dt from 0 to 6 
	integral t^3-6t^2+8t dt from 0 to 6 
	
16.28	Plot[Piecewise[{{5t, 4>t>0}, {20, 10>t>4}}], {t, 0, 10}]
	s=integral_0^4 a*t dt
	a=56*a/(integral_0^4 a*t dt + integral_4^10 4a dt)

16.29	a) 1950-integral_0^4 9.8*t dt
	b) 2130-integral_0^13 9.8*t dt
	c) 1950-integral_0^4 9.8*t dt - integral_4^(45+13) 5 dt > 2130-integral_0^13 9.8*t dt

16.31	integral t/(t^2+1) dt from 0 to 5
	
16.34	integral 80/3-250t dt from 0 to 8/75

16.47	integral 640-64x dx from 0 to 5

16.50	a) 	plot[{4-x,0},{x,0,4}]
		integral 4-x dx from 0 to 4
		1/8 * integral x(4-x) dx from 0 to 4
		1/8 * integral (4-x)^2/2 dx from 0 to 4
	b)	plot[{x^3,0},{x,0,2},{y,0,8}]
		integral x^3 dx from 0 to 2
		1/4 * integral x^4 dx from 0 to 2
		1/8 * integral x^6 dx from 0 to 2