1 files changed, 37 insertions, 0 deletions
diff --git a/3012/CH2/EX2.1/Ex2_1.sce b/3012/CH2/EX2.1/Ex2_1.sce
new file mode 100755
index 000000000..40bfbcdb8
--- /dev/null
+++ b/3012/CH2/EX2.1/Ex2_1.sce
@@ -0,0 +1,37 @@
+// Given:-
+p1 = 3*(10**5)                                // initial pressure of gas in pascal
+v1 = 0.1                                      // initial volumme of gas in meter^3
+v2 = 0.2                                      // final volume of gas in meter^3
+
+// calculations
+// Part (a) i.e. n=1.5
+//constant = p1*(v1**n)                         // p*(v^n) = constant
+constant1 = p1*(v1**1.5) 
+constant2 = p1*(v1**1) 
+constant3 = p1*(v1**0) 
+// function p                     
+function v = p1(v)
+    v = constant1/(v^1.5)
+endfunction
+
+function v = p2(v)
+    v = constant2/(v^1)
+endfunction
+
+function v = p3(v)
+    v = constant3/(v^0)
+endfunction
+
+work1 = intg(v1,v2,p1)                   // integrating pdv from initial to final volume 
+w1 = work1(1)/1000                             // divided by 1000 to convert to KJ
+printf( 'The work done for n=1.5 in KJ is %.2f',w1)
+
+//part(b) i.e. n = 1
+work2 = intg(v1,v2,p2)
+w2 = work2(1)/1000
+printf( 'The work done for n=1 in KJ is %.2f',w2)
+
+//part(c) i.e. n=0
+work3 = intg(v1,v2,p3)
+w3 = work3(1)/1000
+printf( 'The work done for n=0 in KJ is %.2f',w3)