Not able to solve this expression

2 次查看（过去 30 天）
Anshuman S
Anshuman S 2018-10-10
评论： KALYAN ACHARJYA 2018-10-10
I am not able to code this expression , I have tried but it's not giving correct answer ! My try is given below:
if true
sol= zeros(4,4);
s=0;
for i = 1:5
for j = 1:5
r = (j-1)*0.25 ;
fi = 2*pi*(i-1)/4 ;
disp(r)
disp(fi)
disp('r,fi')
for m1 =1:100000
s = s + (2*((r/0.03)^(2*m1))*cos(2*m1*fi)) / ((4*m1*m1 - 1)*(1 + (2*m1*40)/(20*0.03)));
end
sol(i,j)= 1/pi - (1/pi)*s + ((r/0.03)*sin(fi))/(2*(1+ (45/(20*0.03))));
end
end
disp(sol)
end
The r & phi are taken for 5 and 5 points respectively . r = 0, 0,25, 0,50, 0,75 1.0 and simillarly phi = 0, pi/2 , pi , 3pi/2 , 2pi . But my solution is coming very large than the expected value.
  4 个评论
显示 2更早的评论
madhan ravi
madhan ravi 2018-10-10
" for m1 =1:100000 aftermath due to this leads to more values "
Anshuman S
Anshuman S 2018-10-10
it's given in the problem to take many number of terms, m1 shows the number of terms in the series, in the problem it's taken to be infinite

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Anshuman S
Anshuman S 2018-10-10
I have solved the problem now , The r term was not within the bounds. It should be less than r0 ( 0.03 cm ) . thus r = 0.03*(j-1)*0.25 ;
Thanks for your time and help !

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by