Numerical integration with many parameters.

1 vue (au cours des 30 derniers jours)
JACINTA ONWUKA
JACINTA ONWUKA le 4 Juin 2019
Modifié(e) : Torsten le 4 Juin 2019
I will like to have the values of n2 depending on 10 values of Mycp. please i need help with this my code below.
L=1;
T=100;
r=0.03;
I1=0.5;
p=0.005;
epsilon=0.5;
beta=0.1;
rho=0.5;
Mycp=0:1:10;
delta=1-Mycp/100
tau=(1/(beta*(L+delta*p)))*log((L*(I1+delta*p))/(delta*p*(L-I1)));
t05 =(1/(beta*(L+delta*p)))*log((L*(0.05*L+delta*p))/(delta*p*(L-0.05*L)));
I2= (L*delta*p*(exp (beta*(L+delta*p)*t)-1)) ./ (L + delta*p* exp(beta*(L+delta*p)*t));
fun2=@(t,rho,I2,r)rho*I2*exp(-r*t);
n2= integral(@(t)fun2(t,rho,I2,r),t05,tau, 'ArrayValued',1)

Réponse acceptée

Torsten
Torsten le 4 Juin 2019
Modifié(e) : Torsten le 4 Juin 2019
function main
L = 1;
T = 100;
r = 0.03;
I1 = 0.5;
p = 0.005;
epsilon = 0.5;
beta = 0.1;
rho = 0.5;
Mycp = 0:1:10;
n2 = zeros(numel(Mycp),1);
for i = 1:numel(Mycp)
delta = 1-Mycp(i)/100;
tau = (1/(beta*(L+delta*p)))*log((L*(I1+delta*p))/(delta*p*(L-I1)));
t05 =(1/(beta*(L+delta*p)))*log((L*(0.05*L+delta*p))/(delta*p*(L-0.05*L)));
I2= @(t)(L*delta*p*(exp (beta*(L+delta*p)*t)-1)) ./ (L + delta*p* exp(beta*(L+delta*p)*t));
fun2 = @(t)rho*I2(t).*exp(-r*t);
n2(i) = integral(fun2,t05,tau);
end
plot(Mycp,n2)
end

Plus de réponses (1)

darova
darova le 4 Juin 2019
You forgot about element-wise operator
delta=1-Mycp./100;
Look also for VECTORIZE

Catégories

En savoir plus sur MATLAB 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!

Translated by