matlab, интеграл с параметром, график

[fermach]

вопрос по matlab
подскажите, пожалуйста, с помощью каких функций можно нарисовать график такой вот штуки
f(t)=\int_0^t \eta(t-p) exp (-x^2/4p) dp
интеграл не берется явно
спасибо

[photon]

x - в данном случае константа?

[fermach]

здесь можно считать константой
итогом будет график функции двух переменных (t,x).

[photon]

Дайте требуемый диапазон t и х, и я построю Вам этот график

[fermach]

x = 0:20
t = 0:10
интересен не сколько график, а как
спасибо

[photon]

Вечерком нарисовал я эту функцию... но - хотел спать, поэтому:
а) Не перепроверял - доверяй, но проверяй
б) Не писал комментарии: если что-то будет непонятно - добро пожаловать.
Выбросил я эту $\eta$ - потом домножить на константу можно всегда, и получил такую функцию для вычисления интеграла:

function integ_par=integ_par(x,t,number_of_points)
fxt=zeros(1,number_of_points);
integ_par=zeros(length(t),length(x));
for ii=1:length(t)
    p=0:t(ii)/number_of_points:t(ii);
    for jj=1:length(x)
        if ~isempty(p)   
            fxt=(t(ii)-p(2:end)).*exp(-0.25*x(jj)^2./p(2:end));
            integ_par(ii,jj)=0.5*sum((fxt(1:end-1)+fxt(2:end))).*t(ii)/number_of_points;
        end
    end
end

В данном случае вычисление ведется методом трапеций.

Вызов функции и построение графика осуществляется например так:

x=0:.1:20;
t=0:.1:10;
integ=integ_par(x,t,50);
[tt,xx]=meshgrid(t,x);
mesh(t,x,integ)

В результате имеем: