> | 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);# Интеграл от решения (приближенно) |