summaryrefslogtreecommitdiff
path: root/3014/CH8/EX8.9/Ex8_9.sce
diff options
context:
space:
mode:
Diffstat (limited to '3014/CH8/EX8.9/Ex8_9.sce')
-rwxr-xr-x3014/CH8/EX8.9/Ex8_9.sce18
1 files changed, 18 insertions, 0 deletions
diff --git a/3014/CH8/EX8.9/Ex8_9.sce b/3014/CH8/EX8.9/Ex8_9.sce
new file mode 100755
index 000000000..4415a0c05
--- /dev/null
+++ b/3014/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,18 @@
+
+ clc
+// Given that
+m_0 = 9.1e-31 // Mass of electron in kg
+mu_0 = 1.256e-6// SI
+e = 1.6e-19 // Charge on electron in coulomb
+eta_s = 1e28 // superelectron density in no. per cube
+T_1 = 0 // First temp in kelvin
+T_2 = 1 // Second temp in kelvin
+T_c = 3 // Critical temp in kelvin
+
+printf("Example 8.9\n")
+printf("Standard formula used \tlambda_0 = sqrt(m_0/(mu_0*eta_s*e^2))\n")
+lambda_0 = sqrt(m_0/(mu_0*eta_s*e^2))// Calculation of penetration depth at 0K
+lambda_t = lambda_0/sqrt(1-(T_2/T_c)^4) // Calculation of penetration depth at 2K
+
+printf("Penetration depth at %d K is %d angestrom.",T_1,lambda_0*1e10)
+printf("\nPenetration depth at %d K is %f angestrom.\n\n\n",T_2,lambda_t*1e10)