summaryrefslogtreecommitdiff
path: root/377/CH1/EX1.5/1_5.sce
blob: 5a770cfb5e3a7697223bc6490d909745c0295f52 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
disp("T=16*exp(-2αβ/((1+β/α)^2)*((1+α/β)^2)))");
E=1;
V0=4;
pi=3.14;
m0=1.67*(10^-27)*1.6*(10^-13);
h=1.055*10^-34;
disp("β^2=4*pi^2*m0*(V0-E)/(h^2)");
printf('\n');
disp("α^2=4*pi^2*m0*E/(h^2)");
b1=sqrt(2*m0*(V0-E)/(h^2));
printf('\n The value of β is %f*10^14',b1*10^-14);
c=3^0.5;  //say c=β/α
a=1.9*10^14;
b=10^-14;
T=16*(exp(((-2)*a*b))/(((1+c)^2)*((1+(1/c))^2)));
printf('\n The value of T is %fpercent',T*100);