clear all

%변수들

boltzmann_const = 1.38e-23;

planck_const = 6.62e-34;

hbar = planck_const / ( 2 *pi );

transition_temp = 9.2;

gap_energy_at_zero_kelvin = 3.528 / ( 2 * transition_temp * boltzmann_const );

debye_freq = ( 296 * boltzmann_const ) / hbar; %디바이 진동수

a_const = log( ( 1.13 * hbar * debye_freq ) / ( boltzmann_const * transition_temp) );

b_const = ( hbar * debye_freq ) / ( 2 * boltzmann_const);

%%%여기서부터 공식

function fun=gap(C,X)

tanh(1477.92*sqrt(C*C+X*X)/sqrt(C*C+X*X));

end %여기가 문제항

for i=linspace(0.001,1);

X = i;

integral =quad(@gap(C,X),0,1)

delta_val=fsolve(integral,[0,0])-a_const

end

plot(delta,i)


내 매트랩 파일인데


다른게 문제가 아니고


저 X라는 함수를 for문 안에 넣어서


저 범위내로 한번 돌려보고 싶거든


파이썬에서는


for

def


하면되는데 매트랩에서는 없으니까 어떻게 할 방법 없음??