Integration of two functions, one of them is piecewise function
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Consider the function
p(µ, σ, t, x) = exp(-m*t)/x*sqrt(2*pi*(q.^2)*t)*exp(-1/(2*(q.^2)*t)*(log(x)-(m-(q.^2)/2)*t).^2) (0 < x < ∞)
Here let µ = .02, σ = .18 and t = .25, so p is a function on the variable x only. I created a function file pfun.m for this function, having µ, σ, t and x as input.
I need to compute int(C(x)*p(x),x,0,inf)
where C(x) = 0 if x < 1.02
x − 1.02 otherwise
0 commentaires
Réponses (1)
Walter Roberson
le 29 Sep 2015
Modifié(e) : Walter Roberson
le 29 Sep 2015
general_p = @(m, q, t, x) exp(-m.*t) ./ x .*sqrt(2*pi*(q.^2).*t) .* exp(-1/(2*(q.^2).*t) .* (log(x)-(m-(q.^2)/2) .*t ).^2);
p = @(x) general_p(0.02, 0.18, 0.25, x);
C = @(x) (x - 1.02) .* (1.02 <= x);
result = integral(@(x) C(x) .* p(x), 0, inf)
3 commentaires
Voir également
Catégories
En savoir plus sur Assumptions dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!