Matlab problem! Please help...

[Chest]

I'm now writing a program to calculate the input resistance in Matlab code. But i'm faced some problem. May anyone can check it or correct it for me? The equations are

Thanks in advance.

Maltlab codes that i writing (following parts is failure to running)

c=(ko*W/2).*cos(t);

F1 = inline('sin(c).^2.*sin(t).^3.*sin(t)./cos(t).^2');
I1 = quad(F1,0,pi);

I12=sum(sin(c).^2.*sin(t).^3.*sin(t)./cos(t).^2.*besselj(0,(ko*L).*sin(t)));
G1=I1./(120*pi^2);
G12=I12./(120*pi^2);
Rin0=(2.*(G1+G12))^-1;

 

[mdanaie]

it is obvious why you don't find the answere because F1 has two input variables: t and c but in quad you change only one varible from 0 to pi

I reccomend you calculate the integral by using "sum/ delta_t" there is no need to use inline or quad

 

[Chest]

Hi, mdanaie:

It work well after I changing the integral by using "sum/ delta_t". Awesome. Thanks a lot.

Best Regards