The following is the code I used –
q=10e8:(20e9-10e8)/1000:20e9;q=q';s=length(q);epso = 8.8 *10.^-12;epsr=5.7;ro=3.63*(10.^-3);epsl = 2.5;ri=1.13*(10.^-3);copen = (epso/sqrt(epsr)*log(ro/ri));cload= (epsl/sqrt(epsr)*log(ro/ri));t = zeros(1,s);for f = 1:length(q) fun = @(x)((((((besselj(0,(2*pi*q(f))/(3e8))*x*ro)-besselj(0,(2*pi*q(f))/(3e8))*x*ri)).^2)./x)*freload(x,q)); y(f) = copen*abs(integral(fun,0,Inf,'Arrayvalued',true)); t(f)=(1-y(f))/(1+y(f));endplot (q,t);
where the function freload is as follows –
function f = freload(x,q)epsl=3.4-0.001j;f=(1./sqrt(epsl-(x.^2))).*(1+exp(-2*1i*2*pi*q./(3e8).*sqrt(epsl-(x.^2))))./((1-exp(-2*1i*2*pi*q/(3e8).*sqrt(epsl-(x.^2)))));end
The error i get is as follows – Subscripted assignment dimension mismatch.
Error in numint (line 14) y(f) = copen*abs(integral(fun,0,Inf,'arrayvalued',true));
Any help or suggestion would be greatly appreciated. Thanks
Best Answer