summaryrefslogtreecommitdiff
path: root/2453/CH9
diff options
context:
space:
mode:
Diffstat (limited to '2453/CH9')
-rwxr-xr-x2453/CH9/EX9.1/9_1.sce6
-rwxr-xr-x2453/CH9/EX9.2/9_2.sce11
-rwxr-xr-x2453/CH9/EX9.3/9_3.sce6
-rwxr-xr-x2453/CH9/EX9.4/9_4.sce9
4 files changed, 32 insertions, 0 deletions
diff --git a/2453/CH9/EX9.1/9_1.sce b/2453/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..f64067a45
--- /dev/null
+++ b/2453/CH9/EX9.1/9_1.sce
@@ -0,0 +1,6 @@
+//To calculate the critical field
+Tc = 3.7; //critical temperature, K
+Hc_0 = 0.0306; //critical field, T
+T = 2; //temperature, K
+Hc_2 = Hc_0*(1-(T/Tc)^2); //critical field, T
+printf("critical field at 2K is %5.5f T",Hc_2);
diff --git a/2453/CH9/EX9.2/9_2.sce b/2453/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..7ef34b6c7
--- /dev/null
+++ b/2453/CH9/EX9.2/9_2.sce
@@ -0,0 +1,11 @@
+//To calculate the critical current
+T = 4.2; //temperature, K
+d = 1; //diameter, mm
+d = d*10^-3; //diameter, m
+Tc = 7.18; //critical temperature, K
+H0 = 6.5*10^4; //critical field, A/m
+Hc = H0*(1-(T/Tc)^2); //critical field at 2K, A/m
+ic = %pi*d*Hc; //critical current, A
+printf("critical current for lead is %5.2f A",ic);
+
+//answer given in the book is wrong
diff --git a/2453/CH9/EX9.3/9_3.sce b/2453/CH9/EX9.3/9_3.sce
new file mode 100755
index 000000000..9a01d99fc
--- /dev/null
+++ b/2453/CH9/EX9.3/9_3.sce
@@ -0,0 +1,6 @@
+//To calculate the penetration depth
+lamda_T = 750; //penetration depth of mercury, Armstrong
+T = 3.5; //temperature, K
+Tc = 4.12; //critical temperarure, K
+lamda_0 = lamda_T*((1-(T/Tc)^4))^(1/2); //penetration depth, Armstrong
+printf("penetration depth at 0K is %d armstrong",lamda_0);
diff --git a/2453/CH9/EX9.4/9_4.sce b/2453/CH9/EX9.4/9_4.sce
new file mode 100755
index 000000000..9876ff179
--- /dev/null
+++ b/2453/CH9/EX9.4/9_4.sce
@@ -0,0 +1,9 @@
+//To calculate the critical temperature
+T1 = 3; //temperature, K
+T2 = 7.1; //temperature, K
+lamda_T1 = 396; //penetration depth, armstrong
+lamda_T2 = 1730; //penetration depth, armstrong
+A = (((lamda_T2/lamda_T1)^2)*T2^4) - T1^4;
+B = ((lamda_T2/lamda_T1)^2)-1;
+Tc = (A/B)^(1/4); //critical temperature, K
+printf("critical temperature for lead is %5.3f K",Tc);