1 files changed, 20 insertions, 0 deletions
diff --git a/812/CH4/EX4.02/4_02.sce b/812/CH4/EX4.02/4_02.sce
new file mode 100755
index 000000000..a209b74e0
--- /dev/null
+++ b/812/CH4/EX4.02/4_02.sce
@@ -0,0 +1,20 @@
+//Mass flow//
+pathname=get_absolute_file_path('4.02.sce')
+filename=pathname+filesep()+'4.02-data.sci'
+exec(filename)
+//If I=integral of(pV.dA):
+//For system: ICS=Iab+Ibc+Icd+Ida
+//But ICS=0
+
+//For Aab:
+function p=f(y),p=-d*U*w*y^0,endfunction
+IAab=intg(0,t,f)
+
+//For Acd:
+function q=g(y),q=d*U*w*(2*y/t-(y/t)^2),endfunction
+IAcd=intg(0,t,g)
+
+//Mass flow rate across surface bc(in kg/sec):
+mbc=(-IAab-IAcd)/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate across surface bc: %.4f kg/sec\n\n",mbc)