>

restart;

>	a:=0:b:=1.:        # Интервал

>	eq:=diff(y(t),t)+t;# Ур-е

>	dsolve({eq,y(0)=1},y(t));# Можно решить точно

>	int(1-t^2/2,t=a..b);# Интеграл от решения (точно)

>	R:=dsolve({eq,y(0)=1},y(t),numeric,output=listprocedure);# Численно интегрируем

>	Y:=eval(y(t),R);

>	Y(0); #Проверка

>	plot(Y(t),t=a..b); # График

>	N:=100:#Интегрируем методом прямоугольников (можно Симпсоном при желании)

>	dx:=(b-a)/N:

>	INT:=add(dx*Y(i/N),i=1..N);# Интеграл от решения (приближенно)