Solving a system of two integral equations numerically
1 回表示 (過去 30 日間)
古いコメントを表示
I am trying to solve a system of two equations for
given that
and letting
,
vary from 0.1 to 1 in steps of 0.1.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581966/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581971/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581976/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581981/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581986/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581991/image.png)
I think I need to do a double for loop that numerically approximates the integral and then solves for appropriate
values.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/581996/image.png)
I did make an attempt but it is not going so well.
count= 1; %counter
svals1 = zeros(10,20); %an array
svals2 = zeros(10,20); %an array
g = 0.65; %gamma
d1= 1; %delta1
b1 = 1.1; %beta1
a11= 1; %alpha_11
a12 = 1; %alpha_12
a21 = 1; %alpha_21
a22 = 1; %alpha_22
syms x;
syms S1;
syms S2;
for d2= 0:0.1:count
for b2 = 0:0.1:count
svals1(d2,b2,x,S1,S2) = vpasolve( (b1.*g)/d1.*integral( exp(- g.*x).* (1-exp(-d1.* x)).* exp(-(a11.*(b1.*S1)/d1 (1-exp(-d1.* x) ) +a12.* (b2.*S2)/d2 (1-exp(-d2.* x)))),0,Inf)-1,S1)
svals2(d2,b2,x,S1,S2) = vpasolve( (b2.*g)/d2.*integral( exp(- g.* x).* (1-exp(-d2.* x)) .*exp(-(a22.*(b2.*S2)/d2 (1-exp(-d2.* x) ) +a21.* (b1.*S1)/d1 (1-exp(-d1.* x)))),0,Inf)-1,S2)
disp(d2,b2,svals1,svals2)
end
end
0 件のコメント
回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Calculus についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!