1 files changed, 26 insertions, 0 deletions
diff --git a/147/CH3/EX3.3/Example3_3.sce b/147/CH3/EX3.3/Example3_3.sce
new file mode 100644
index 000000000..c81b4d6c2
--- /dev/null
+++ b/147/CH3/EX3.3/Example3_3.sce
@@ -0,0 +1,26 @@
+close();
+clear;
+clc;
+//resistance 'R', rate of electrical energy 'rate', time of operation 't', time period 'T'
+T = 0.01; //s
+R = 20; //ohm
+rate = 0.06; //$ per kWh
+time = 24; //hours
+//from graph, voltage as function of time 'v'
+function V = v(t)
+  V = 10^4 * t; //V
+  endfunction
+//current as a function of time 'i'
+function I = i(t)
+  I = v(t)/R; //A
+endfunction
+//power as afunction of time 'p'
+function power = p(t)
+  power = v(t)*i(t); //W
+  endfunction
+//power 'P'
+P = 1/T * intg(0,T,p)/1000; //kW
+energy = P*time; //kWh
+//cost of electrical energy 'cost'
+cost = rate*energy; //$
+mprintf("Electrical energy costs $ %0.2f for %d hours",cost,time);