clc; clear;
%case1
syms v w theta Xrx Yrx Zrx y d %v는 대문자파이 w는 소문자 파이
d = 0.15;
r1 = 0.5*d; % rx반지름
r2 = 0.075; % tx반지름
y = 0;

Ntx = 18; % tx권선수
Nrx = 18; % rx권선수
Ltx = 217.5*10^-6 % 인덕터 값들
Lrx = 80.3*10^-6
u0 = 4*pi*10^-7 %자유공간 투자율
T = [1 0 0; 0 cos(theta) -sin(theta); 0 sin(theta) cos(theta)] % 변환행렬
RX = T*[r2*cos(w);r2*sin(w);0]+[0;y;d] % 변환 한 공식 좌표
TX = [r1*cos(v); r1*sin(v); 0]
dltx = diff(TX,v) % TX v에 대해 미분
dlrx = diff(RX,w) % RX w에 대해 미분
R = ((RX(1)-TX(1))^2 + (RX(2)-TX(2))^2 + (RX(3)-TX(3))^2)^(1/2) %좌표사이의 거리
dltxrx = @(v,w)dltx'*dlrx % 상호인덕턴스 적분에 쓰일 값
eqn = dltxrx/R;  % 상호인덕턴스 적분에 쓰일 값
eqn1 = int(eqn,v,0,2*pi)
eqn2 = int(eqn1,w,0,2*pi)
eqn3 = ((Ntx*Nrx*u0)/(4*pi*(Ltx*Lrx)^(1/2))) % k(theta)에서 알파1 계산
k(theta) = eqn2*eqn3 % 14번식



여기서 다음 식을 실행해보면 적분값이 그냥


eqn1 =
 
int(((9*sin(conj(v))*sin(w))/1600 + (9*cos(conj(v))*cos(theta)*cos(w))/1600)/(((3*sin(v))/40 - (3*cos(theta)*sin(w))/40)^2 + ((3*cos(v))/40 - (3*cos(w))/40)^2 + ((3*sin(theta)*sin(w))/40 + 3/20)^2)^(1/2), v, 0, 2*pi)
 
 
eqn2 =
 
int(int(((9*sin(conj(v))*sin(w))/1600 + (9*cos(conj(v))*cos(theta)*cos(w))/1600)/(((3*sin(v))/40 - (3*cos(theta)*sin(w))/40)^2 + ((3*cos(v))/40 - (3*cos(w))/40)^2 + ((3*sin(theta)*sin(w))/40 + 3/20)^2)^(1/2), v, 0, 2*pi), w, 0, 2*pi)
 

이대로 int값 살아있는 체로 나오네요.. 왜이렇게 되는 건가요

??