Not able to solve this expression

2 次查看(过去 30 天)
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 个评论
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 CenterFile Exchange 中查找有关 Logical 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by