Untitled

X = @(t) sin(sqrt(t)) .* exp(-2 .* t.^2); % Funzione da integrare
concentrazione = @(t) arrayfun(@(T) integral(X, 0, T), t); % Definire l'integrale dipendente da t
t_val = 1; % Ad esempio, t = 1
disp(concentrazione(t_val));
t_vals = linspace(0, 2, 100); % Valori di t tra 0 e 2
concentrazioni = concentrazione(t_vals); % Calcolare per ogni valore di t
disp(concentrazioni);