summaryrefslogtreecommitdiff
path: root/728
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /728
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '728')
-rwxr-xr-x728/CH3/EX3.1/EX3_1.sce10
-rwxr-xr-x728/CH3/EX3.1/Ex3_1.txt11
-rwxr-xr-x728/CH3/EX3.1/Ex3_1_ans.txt3
-rwxr-xr-x728/CH3/EX3.10/Ex3_10.sce22
-rwxr-xr-x728/CH3/EX3.10/Ex3_10.txt26
-rwxr-xr-x728/CH3/EX3.10/Ex3_10_ans.txt12
-rwxr-xr-x728/CH3/EX3.2/Ex3_2.sce31
-rwxr-xr-x728/CH3/EX3.2/Ex3_2.txt34
-rwxr-xr-x728/CH3/EX3.2/Ex3_2_ans.txt16
-rwxr-xr-x728/CH3/EX3.3/Ex3_3.sce9
-rwxr-xr-x728/CH3/EX3.3/Ex3_3.txt13
-rwxr-xr-x728/CH3/EX3.3/Ex3_3_ans.txt4
-rwxr-xr-x728/CH3/EX3.4/Ex3_4.sce22
-rwxr-xr-x728/CH3/EX3.4/Ex3_4.txt26
-rwxr-xr-x728/CH3/EX3.4/Ex3_4_ans.txt12
-rwxr-xr-x728/CH3/EX3.5/Ex3_5.sce17
-rwxr-xr-x728/CH3/EX3.5/Ex3_5.txt20
-rwxr-xr-x728/CH3/EX3.5/Ex3_5_ans.txt16
-rwxr-xr-x728/CH3/EX3.6/Ex3_6.sce15
-rwxr-xr-x728/CH3/EX3.6/Ex3_6.txt19
-rwxr-xr-x728/CH3/EX3.6/Ex3_6_ans.txt8
-rwxr-xr-x728/CH3/EX3.7/Ex3_7.sce15
-rwxr-xr-x728/CH3/EX3.7/Ex3_7.txt18
-rwxr-xr-x728/CH3/EX3.7/Ex3_7_ans.txt10
-rwxr-xr-x728/CH3/EX3.8/Ex3_8.sce24
-rwxr-xr-x728/CH3/EX3.8/Ex3_8.txt25
-rwxr-xr-x728/CH3/EX3.8/Ex3_8_ans.txt28
-rwxr-xr-x728/CH3/EX3.9/Ex3_9.sce18
-rwxr-xr-x728/CH3/EX3.9/Ex3_9.txt22
-rwxr-xr-x728/CH3/EX3.9/Ex3_9_ans.txt12
-rwxr-xr-x728/CH4/EX4.1/Ex4_1.sce18
-rwxr-xr-x728/CH4/EX4.1/Ex4_1.txt18
-rwxr-xr-x728/CH4/EX4.1/Ex4_1_ans.txt15
-rwxr-xr-x728/CH4/EX4.10/Ex4_10.sce15
-rwxr-xr-x728/CH4/EX4.10/Ex4_10.txt15
-rwxr-xr-x728/CH4/EX4.10/Ex4_10_ans.txt13
-rwxr-xr-x728/CH4/EX4.11/Ex4_11.sce29
-rwxr-xr-x728/CH4/EX4.11/Ex4_11.txt29
-rwxr-xr-x728/CH4/EX4.11/Ex4_11_ans.txt6
-rwxr-xr-x728/CH4/EX4.12/Ex4_11.sce29
-rwxr-xr-x728/CH4/EX4.12/Ex4_11.txt29
-rwxr-xr-x728/CH4/EX4.12/Ex4_11_ans.txt6
-rwxr-xr-x728/CH4/EX4.13/Ex4_12.sce16
-rwxr-xr-x728/CH4/EX4.13/Ex4_12.txt16
-rwxr-xr-x728/CH4/EX4.13/Ex4_12_ans.txt12
-rwxr-xr-x728/CH4/EX4.14/Ex4_14.sce22
-rwxr-xr-x728/CH4/EX4.14/Ex4_14.txt22
-rwxr-xr-x728/CH4/EX4.14/Ex4_14_ans.txt12
-rwxr-xr-x728/CH4/EX4.15/Ex4_15.sce45
-rwxr-xr-x728/CH4/EX4.15/Ex4_15.txt45
-rwxr-xr-x728/CH4/EX4.15/Ex4_15_ans.txt20
-rwxr-xr-x728/CH4/EX4.16/Ex4_16.sce16
-rwxr-xr-x728/CH4/EX4.16/Ex4_16.txt16
-rwxr-xr-x728/CH4/EX4.16/Ex4_16_ans.txt4
-rwxr-xr-x728/CH4/EX4.17/EX4_17_ans.txt8
-rwxr-xr-x728/CH4/EX4.17/Ex4_17.sce14
-rwxr-xr-x728/CH4/EX4.17/Ex4_17.txt14
-rwxr-xr-x728/CH4/EX4.18/Ex4_18.sce20
-rwxr-xr-x728/CH4/EX4.18/Ex4_18.txt20
-rwxr-xr-x728/CH4/EX4.18/Ex4_18_ans.txt2
-rwxr-xr-x728/CH4/EX4.19/Ex4_19.sce21
-rwxr-xr-x728/CH4/EX4.19/Ex4_19.txt21
-rwxr-xr-x728/CH4/EX4.19/Ex4_19_ans.txt5
-rwxr-xr-x728/CH4/EX4.2/Ex4_2.sce25
-rwxr-xr-x728/CH4/EX4.2/Ex4_2.txt25
-rwxr-xr-x728/CH4/EX4.2/Ex4_2_ans.txt19
-rwxr-xr-x728/CH4/EX4.20/Ex4_20.sce15
-rwxr-xr-x728/CH4/EX4.20/Ex4_20.txt15
-rwxr-xr-x728/CH4/EX4.20/Ex4_20_ans.txt4
-rwxr-xr-x728/CH4/EX4.21/Ex4_21.sce15
-rwxr-xr-x728/CH4/EX4.21/Ex4_21.txt15
-rwxr-xr-x728/CH4/EX4.21/Ex4_21_ans.txt5
-rwxr-xr-x728/CH4/EX4.22/Ex4_22.sce15
-rwxr-xr-x728/CH4/EX4.22/Ex4_22.txt15
-rwxr-xr-x728/CH4/EX4.22/Ex4_22_ans.txt5
-rwxr-xr-x728/CH4/EX4.23/Ex4_23.sce13
-rwxr-xr-x728/CH4/EX4.23/Ex4_23.txt13
-rwxr-xr-x728/CH4/EX4.23/Ex4_23_ans.txt4
-rwxr-xr-x728/CH4/EX4.24/Ex4_24.sce15
-rwxr-xr-x728/CH4/EX4.24/Ex4_24.txt15
-rwxr-xr-x728/CH4/EX4.24/Ex4_24_ans.txt5
-rwxr-xr-x728/CH4/EX4.3/EX4_3_ans.txt3
-rwxr-xr-x728/CH4/EX4.3/Ex4_3.sce12
-rwxr-xr-x728/CH4/EX4.3/Ex4_3.txt12
-rwxr-xr-x728/CH4/EX4.4/Ex4_4.sce13
-rwxr-xr-x728/CH4/EX4.4/Ex4_4.txt13
-rwxr-xr-x728/CH4/EX4.4/Ex4_4_ans.txt9
-rwxr-xr-x728/CH4/EX4.5/Ex4_5.sce25
-rwxr-xr-x728/CH4/EX4.5/Ex4_5.txt25
-rwxr-xr-x728/CH4/EX4.5/Ex4_5_ans.txt28
-rwxr-xr-x728/CH4/EX4.6/Ex4_6.sce24
-rwxr-xr-x728/CH4/EX4.6/Ex4_6.txt24
-rwxr-xr-x728/CH4/EX4.6/Ex4_6_ans.txt8
-rwxr-xr-x728/CH4/EX4.7/Ex4_7.sce12
-rwxr-xr-x728/CH4/EX4.7/Ex4_7.txt12
-rwxr-xr-x728/CH4/EX4.7/Ex4_7_ans.txt4
-rwxr-xr-x728/CH4/EX4.8/Ex4_8.sce17
-rwxr-xr-x728/CH4/EX4.8/Ex4_8.txt17
-rwxr-xr-x728/CH4/EX4.8/Ex4_8_ans.txt15
-rwxr-xr-x728/CH4/EX4.9/EX4_9_ans.txt15
-rwxr-xr-x728/CH4/EX4.9/Ex4_9.sce24
-rwxr-xr-x728/CH4/EX4.9/Ex4_9.txt24
-rwxr-xr-x728/CH5/EX5.1/Ex5_1.sce12
-rwxr-xr-x728/CH5/EX5.1/Ex5_1.txt12
-rwxr-xr-x728/CH5/EX5.1/Ex5_1_ans.txt3
-rwxr-xr-x728/CH5/EX5.2/EX5_2.sce16
-rwxr-xr-x728/CH5/EX5.2/Ex5_2.txt16
-rwxr-xr-x728/CH5/EX5.2/Ex5_2_ans.txt6
-rwxr-xr-x728/CH5/EX5.3/Ex5_3.sce19
-rwxr-xr-x728/CH5/EX5.3/Ex5_3.txt19
-rwxr-xr-x728/CH5/EX5.3/Ex5_3_ans.txt4
-rwxr-xr-x728/CH5/EX5.4/Ex5_4.sce16
-rwxr-xr-x728/CH5/EX5.4/Ex5_4.txt16
-rwxr-xr-x728/CH5/EX5.4/Ex5_4_ans.txt4
-rwxr-xr-x728/CH6/EX6.10/Ex6_10.sce17
-rwxr-xr-x728/CH6/EX6.10/Ex6_10.txt17
-rwxr-xr-x728/CH6/EX6.10/Ex6_10_ans.txt5
-rwxr-xr-x728/CH6/EX6.11/Ex6_11.sce13
-rwxr-xr-x728/CH6/EX6.11/Ex6_11.txt13
-rwxr-xr-x728/CH6/EX6.11/Ex6_11_ans.txt12
-rwxr-xr-x728/CH6/EX6.12/Ex6_12.sce10
-rwxr-xr-x728/CH6/EX6.12/Ex6_12.txt10
-rwxr-xr-x728/CH6/EX6.12/Ex6_12_ans.txt5
-rwxr-xr-x728/CH6/EX6.13/Ex6_13.sce9
-rwxr-xr-x728/CH6/EX6.13/Ex6_13.txt9
-rwxr-xr-x728/CH6/EX6.13/Ex6_13_ans.txt5
-rwxr-xr-x728/CH6/EX6.2/Ex6_2.sce13
-rwxr-xr-x728/CH6/EX6.2/Ex6_2.txt13
-rwxr-xr-x728/CH6/EX6.2/Ex6_2_ans.txt5
-rwxr-xr-x728/CH6/EX6.3/Ex6_3.sce28
-rwxr-xr-x728/CH6/EX6.3/Ex6_3.txt28
-rwxr-xr-x728/CH6/EX6.3/Ex6_3_ans.txt8
-rwxr-xr-x728/CH6/EX6.4/Ex6_4.sce14
-rwxr-xr-x728/CH6/EX6.4/Ex6_4.txt14
-rwxr-xr-x728/CH6/EX6.4/Ex6_4_ans.txt13
-rwxr-xr-x728/CH6/EX6.5/Ex6_5.sce13
-rwxr-xr-x728/CH6/EX6.5/Ex6_5.txt13
-rwxr-xr-x728/CH6/EX6.5/Ex6_5_ans.txt9
-rwxr-xr-x728/CH6/EX6.6/Ex6_6.sce20
-rwxr-xr-x728/CH6/EX6.6/Ex6_6.txt20
-rwxr-xr-x728/CH6/EX6.6/Ex6_6_ans.txt16
-rwxr-xr-x728/CH6/EX6.7/Ex6_7.sce13
-rwxr-xr-x728/CH6/EX6.7/Ex6_7.txt13
-rwxr-xr-x728/CH6/EX6.7/Ex6_7_ans.txt6
-rwxr-xr-x728/CH6/EX6.9/Ex6_9.sce20
-rwxr-xr-x728/CH6/EX6.9/Ex6_9.txt20
-rwxr-xr-x728/CH6/EX6.9/Ex6_9_ans.txt7
-rwxr-xr-x728/CH7/EX7.1/Ex7_1.sce17
-rwxr-xr-x728/CH7/EX7.1/Ex7_1.txt17
-rwxr-xr-x728/CH7/EX7.1/Ex7_1_ans.txt5
-rwxr-xr-x728/CH7/EX7.2/Ex7_2.sce11
-rwxr-xr-x728/CH7/EX7.2/Ex7_2.txt11
-rwxr-xr-x728/CH7/EX7.2/Ex7_2_ans.txt5
-rwxr-xr-x728/CH7/EX7.3/Ex7_3.sce11
-rwxr-xr-x728/CH7/EX7.3/Ex7_3.txt11
-rwxr-xr-x728/CH7/EX7.3/Ex7_3_ans.txt4
-rwxr-xr-x728/CH7/EX7.4/Ex7_4.sce13
-rwxr-xr-x728/CH7/EX7.4/Ex7_4.txt13
-rwxr-xr-x728/CH7/EX7.4/Ex7_4_ans.txt5
-rwxr-xr-x728/CH8/EX8.1/Ex8_1.sce26
-rwxr-xr-x728/CH8/EX8.1/Ex8_1.txt26
-rwxr-xr-x728/CH8/EX8.1/Ex8_1_ans.txt24
-rwxr-xr-x728/CH8/EX8.10/Ex8_10.sce23
-rwxr-xr-x728/CH8/EX8.10/Ex8_10.txt23
-rwxr-xr-x728/CH8/EX8.10/Ex8_10_ans.txt16
-rwxr-xr-x728/CH8/EX8.11/Ex8_11.sce24
-rwxr-xr-x728/CH8/EX8.11/Ex8_11.txt24
-rwxr-xr-x728/CH8/EX8.11/Ex8_11_ans.txt24
-rwxr-xr-x728/CH8/EX8.12/Ex8_12.sce12
-rwxr-xr-x728/CH8/EX8.12/Ex8_12.txt12
-rwxr-xr-x728/CH8/EX8.12/Ex8_12_ans.txt4
-rwxr-xr-x728/CH8/EX8.13/Ex8_13.sce29
-rwxr-xr-x728/CH8/EX8.13/Ex8_13.txt29
-rwxr-xr-x728/CH8/EX8.13/Ex8_13_ans.txt12
-rwxr-xr-x728/CH8/EX8.14/Ex8_14.sce17
-rwxr-xr-x728/CH8/EX8.14/Ex8_14.txt17
-rwxr-xr-x728/CH8/EX8.14/Ex8_14_ans.txt12
-rwxr-xr-x728/CH8/EX8.15/Ex8_15.sce16
-rwxr-xr-x728/CH8/EX8.15/Ex8_15.txt16
-rwxr-xr-x728/CH8/EX8.15/Ex8_15_ans.txt12
-rwxr-xr-x728/CH8/EX8.16/Ex8_16.sce20
-rwxr-xr-x728/CH8/EX8.16/Ex8_16.txt20
-rwxr-xr-x728/CH8/EX8.16/Ex8_16_ans.txt8
-rwxr-xr-x728/CH8/EX8.2/Ex8_2.sce16
-rwxr-xr-x728/CH8/EX8.2/Ex8_2.txt16
-rwxr-xr-x728/CH8/EX8.2/Ex8_2_ans.txt12
-rwxr-xr-x728/CH8/EX8.3/Ex8_3.sce16
-rwxr-xr-x728/CH8/EX8.3/Ex8_3.txt16
-rwxr-xr-x728/CH8/EX8.3/Ex8_3_ans.txt11
-rwxr-xr-x728/CH8/EX8.4/Ex8_4.sce23
-rwxr-xr-x728/CH8/EX8.4/Ex8_4.txt23
-rwxr-xr-x728/CH8/EX8.4/Ex8_4_ans.txt16
-rwxr-xr-x728/CH8/EX8.5/Ex8_5.sce28
-rwxr-xr-x728/CH8/EX8.5/Ex8_5.txt28
-rwxr-xr-x728/CH8/EX8.5/Ex8_5_ans.txt16
-rwxr-xr-x728/CH8/EX8.6/Ex8_6.sce26
-rwxr-xr-x728/CH8/EX8.6/Ex8_6.txt26
-rwxr-xr-x728/CH8/EX8.6/Ex8_6_ans.txt12
-rwxr-xr-x728/CH8/EX8.7/Ex8_7.sce16
-rwxr-xr-x728/CH8/EX8.7/Ex8_7.txt16
-rwxr-xr-x728/CH8/EX8.7/Ex8_7_ans.txt12
-rwxr-xr-x728/CH8/EX8.8/Ex8_8.sce20
-rwxr-xr-x728/CH8/EX8.8/Ex8_8.txt20
-rwxr-xr-x728/CH8/EX8.8/Ex8_8_ans.txt13
-rwxr-xr-x728/CH8/EX8.9/EX8_9_ans.txt8
-rwxr-xr-x728/CH8/EX8.9/Ex8_9.sce16
-rwxr-xr-x728/CH8/EX8.9/Ex8_9.txt16
-rwxr-xr-x728/CH9/EX9.1/Ex9_1.sce9
-rwxr-xr-x728/CH9/EX9.1/Ex9_1.txt9
-rwxr-xr-x728/CH9/EX9.1/Ex9_1_ans.txt5
-rwxr-xr-x728/CH9/EX9.10/Ex9_10.sce11
-rwxr-xr-x728/CH9/EX9.10/Ex9_10.txt11
-rwxr-xr-x728/CH9/EX9.10/Ex9_10_ans.txt9
-rwxr-xr-x728/CH9/EX9.11/Ex9_11.sce14
-rwxr-xr-x728/CH9/EX9.11/Ex9_11.txt14
-rwxr-xr-x728/CH9/EX9.11/Ex9_11_ans.txt8
-rwxr-xr-x728/CH9/EX9.12/Ex9_12.sce18
-rwxr-xr-x728/CH9/EX9.12/Ex9_12.txt18
-rwxr-xr-x728/CH9/EX9.12/Ex9_12_ans.txt11
-rwxr-xr-x728/CH9/EX9.13/Ex9_13.sce32
-rwxr-xr-x728/CH9/EX9.13/Ex9_13.txt32
-rwxr-xr-x728/CH9/EX9.13/Ex9_13_ans.txt16
-rwxr-xr-x728/CH9/EX9.2/Ex9_2.sce10
-rwxr-xr-x728/CH9/EX9.2/Ex9_2.txt10
-rwxr-xr-x728/CH9/EX9.2/Ex9_2_ans.txt4
-rwxr-xr-x728/CH9/EX9.3/Ex9_3.sce13
-rwxr-xr-x728/CH9/EX9.3/Ex9_3.txt13
-rwxr-xr-x728/CH9/EX9.3/Ex9_3_ans.txt8
-rwxr-xr-x728/CH9/EX9.4/Ex9_4.sce17
-rwxr-xr-x728/CH9/EX9.4/Ex9_4.txt17
-rwxr-xr-x728/CH9/EX9.4/Ex9_4_ans.txt11
-rwxr-xr-x728/CH9/EX9.5/Ex9_5.sce10
-rwxr-xr-x728/CH9/EX9.5/Ex9_5.txt10
-rwxr-xr-x728/CH9/EX9.5/Ex9_5_ans.txt4
-rwxr-xr-x728/CH9/EX9.6/Ex9_6.sce9
-rwxr-xr-x728/CH9/EX9.6/Ex9_6.txt9
-rwxr-xr-x728/CH9/EX9.6/Ex9_6_ans.txt5
-rwxr-xr-x728/CH9/EX9.7/Ex9_7.sce9
-rwxr-xr-x728/CH9/EX9.7/Ex9_7.txt9
-rwxr-xr-x728/CH9/EX9.7/Ex9_7_ans.txt4
-rwxr-xr-x728/CH9/EX9.8/EX9_8.txt12
-rwxr-xr-x728/CH9/EX9.8/Ex9_8.sce12
-rwxr-xr-x728/CH9/EX9.8/Ex9_8_ans.txt9
-rwxr-xr-x728/CH9/EX9.9/Ex9_9.sce15
-rwxr-xr-x728/CH9/EX9.9/Ex9_9.txt15
-rwxr-xr-x728/CH9/EX9.9/Ex9_9_ans.txt8
246 files changed, 3710 insertions, 0 deletions
diff --git a/728/CH3/EX3.1/EX3_1.sce b/728/CH3/EX3.1/EX3_1.sce
new file mode 100755
index 000000000..362dd1ce9
--- /dev/null
+++ b/728/CH3/EX3.1/EX3_1.sce
@@ -0,0 +1,10 @@
+//Caption:Program to find value of terminating impedance of lossless transmission line.
+//Exa:3.1
+clc;
+clear;
+close;
+//Given:
+Z_ch=100;//in ohms
+S=5;//VSWR (unitless)
+Z=Z_ch*S;
+printf("\n\n\t The terminating impedence = %d ohms ",Z); \ No newline at end of file
diff --git a/728/CH3/EX3.1/Ex3_1.txt b/728/CH3/EX3.1/Ex3_1.txt
new file mode 100755
index 000000000..f2235deb1
--- /dev/null
+++ b/728/CH3/EX3.1/Ex3_1.txt
@@ -0,0 +1,11 @@
+//Caption:Program to find value of terminating impedance of lossless transmission line.
+//Exa:3.1
+clc;
+clear;
+close;
+//Given:
+Z_ch=100;//in ohms
+S=5;//VSWR (unitless)
+Z=Z_ch*S;
+printf("\n\n\t The terminating impedence = %d ohms ",Z);
+
diff --git a/728/CH3/EX3.1/Ex3_1_ans.txt b/728/CH3/EX3.1/Ex3_1_ans.txt
new file mode 100755
index 000000000..c1e309d8e
--- /dev/null
+++ b/728/CH3/EX3.1/Ex3_1_ans.txt
@@ -0,0 +1,3 @@
+
+ The terminating impedence = 500 ohms
+
diff --git a/728/CH3/EX3.10/Ex3_10.sce b/728/CH3/EX3.10/Ex3_10.sce
new file mode 100755
index 000000000..6c0280f50
--- /dev/null
+++ b/728/CH3/EX3.10/Ex3_10.sce
@@ -0,0 +1,22 @@
+//Caption:Calculate the charcterstic impedance & phase velocity
+//Exa:3.10
+clc;
+clear;
+close;
+e=2.718;
+f=1000;//in Hz
+l=10000;//in meters
+Z_sc=(2631+%i*1289);//in ohms
+Z_oc=(221-%i*137);//in ohms
+Z_o=[Z_sc*Z_oc]^0.5;
+Z_mag=[real(Z_o)^2+imag(Z_o)^2]^0.5;
+Z_ang=[atan((imag(Z_o))/real(Z_o))]*180/%pi;
+disp(Z_mag,"Characteristic impedance (in ohms) =");
+disp(Z_ang,"Angle (in degrees) =");
+x=[(Z_oc/Z_sc)^0.5];
+//x=tanh(v*l)
+//As, tanh(t)=[e^t-e^-t]/[e^t+e^-t]
+v=(261+%i*2988)/l;
+a=real(v);
+b=imag(v);
+disp(2*%pi*f/b,"Phase velocity (in meter per sec.) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.10/Ex3_10.txt b/728/CH3/EX3.10/Ex3_10.txt
new file mode 100755
index 000000000..eb7d10e35
--- /dev/null
+++ b/728/CH3/EX3.10/Ex3_10.txt
@@ -0,0 +1,26 @@
+//Caption:Calculate the charcterstic impedance & phase velocity
+//Exa:3.10
+clc;
+clear;
+close;
+e=2.718;
+f=1000;//in Hz
+l=10000;//in meters
+Z_sc=(2631+%i*1289);//in ohms
+Z_oc=(221-%i*137);//in ohms
+Z_o=[Z_sc*Z_oc]^0.5;
+Z_mag=[real(Z_o)^2+imag(Z_o)^2]^0.5;
+Z_ang=[atan((imag(Z_o))/real(Z_o))]*180/%pi;
+disp(Z_mag,"Characteristic impedance (in ohms) =");
+disp(Z_ang,"Angle (in degrees) =");
+x=[(Z_oc/Z_sc)^0.5];
+//x=tanh(v*l)
+//As, tanh(t)=[e^t-e^-t]/[e^t+e^-t]
+v=(261+%i*2988)/l;
+a=real(v);
+b=imag(v);
+disp(2*%pi*f/b,"Phase velocity (in meter per sec.) =");
+
+
+
+
diff --git a/728/CH3/EX3.10/Ex3_10_ans.txt b/728/CH3/EX3.10/Ex3_10_ans.txt
new file mode 100755
index 000000000..e19688f2b
--- /dev/null
+++ b/728/CH3/EX3.10/Ex3_10_ans.txt
@@ -0,0 +1,12 @@
+Characteristic impedance (in ohms) =
+
+ 872.81286
+
+ Angle (in degrees) =
+
+ - 2.846825
+
+ Phase velocity (in meter per sec.) =
+
+ 21028.063
+
diff --git a/728/CH3/EX3.2/Ex3_2.sce b/728/CH3/EX3.2/Ex3_2.sce
new file mode 100755
index 000000000..aa91a23f5
--- /dev/null
+++ b/728/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,31 @@
+//Caption:Calculate the charcteristic impedance, attenuation constant, phase constant of transmission line Calculate power delivered to load , if line length=500 km.
+//Exa:3.2
+clc;
+clear;
+close;
+e=2.718;
+//Given:
+R=8;//in ohm/kilometer
+L=2*10^-3;//in henry/kilometer
+C=0.002*10^-6;//in farad/kilometer
+G=0.07*10^-6;//second/kilometer
+f=2000;//in hertz
+//Since [w=2*(pi)*f] & [Zch={(R+jwL)/(G+jwC)}^0.5]
+w=2*%pi*f;//in radians
+Z_ch={(R+%i*w*L)/(G+%i*w*C)}^0.5;//computing characteristic impedance
+disp(Z_ch,"Characteristic impedance (in ohms) =");
+y={(R+%i*w*L)*(G+%i*w*C)}^0.5;
+a=real(y);//atteneuation constant
+b=imag(y);//phase constant
+disp(a,"Atteneuation constant (in NP/km) =");
+disp(b,"Phase constant (in radian/km) =");
+V_in=2;//in volts
+l=500;//in kilometers
+Z_in=Z_ch;//Since line terminated at its char. imped. so, Z_in=Z_ch=Z(load)
+I_s=V_in/Z_in;
+Imag=[[{{real(I_s)}^2}+{{imag(I_s)}^2}]^0.5]*10^3;//in milliampere
+Iang=atan(imag(I_s)/real(I_s))*(180/%pi);//in degrees
+I=Imag*e^-1.99;//I=Is*e^-yl
+//P(power delivered)=I*I*REAL(Z_ch)
+P=I*I*real(Z_ch);
+disp(P,"Power delivered to load (in microwatt =)"); \ No newline at end of file
diff --git a/728/CH3/EX3.2/Ex3_2.txt b/728/CH3/EX3.2/Ex3_2.txt
new file mode 100755
index 000000000..d08f8b69f
--- /dev/null
+++ b/728/CH3/EX3.2/Ex3_2.txt
@@ -0,0 +1,34 @@
+//Caption:Calculate the charcteristic impedance, attenuation constant, phase constant of transmission line Calculate power delivered to load , if line length=500 km.
+//Exa:3.2
+clc;
+clear;
+close;
+e=2.718;
+//Given:
+R=8;//in ohm/kilometer
+L=2*10^-3;//in henry/kilometer
+C=0.002*10^-6;//in farad/kilometer
+G=0.07*10^-6;//second/kilometer
+f=2000;//in hertz
+//Since [w=2*(pi)*f] & [Zch={(R+jwL)/(G+jwC)}^0.5]
+w=2*%pi*f;//in radians
+Z_ch={(R+%i*w*L)/(G+%i*w*C)}^0.5;//computing characteristic impedance
+disp(Z_ch,"Characteristic impedance (in ohms) =");
+y={(R+%i*w*L)*(G+%i*w*C)}^0.5;
+a=real(y);//atteneuation constant
+b=imag(y);//phase constant
+disp(a,"Atteneuation constant (in NP/km) =");
+disp(b,"Phase constant (in radian/km) =");
+V_in=2;//in volts
+l=500;//in kilometers
+Z_in=Z_ch;//Since line terminated at its char. imped. so, Z_in=Z_ch=Z(load)
+I_s=V_in/Z_in;
+Imag=[[{{real(I_s)}^2}+{{imag(I_s)}^2}]^0.5]*10^3;//in milliampere
+Iang=atan(imag(I_s)/real(I_s))*(180/%pi);//in degrees
+I=Imag*e^-1.99;//I=Is*e^-yl
+//P(power delivered)=I*I*REAL(Z_ch)
+P=I*I*real(Z_ch);
+disp(P,"Power delivered to load (in microwatt =)");
+
+
+
diff --git a/728/CH3/EX3.2/Ex3_2_ans.txt b/728/CH3/EX3.2/Ex3_2_ans.txt
new file mode 100755
index 000000000..6ddfa8e0d
--- /dev/null
+++ b/728/CH3/EX3.2/Ex3_2_ans.txt
@@ -0,0 +1,16 @@
+ Characteristic impedance (in ohms) =
+
+ 1012.5002 - 155.81342i
+
+ Atteneuation constant (in NP/km) =
+
+ 0.0039869
+
+ Phase constant (in radian/km) =
+
+ 0.0254360
+
+ Power delivered to load (in microwatt =)
+
+ 72.141799
+
diff --git a/728/CH3/EX3.3/Ex3_3.sce b/728/CH3/EX3.3/Ex3_3.sce
new file mode 100755
index 000000000..83bf5687e
--- /dev/null
+++ b/728/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,9 @@
+//Caption:Calculate phase velocity of the wave that propogates on line as given in example 3.2
+//Exa:3.3
+clc;
+clear;
+close;
+w=4*%pi*10^3;//in rad/sec
+b=0.02543;//in rad/km
+V_p=w/b;// phase velocity
+disp(V_p,"Phase velocity (in km/sec) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.3/Ex3_3.txt b/728/CH3/EX3.3/Ex3_3.txt
new file mode 100755
index 000000000..654b86c01
--- /dev/null
+++ b/728/CH3/EX3.3/Ex3_3.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate phase velocity of the wave that propogates on line as given in example 3.2
+//Exa:3.3
+clc;
+clear;
+close;
+w=4*%pi*10^3;//in rad/sec
+b=0.02543;//in rad/km
+V_p=w/b;// phase velocity
+disp(V_p,"Phase velocity (in km/sec) =");// But answer in book is misprinted as 494.22
+
+
+
+
diff --git a/728/CH3/EX3.3/Ex3_3_ans.txt b/728/CH3/EX3.3/Ex3_3_ans.txt
new file mode 100755
index 000000000..2b89db88f
--- /dev/null
+++ b/728/CH3/EX3.3/Ex3_3_ans.txt
@@ -0,0 +1,4 @@
+Phase velocity (in km/sec) =
+
+ 494155.35
+
diff --git a/728/CH3/EX3.4/Ex3_4.sce b/728/CH3/EX3.4/Ex3_4.sce
new file mode 100755
index 000000000..e606c7846
--- /dev/null
+++ b/728/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,22 @@
+//Caption:Calculate (a)-Current drawn from generator. (b)-Magnitune & phase of load current. (c)-Power delivered to load.
+//Exa: 3.4
+clc;
+clear;
+close;
+f=37.5*10^6;//frequency(in hertz)
+wl=(3*10^8)/f;//wavelength (in meters)
+Z_l=100;//in ohms
+Z_o=200;//in ohms
+l=5*wl/4;//length of line (in meters)
+b=2*%pi/wl;
+//At generator end,
+Z_i=Z_o*(Z_l+%i*Z_o*tan(b*l))/(Z_o+%i*Z_l*tan(b*l));
+V_s=200*Z_i/(200+Z_i);
+I_s=200/(200+Z_i);
+disp(real(I_s),"Current drawn from generator(in amps) =");
+//for a lossless line , P(avg)*I_input=P(avg)*I_load
+P_avg=V_s*I_s;//in watts
+disp(real(P_avg),"Power delivered to load (in watts) =)");
+//Real(Vs*Is)=Real(Vs*I_load)
+I_load=(P_avg/Z_l)^0.5;//in amps
+disp(real(I_load),"Current flowing in load (in amps) =)"); \ No newline at end of file
diff --git a/728/CH3/EX3.4/Ex3_4.txt b/728/CH3/EX3.4/Ex3_4.txt
new file mode 100755
index 000000000..1b320a4bb
--- /dev/null
+++ b/728/CH3/EX3.4/Ex3_4.txt
@@ -0,0 +1,26 @@
+//Caption:Calculate (a)-Current drawn from generator. (b)-Magnitune & phase of load current. (c)-Power delivered to load.
+//Exa: 3.4
+clc;
+clear;
+close;
+f=37.5*10^6;//frequency(in hertz)
+wl=(3*10^8)/f;//wavelength (in meters)
+Z_l=100;//in ohms
+Z_o=200;//in ohms
+l=5*wl/4;//length of line (in meters)
+b=2*%pi/wl;
+//At generator end,
+Z_i=Z_o*(Z_l+%i*Z_o*tan(b*l))/(Z_o+%i*Z_l*tan(b*l));
+V_s=200*Z_i/(200+Z_i);
+I_s=200/(200+Z_i);
+disp(real(I_s),"Current drawn from generator(in amps) =");
+//for a lossless line , P(avg)*I_input=P(avg)*I_load
+P_avg=V_s*I_s;//in watts
+disp(real(P_avg),"Power delivered to load (in watts) =)");
+//Real(Vs*Is)=Real(Vs*I_load)
+I_load=(P_avg/Z_l)^0.5;//in amps
+disp(real(I_load),"Current flowing in load (in amps) =)");
+
+
+
+
diff --git a/728/CH3/EX3.4/Ex3_4_ans.txt b/728/CH3/EX3.4/Ex3_4_ans.txt
new file mode 100755
index 000000000..8f47b7498
--- /dev/null
+++ b/728/CH3/EX3.4/Ex3_4_ans.txt
@@ -0,0 +1,12 @@
+Current drawn from generator(in amps) =
+
+ 0.3333333
+
+ Power delivered to load (in watts) =)
+
+ 44.444444
+
+ Current flowing in load (in amps) =)
+
+ 0.6666667
+
diff --git a/728/CH3/EX3.5/Ex3_5.sce b/728/CH3/EX3.5/Ex3_5.sce
new file mode 100755
index 000000000..75613a3ca
--- /dev/null
+++ b/728/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate VSWR & reflection coefficient.
+//Exa: 3.5
+clc;
+clear;
+close;
+Z_o=50;//in ohms
+f=300*10^6;//in Hz
+Z_l=50+%i*50;//in ohms
+wl=(3*10^8)/f;//wavelength(in meters)
+P=[(Z_l-Z_o)/(Z_l+Z_o)];
+P_mag={(real(P)^2)+(imag(P)^2)}^0.5;
+P_ang=atan(imag(P)/real(P))*180/%pi;//in degrees
+S={1+P_mag}/{1-P_mag};
+disp(P,"Reflection coefficient =");
+disp(P_mag,"Magnitude of reflection coeffcient =");
+disp(P_ang,"Angle (in degree) =");
+disp(S,"VSWR ="); \ No newline at end of file
diff --git a/728/CH3/EX3.5/Ex3_5.txt b/728/CH3/EX3.5/Ex3_5.txt
new file mode 100755
index 000000000..b822a5976
--- /dev/null
+++ b/728/CH3/EX3.5/Ex3_5.txt
@@ -0,0 +1,20 @@
+//Caption:Calculate VSWR & reflection coefficient.
+//Exa: 3.5
+clc;
+clear;
+close;
+Z_o=50;//in ohms
+f=300*10^6;//in Hz
+Z_l=50+%i*50;//in ohms
+wl=(3*10^8)/f;//wavelength(in meters)
+P=[(Z_l-Z_o)/(Z_l+Z_o)];
+P_mag={(real(P)^2)+(imag(P)^2)}^0.5;
+P_ang=atan(imag(P)/real(P))*180/%pi;//in degrees
+S={1+P_mag}/{1-P_mag};
+disp(P,"Reflection coefficient =");
+disp(P_mag,"Magnitude of reflection coeffcient =");
+disp(P_ang,"Angle (in degree) =");
+disp(S,"VSWR =");
+
+
+
diff --git a/728/CH3/EX3.5/Ex3_5_ans.txt b/728/CH3/EX3.5/Ex3_5_ans.txt
new file mode 100755
index 000000000..98df14638
--- /dev/null
+++ b/728/CH3/EX3.5/Ex3_5_ans.txt
@@ -0,0 +1,16 @@
+Reflection coefficient =
+
+ 0.2 + 0.4i
+
+ Magnitude of reflection coeffcient =
+
+ 0.4472136
+
+ Angle (in degree) =
+
+ 63.434949
+
+ VSWR =
+
+ 2.618034
+
diff --git a/728/CH3/EX3.6/Ex3_6.sce b/728/CH3/EX3.6/Ex3_6.sce
new file mode 100755
index 000000000..3e2c5adf1
--- /dev/null
+++ b/728/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,15 @@
+//Caption:Determine point of attachment & length of stub.
+//Exa 3.6
+clc;
+clear;
+close;
+Z_l=100;//in ohms
+Z_o=600;//in ohms
+f=100*10^6;//in Hz
+wl=(3*10^8)/f;
+//Position of stub is :
+m=((Z_l*Z_o)/(Z_l-Z_o))^0.5;
+pos={wl/(2*%pi)}*atan((Z_l/Z_o)^0.5);//in meters
+l={wl/(2*%pi)}*{atan(m)};//in meters
+disp(pos,"Position of stub (in meters) =");
+disp(abs(l),"Length of stub (in meters) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.6/Ex3_6.txt b/728/CH3/EX3.6/Ex3_6.txt
new file mode 100755
index 000000000..eba4f5cab
--- /dev/null
+++ b/728/CH3/EX3.6/Ex3_6.txt
@@ -0,0 +1,19 @@
+//Caption:Determine point of attachment & length of stub.
+//Exa 3.6
+clc;
+clear;
+close;
+Z_l=100;//in ohms
+Z_o=600;//in ohms
+f=100*10^6;//in Hz
+wl=(3*10^8)/f;
+//Position of stub is :
+m=((Z_l*Z_o)/(Z_l-Z_o))^0.5;
+pos={wl/(2*%pi)}*atan((Z_l/Z_o)^0.5);//in meters
+l={wl/(2*%pi)}*{atan(m)};//in meters
+disp(pos,"Position of stub (in meters) =");
+disp(abs(l),"Length of stub (in meters) =");// But answer in book is misprinted as 1.28 m
+
+
+
+
diff --git a/728/CH3/EX3.6/Ex3_6_ans.txt b/728/CH3/EX3.6/Ex3_6_ans.txt
new file mode 100755
index 000000000..35f626f8e
--- /dev/null
+++ b/728/CH3/EX3.6/Ex3_6_ans.txt
@@ -0,0 +1,8 @@
+Position of stub (in meters) =
+
+ 0.1850638
+
+ Length of stub (in meters) =
+
+ 0.7512725
+
diff --git a/728/CH3/EX3.7/Ex3_7.sce b/728/CH3/EX3.7/Ex3_7.sce
new file mode 100755
index 000000000..0a0574dec
--- /dev/null
+++ b/728/CH3/EX3.7/Ex3_7.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate terminating impedance.
+//Exa: 3.7
+clc;
+clear;
+close;
+Z_o=50;
+S=3.2;
+X_min=0.23;//in terms of wavelength(wl))
+//So :
+Z_l=Z_o*[[1-%i*S*tan(2*%pi*X_min)]/[S-%i*tan(2*%pi*X_min)]];//in ohms
+Z_lmag=[(real(Z_l)^2)+(imag(Z_l)^2)]^0.5;
+Z_lang=atan(imag(Z_l)/real(Z_l));
+disp("The load impedance");
+disp(Z_lmag,"magnitude (in ohms) =");
+disp(Z_lang*180/%pi,"angle (in degrees) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.7/Ex3_7.txt b/728/CH3/EX3.7/Ex3_7.txt
new file mode 100755
index 000000000..b22021c06
--- /dev/null
+++ b/728/CH3/EX3.7/Ex3_7.txt
@@ -0,0 +1,18 @@
+//Caption:Calculate terminating impedance.
+//Exa: 3.7
+clc;
+clear;
+close;
+Z_o=50;
+S=3.2;
+X_min=0.23;//in terms of wavelength(wl))
+//So :
+Z_l=Z_o*[[1-%i*S*tan(2*%pi*X_min)]/[S-%i*tan(2*%pi*X_min)]];//in ohms
+Z_lmag=[(real(Z_l)^2)+(imag(Z_l)^2)]^0.5;
+Z_lang=atan(imag(Z_l)/real(Z_l));
+disp("The load impedance");
+disp(Z_lmag,"magnitude (in ohms) =");// But answer in book is mispinted as 14.9
+disp(Z_lang*180/%pi,"angle (in degrees) =");
+
+
+
diff --git a/728/CH3/EX3.7/Ex3_7_ans.txt b/728/CH3/EX3.7/Ex3_7_ans.txt
new file mode 100755
index 000000000..2cf9a6bcc
--- /dev/null
+++ b/728/CH3/EX3.7/Ex3_7_ans.txt
@@ -0,0 +1,10 @@
+The load impedance
+
+ magnitude (in ohms) =
+
+ 148.45323
+
+ angle (in degrees) =
+
+ - 19.750474
+
diff --git a/728/CH3/EX3.8/Ex3_8.sce b/728/CH3/EX3.8/Ex3_8.sce
new file mode 100755
index 000000000..dd261e896
--- /dev/null
+++ b/728/CH3/EX3.8/Ex3_8.sce
@@ -0,0 +1,24 @@
+//Caption:Determine:(a)VSWR; (b)Position of 1st Vmin & Vmax; (c)Vmin & Vmax; (d)Impedance at Vmin & Vmax.
+//Exa: 3.8
+clc;
+clear;
+close;
+Z_o=50;//in ohms
+Z_l=100;//in ohms
+f=300*10^3;//in Hz
+P_l=50*10^-3;//in watts
+wl=(3*10^8)/f;
+p=(Z_l-Z_o)/(Z_l+Z_o);
+S=(1+abs(p))/(1-abs(p));
+disp(S,"VSWR =");
+//Since real Zl > Zo ,
+pos=wl/4;
+disp("First Vmax is located --->at the load ");
+disp("First Vmin is located at --->(wavelength/4)= ");
+disp(pos,"(in meters)");
+V_max=(P_l*Z_l)^0.5;
+V_min=V_max/S;
+disp(V_max,"Vmax (in volts) =");
+disp(V_min,"Vmin (in volts) =");
+disp(Z_o/S,"Zin at Vmin (in ohms) =:");
+disp(Z_o*S,"Zin at Vmax (in ohms) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.8/Ex3_8.txt b/728/CH3/EX3.8/Ex3_8.txt
new file mode 100755
index 000000000..602791344
--- /dev/null
+++ b/728/CH3/EX3.8/Ex3_8.txt
@@ -0,0 +1,25 @@
+//Caption:Determine:(a)VSWR; (b)Position of 1st Vmin & Vmax; (c)Vmin & Vmax; (d)Impedance at Vmin & Vmax.
+//Exa: 3.8
+clc;
+clear;
+close;
+Z_o=50;//in ohms
+Z_l=100;//in ohms
+f=300*10^3;//in Hz
+P_l=50*10^-3;//in watts
+wl=(3*10^8)/f;
+p=(Z_l-Z_o)/(Z_l+Z_o);
+S=(1+abs(p))/(1-abs(p));
+disp(S,"VSWR =");
+//Since real Zl > Zo ,
+pos=wl/4;
+disp("First Vmax is located --->at the load ");
+disp("First Vmin is located at --->(wavelength/4)= ");
+disp(pos,"(in meters)");
+V_max=(P_l*Z_l)^0.5;
+V_min=V_max/S;
+disp(V_max,"Vmax (in volts) =");
+disp(V_min,"Vmin (in volts) =");
+disp(Z_o/S,"Zin at Vmin (in ohms) =:");
+disp(Z_o*S,"Zin at Vmax (in ohms) =");
+
diff --git a/728/CH3/EX3.8/Ex3_8_ans.txt b/728/CH3/EX3.8/Ex3_8_ans.txt
new file mode 100755
index 000000000..1b503114b
--- /dev/null
+++ b/728/CH3/EX3.8/Ex3_8_ans.txt
@@ -0,0 +1,28 @@
+VSWR =
+
+ 2.
+
+ First Vmax is located --->at the load
+
+ First Vmin is located at --->(wavelength/4)=
+
+ (in meters)
+
+ 250.
+
+ Vmax (in volts) =
+
+ 2.236068
+
+ Vmin (in volts) =
+
+ 1.118034
+
+ Zin at Vmin (in ohms) =:
+
+ 25.
+
+ Zin at Vmax (in ohms) =
+
+ 100.
+
diff --git a/728/CH3/EX3.9/Ex3_9.sce b/728/CH3/EX3.9/Ex3_9.sce
new file mode 100755
index 000000000..8bf0a5429
--- /dev/null
+++ b/728/CH3/EX3.9/Ex3_9.sce
@@ -0,0 +1,18 @@
+//Caption:Determine in dB: (a)-reflection loss, (b)-transmission line (c)-return loss.
+//Exa: 3.9
+clc;
+clear;
+close;
+Z_o=600;//in ohm
+Z_s=50;//in ohm
+l=200;//in meter
+Z_l=500;//in ohm
+p=(Z_l-Z_o)/(Z_l+Z_o);
+ref_los=10*(log(1/(1-(abs(p))^2)))/(log(10));//in dB
+disp(ref_los,"Reflection loss (in dB) =");
+//attenuation loss= 0 dB
+//Transmisson loss = (attenuation loss)+(reflection loss) = (reflection loss)
+tran_los=ref_los;
+disp(tran_los,"Transmisson loss (in dB) =");
+ret_los=10*((log(abs(p)))/(log(10)));
+disp(ret_los,"Return loss(in dB) ="); \ No newline at end of file
diff --git a/728/CH3/EX3.9/Ex3_9.txt b/728/CH3/EX3.9/Ex3_9.txt
new file mode 100755
index 000000000..c499509b6
--- /dev/null
+++ b/728/CH3/EX3.9/Ex3_9.txt
@@ -0,0 +1,22 @@
+//Caption:Determine in dB: (a)-reflection loss, (b)-transmission line (c)-return loss.
+//Exa: 3.9
+clc;
+clear;
+close;
+Z_o=600;//in ohm
+Z_s=50;//in ohm
+l=200;//in meter
+Z_l=500;//in ohm
+p=(Z_l-Z_o)/(Z_l+Z_o);
+ref_los=10*(log(1/(1-(abs(p))^2)))/(log(10));//in dB
+disp(ref_los,"Reflection loss (in dB) =");
+//attenuation loss= 0 dB
+//Transmisson loss = (attenuation loss)+(reflection loss) = (reflection loss)
+tran_los=ref_los;
+disp(tran_los,"Transmisson loss (in dB) =");
+ret_los=10*((log(abs(p)))/(log(10)));
+disp(ret_los,"Return loss(in dB) =");
+
+
+
+
diff --git a/728/CH3/EX3.9/Ex3_9_ans.txt b/728/CH3/EX3.9/Ex3_9_ans.txt
new file mode 100755
index 000000000..1958f2c98
--- /dev/null
+++ b/728/CH3/EX3.9/Ex3_9_ans.txt
@@ -0,0 +1,12 @@
+Reflection loss (in dB) =
+
+ 0.0360412
+
+ Transmisson loss (in dB) =
+
+ 0.0360412
+
+ Return loss(in dB) =
+
+ - 10.413927
+
diff --git a/728/CH4/EX4.1/Ex4_1.sce b/728/CH4/EX4.1/Ex4_1.sce
new file mode 100755
index 000000000..22b009f8b
--- /dev/null
+++ b/728/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,18 @@
+//Caption:Calculate (i)-inductance per unit length,(ii)-capacitance per unit length,(iii)-charcteristic impedance,(iv)-velocity of propagation
+//Exa:4.1
+clc;
+clear;
+close;
+//Given:
+d=0.49;//in cm
+D=1.1;//in cm
+e_r=2.3;
+c=3*10^8;//in meter/second
+L=2*(10^-7)*log(D/d);//in Henry/meter
+C=55.56*(10^-12)*(e_r)/log(D/d);//in farad/meter
+R_o=(60/sqrt(e_r)) *log(D/d);//in ohms
+v=c/sqrt(e_r);//in meter/second
+disp(L,'Inductance per unit length(in H/m) =');
+disp(C,'Capacitance per unit length(in F/m) =');
+disp(R_o,'Characteristic Impedance (in ohms) =');
+disp(v,'Velocity of propagation (in m/s)='); \ No newline at end of file
diff --git a/728/CH4/EX4.1/Ex4_1.txt b/728/CH4/EX4.1/Ex4_1.txt
new file mode 100755
index 000000000..407e14ce5
--- /dev/null
+++ b/728/CH4/EX4.1/Ex4_1.txt
@@ -0,0 +1,18 @@
+//Caption:Calculate (i)-inductance per unit length,(ii)-capacitance per unit length,(iii)-charcteristic impedance,(iv)-velocity of propagation
+//Exa:4.1
+clc;
+clear;
+close;
+//Given:
+d=0.49;//in cm
+D=1.1;//in cm
+e_r=2.3;
+c=3*10^8;//in meter/second
+L=2*(10^-7)*log(D/d);//in Henry/meter
+C=55.56*(10^-12)*(e_r)/log(D/d);//in farad/meter
+R_o=(60/sqrt(e_r)) *log(D/d);//in ohms
+v=c/sqrt(e_r);//in meter/second
+disp(L,'Inductance per unit length(in H/m) =');
+disp(C,'Capacitance per unit length(in F/m) =');
+disp(R_o,'Characteristic Impedance (in ohms) =');
+disp(v,'Velocity of propagation (in m/s)=');
diff --git a/728/CH4/EX4.1/Ex4_1_ans.txt b/728/CH4/EX4.1/Ex4_1_ans.txt
new file mode 100755
index 000000000..2968fc1d7
--- /dev/null
+++ b/728/CH4/EX4.1/Ex4_1_ans.txt
@@ -0,0 +1,15 @@
+Inductance per unit length(in H/m) =
+
+ 0.0000002
+
+ Capacitance per unit length(in F/m) =
+
+ 1.580D-10
+
+ Characteristic Impedance (in ohms) =
+
+ 31.992879
+
+ Velocity of propagation (in m/s)=
+
+ 1.978D+08
diff --git a/728/CH4/EX4.10/Ex4_10.sce b/728/CH4/EX4.10/Ex4_10.sce
new file mode 100755
index 000000000..3c1fc0b5c
--- /dev/null
+++ b/728/CH4/EX4.10/Ex4_10.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate (i)-required size of guide,(ii)-frequencies that can be used for this mode of propagation
+//Exa:4.10
+clc;
+clear;
+close;
+wl_c=10;//in cm
+c=3*10^10;//in cm/s
+r=wl_c/(2*%pi/1.841);//in cm
+area=%pi*r^2;//in sq. cm
+f_c=c/wl_c;
+disp(r,'Radius of circular waveguide(in cm) =');
+disp(area,'Area of cross-section of circular waveguide(in cm) =');
+disp('Frequency above');
+disp(f_c);
+disp('can be propagated'); \ No newline at end of file
diff --git a/728/CH4/EX4.10/Ex4_10.txt b/728/CH4/EX4.10/Ex4_10.txt
new file mode 100755
index 000000000..4815e1ab4
--- /dev/null
+++ b/728/CH4/EX4.10/Ex4_10.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate (i)-required size of guide,(ii)-frequencies that can be used for this mode of propagation
+//Exa:4.10
+clc;
+clear;
+close;
+wl_c=10;//in cm
+c=3*10^10;//in cm/s
+r=wl_c/(2*%pi/1.841);//in cm
+area=%pi*r^2;//in sq. cm
+f_c=c/wl_c;
+disp(r,'Radius of circular waveguide(in cm) =');
+disp(area,'Area of cross-section of circular waveguide(in cm) =');
+disp('Frequency above');
+disp(f_c);
+disp('can be propagated');
diff --git a/728/CH4/EX4.10/Ex4_10_ans.txt b/728/CH4/EX4.10/Ex4_10_ans.txt
new file mode 100755
index 000000000..29e96150f
--- /dev/null
+++ b/728/CH4/EX4.10/Ex4_10_ans.txt
@@ -0,0 +1,13 @@
+ Radius of circular waveguide(in cm) =
+
+ 2.9300425
+
+ Area of cross-section of circular waveguide(in cm) =
+
+ 26.971041
+
+ Frequency above
+
+ 3.000D+09
+
+ can be propagated
diff --git a/728/CH4/EX4.11/Ex4_11.sce b/728/CH4/EX4.11/Ex4_11.sce
new file mode 100755
index 000000000..1f3099ebe
--- /dev/null
+++ b/728/CH4/EX4.11/Ex4_11.sce
@@ -0,0 +1,29 @@
+//Caption:Find all modes that can propagate at 5000MHz.
+//Exa:4_11
+clc;
+clear;
+close;
+a=4;//in cm
+b=3;//in cm
+f=5*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+//For TE waves:
+wl_c_TE01=2*b;//for TE01
+wl_c_TE10=2*a;//for TE10
+wl_c_TE11=2*a*b/sqrt(a^2+b^2);//for TE11
+if(wl_c_TE01>wl_o)
+ disp('TE01 can propagate');
+else
+ disp('TE01 cannot propagate');
+end
+if(wl_c_TE10>wl_o)
+ disp('TE10 can propagate');
+else
+ disp('TE10 cannot propagate');
+end
+if(wl_c_TE11>wl_o)
+ disp('TE11 can propagate');
+else
+ disp('TE11 cannot propagate');
+end \ No newline at end of file
diff --git a/728/CH4/EX4.11/Ex4_11.txt b/728/CH4/EX4.11/Ex4_11.txt
new file mode 100755
index 000000000..b9df1fc2e
--- /dev/null
+++ b/728/CH4/EX4.11/Ex4_11.txt
@@ -0,0 +1,29 @@
+//Caption:Find all modes that can propagate at 5000MHz.
+//Exa:4_11
+clc;
+clear;
+close;
+a=4;//in cm
+b=3;//in cm
+f=5*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+//For TE waves:
+wl_c_TE01=2*b;//for TE01
+wl_c_TE10=2*a;//for TE10
+wl_c_TE11=2*a*b/sqrt(a^2+b^2);//for TE11
+if(wl_c_TE01>wl_o)
+ disp('TE01 can propagate');
+else
+ disp('TE01 cannot propagate');
+end
+if(wl_c_TE10>wl_o)
+ disp('TE10 can propagate');
+else
+ disp('TE10 cannot propagate');
+end
+if(wl_c_TE11>wl_o)
+ disp('TE11 can propagate');
+else
+ disp('TE11 cannot propagate');
+end
diff --git a/728/CH4/EX4.11/Ex4_11_ans.txt b/728/CH4/EX4.11/Ex4_11_ans.txt
new file mode 100755
index 000000000..abfe71cb6
--- /dev/null
+++ b/728/CH4/EX4.11/Ex4_11_ans.txt
@@ -0,0 +1,6 @@
+ TE01 cannot propagate
+
+ TE10 can propagate
+
+ TE11 cannot propagate
+
diff --git a/728/CH4/EX4.12/Ex4_11.sce b/728/CH4/EX4.12/Ex4_11.sce
new file mode 100755
index 000000000..1f3099ebe
--- /dev/null
+++ b/728/CH4/EX4.12/Ex4_11.sce
@@ -0,0 +1,29 @@
+//Caption:Find all modes that can propagate at 5000MHz.
+//Exa:4_11
+clc;
+clear;
+close;
+a=4;//in cm
+b=3;//in cm
+f=5*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+//For TE waves:
+wl_c_TE01=2*b;//for TE01
+wl_c_TE10=2*a;//for TE10
+wl_c_TE11=2*a*b/sqrt(a^2+b^2);//for TE11
+if(wl_c_TE01>wl_o)
+ disp('TE01 can propagate');
+else
+ disp('TE01 cannot propagate');
+end
+if(wl_c_TE10>wl_o)
+ disp('TE10 can propagate');
+else
+ disp('TE10 cannot propagate');
+end
+if(wl_c_TE11>wl_o)
+ disp('TE11 can propagate');
+else
+ disp('TE11 cannot propagate');
+end \ No newline at end of file
diff --git a/728/CH4/EX4.12/Ex4_11.txt b/728/CH4/EX4.12/Ex4_11.txt
new file mode 100755
index 000000000..b9df1fc2e
--- /dev/null
+++ b/728/CH4/EX4.12/Ex4_11.txt
@@ -0,0 +1,29 @@
+//Caption:Find all modes that can propagate at 5000MHz.
+//Exa:4_11
+clc;
+clear;
+close;
+a=4;//in cm
+b=3;//in cm
+f=5*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+//For TE waves:
+wl_c_TE01=2*b;//for TE01
+wl_c_TE10=2*a;//for TE10
+wl_c_TE11=2*a*b/sqrt(a^2+b^2);//for TE11
+if(wl_c_TE01>wl_o)
+ disp('TE01 can propagate');
+else
+ disp('TE01 cannot propagate');
+end
+if(wl_c_TE10>wl_o)
+ disp('TE10 can propagate');
+else
+ disp('TE10 cannot propagate');
+end
+if(wl_c_TE11>wl_o)
+ disp('TE11 can propagate');
+else
+ disp('TE11 cannot propagate');
+end
diff --git a/728/CH4/EX4.12/Ex4_11_ans.txt b/728/CH4/EX4.12/Ex4_11_ans.txt
new file mode 100755
index 000000000..abfe71cb6
--- /dev/null
+++ b/728/CH4/EX4.12/Ex4_11_ans.txt
@@ -0,0 +1,6 @@
+ TE01 cannot propagate
+
+ TE10 can propagate
+
+ TE11 cannot propagate
+
diff --git a/728/CH4/EX4.13/Ex4_12.sce b/728/CH4/EX4.13/Ex4_12.sce
new file mode 100755
index 000000000..cc6ebf8c1
--- /dev/null
+++ b/728/CH4/EX4.13/Ex4_12.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-cutoff wavelength,(ii)-cutoff frequency,(iii)-wavelength in guide
+//Exa:4.12
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+d=4;//in cm
+r=d/2;//in cm
+wl_c=2*%pi*r/1.841;//in cm
+f_c=c/wl_c;
+f_signal=5*10^9;//in Hz
+wl_o=c/f_signal;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+disp(wl_c,'Cut-off wavelength (in cm) =');
+disp(f_c,'Cut-off frequency (in Hz) =');
+disp(wl_g,'Guide wavelength (in cm) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.13/Ex4_12.txt b/728/CH4/EX4.13/Ex4_12.txt
new file mode 100755
index 000000000..65db64e03
--- /dev/null
+++ b/728/CH4/EX4.13/Ex4_12.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-cutoff wavelength,(ii)-cutoff frequency,(iii)-wavelength in guide
+//Exa:4.12
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+d=4;//in cm
+r=d/2;//in cm
+wl_c=2*%pi*r/1.841;//in cm
+f_c=c/wl_c;
+f_signal=5*10^9;//in Hz
+wl_o=c/f_signal;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+disp(wl_c,'Cut-off wavelength (in cm) =');
+disp(f_c,'Cut-off frequency (in Hz) =');
+disp(wl_g,'Guide wavelength (in cm) =');
diff --git a/728/CH4/EX4.13/Ex4_12_ans.txt b/728/CH4/EX4.13/Ex4_12_ans.txt
new file mode 100755
index 000000000..fc4c9ed0d
--- /dev/null
+++ b/728/CH4/EX4.13/Ex4_12_ans.txt
@@ -0,0 +1,12 @@
+Cut-off wavelength (in cm) =
+
+ 6.8258396
+
+ Cut-off frequency (in Hz) =
+
+ 4.395D+09
+
+ Guide wavelength (in cm) =
+
+ 12.583938
+
diff --git a/728/CH4/EX4.14/Ex4_14.sce b/728/CH4/EX4.14/Ex4_14.sce
new file mode 100755
index 000000000..716271d08
--- /dev/null
+++ b/728/CH4/EX4.14/Ex4_14.sce
@@ -0,0 +1,22 @@
+//Caption:Calculate (i)-guide wavelength,(ii)-phase constant,(iii)-phase velocity for dominant mode
+//Exa:4.14
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+a=5;//in cm
+b=2.5;//in cm
+wl_o=4.5;//in cm
+//For TE10 mode:
+wl_c=2*a;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+V_p=c/sqrt(1-(wl_o/wl_c)^2);
+w=2*%pi*c/wl_o;
+w_c=2*%pi*c/wl_c;
+b=sqrt(w^2-w_c^2)/c;
+disp(wl_g,'Guide wavelength (in cm) =');
+disp(b,'Phase constant =');
+disp(V_p,'Phase velocity (in cm/s) =');
+
+//answer in book is wrongly written as guide wavelength =7.803 cm
+//answer in book is wrongly written as Phase velocity = 5.22*10^10 cm/s \ No newline at end of file
diff --git a/728/CH4/EX4.14/Ex4_14.txt b/728/CH4/EX4.14/Ex4_14.txt
new file mode 100755
index 000000000..1c6d6a884
--- /dev/null
+++ b/728/CH4/EX4.14/Ex4_14.txt
@@ -0,0 +1,22 @@
+//Caption:Calculate (i)-guide wavelength,(ii)-phase constant,(iii)-phase velocity for dominant mode
+//Exa:4.14
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+a=5;//in cm
+b=2.5;//in cm
+wl_o=4.5;//in cm
+//For TE10 mode:
+wl_c=2*a;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+V_p=c/sqrt(1-(wl_o/wl_c)^2);
+w=2*%pi*c/wl_o;
+w_c=2*%pi*c/wl_c;
+b=sqrt(w^2-w_c^2)/c;
+disp(wl_g,'Guide wavelength (in cm) =');
+disp(b,'Phase constant =');
+disp(V_p,'Phase velocity (in cm/s) =');
+
+//answer in book is wrongly written as guide wavelength =7.803 cm
+//answer in book is wrongly written as Phase velocity = 5.22*10^10 cm/s
diff --git a/728/CH4/EX4.14/Ex4_14_ans.txt b/728/CH4/EX4.14/Ex4_14_ans.txt
new file mode 100755
index 000000000..a203e377e
--- /dev/null
+++ b/728/CH4/EX4.14/Ex4_14_ans.txt
@@ -0,0 +1,12 @@
+
+ Guide wavelength (in cm) =
+
+ 5.0390326
+
+ Phase constant =
+
+ 1.2469031
+
+ Phase velocity (in cm/s) =
+
+ 3.359D+10
diff --git a/728/CH4/EX4.15/Ex4_15.sce b/728/CH4/EX4.15/Ex4_15.sce
new file mode 100755
index 000000000..fc616db2b
--- /dev/null
+++ b/728/CH4/EX4.15/Ex4_15.sce
@@ -0,0 +1,45 @@
+//Caption:Calculate what modes propagate at free space wavelength of (i)10 cm,(ii)5 cm
+//Exa:4.15
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+wl_c_TE10=16;//Critical wavelength of TE10
+wl_c_TM11=7.16;//Critical wavelength of TM11
+wl_c_TM21=5.6;//Critical wavelength of TM21
+//For (i): 10 cm
+wl_o=10;//in cm
+disp(wl_o,'For free space wavelength (in cm) =');
+if(wl_c_TE10>wl_o)
+ disp(' TE10 can propagate');
+else
+ disp(' TE10 cannot propagate');
+end
+if(wl_c_TM11>wl_o)
+ disp(' TM11 can propagate');
+else
+ disp(' TM11 cannot propagate');
+end
+if(wl_c_TM21>wl_o)
+ disp(' TM21 can propagate');
+else
+ disp(' TM21 cannot propagate');
+end
+//For (ii): 5 cm
+wl_o=5;//in cm
+disp(wl_o,'For free space wavelength (in cm) =');
+if(wl_c_TE10>wl_o)
+ disp(' TE10 can propagate');
+else
+ disp(' TE10 cannot propagate');
+end
+if(wl_c_TM11>wl_o)
+ disp(' TM11 can propagate');
+else
+ disp(' TM11 cannot propagate');
+end
+if(wl_c_TM21>wl_o)
+ disp(' TM21 can propagate');
+else
+ disp(' TM21 cannot propagate');
+end \ No newline at end of file
diff --git a/728/CH4/EX4.15/Ex4_15.txt b/728/CH4/EX4.15/Ex4_15.txt
new file mode 100755
index 000000000..5b737bc43
--- /dev/null
+++ b/728/CH4/EX4.15/Ex4_15.txt
@@ -0,0 +1,45 @@
+//Caption:Calculate what modes propagate at free space wavelength of (i)10 cm,(ii)5 cm
+//Exa:4.15
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+wl_c_TE10=16;//Critical wavelength of TE10
+wl_c_TM11=7.16;//Critical wavelength of TM11
+wl_c_TM21=5.6;//Critical wavelength of TM21
+//For (i): 10 cm
+wl_o=10;//in cm
+disp(wl_o,'For free space wavelength (in cm) =');
+if(wl_c_TE10>wl_o)
+ disp(' TE10 can propagate');
+else
+ disp(' TE10 cannot propagate');
+end
+if(wl_c_TM11>wl_o)
+ disp(' TM11 can propagate');
+else
+ disp(' TM11 cannot propagate');
+end
+if(wl_c_TM21>wl_o)
+ disp(' TM21 can propagate');
+else
+ disp(' TM21 cannot propagate');
+end
+//For (ii): 5 cm
+wl_o=5;//in cm
+disp(wl_o,'For free space wavelength (in cm) =');
+if(wl_c_TE10>wl_o)
+ disp(' TE10 can propagate');
+else
+ disp(' TE10 cannot propagate');
+end
+if(wl_c_TM11>wl_o)
+ disp(' TM11 can propagate');
+else
+ disp(' TM11 cannot propagate');
+end
+if(wl_c_TM21>wl_o)
+ disp(' TM21 can propagate');
+else
+ disp(' TM21 cannot propagate');
+end
diff --git a/728/CH4/EX4.15/Ex4_15_ans.txt b/728/CH4/EX4.15/Ex4_15_ans.txt
new file mode 100755
index 000000000..b5cf9b47a
--- /dev/null
+++ b/728/CH4/EX4.15/Ex4_15_ans.txt
@@ -0,0 +1,20 @@
+
+ For free space wavelength (in cm) =
+
+ 10.
+
+ TE10 can propagate
+
+ TM11 cannot propagate
+
+ TM21 cannot propagate
+
+ For free space wavelength (in cm) =
+
+ 5.
+
+ TE10 can propagate
+
+ TM11 can propagate
+
+ TM21 can propagate
diff --git a/728/CH4/EX4.16/Ex4_16.sce b/728/CH4/EX4.16/Ex4_16.sce
new file mode 100755
index 000000000..f930728d2
--- /dev/null
+++ b/728/CH4/EX4.16/Ex4_16.sce
@@ -0,0 +1,16 @@
+//Caption:Determine the charcteristic wave impedance
+//Exa:4.16
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=10*10^9;//in Hz
+a=3;//in cm
+b=2;//in cm
+n=120*%pi;
+wl_o=c/f;
+wl_c=2*a*b/sqrt(a^2+b^2);
+Z_TM=n*sqrt(1-(wl_o/wl_c)^2);
+disp(Z_TM,'Characteristic impedance (in ohms) =');
+
+//answer in book is wrongly written as 61.618 ohms \ No newline at end of file
diff --git a/728/CH4/EX4.16/Ex4_16.txt b/728/CH4/EX4.16/Ex4_16.txt
new file mode 100755
index 000000000..5031bff93
--- /dev/null
+++ b/728/CH4/EX4.16/Ex4_16.txt
@@ -0,0 +1,16 @@
+//Caption:Determine the charcteristic wave impedance
+//Exa:4.16
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=10*10^9;//in Hz
+a=3;//in cm
+b=2;//in cm
+n=120*%pi;
+wl_o=c/f;
+wl_c=2*a*b/sqrt(a^2+b^2);
+Z_TM=n*sqrt(1-(wl_o/wl_c)^2);
+disp(Z_TM,'Characteristic impedance (in ohms) =');
+
+//answer in book is wrongly written as 61.618 ohms
diff --git a/728/CH4/EX4.16/Ex4_16_ans.txt b/728/CH4/EX4.16/Ex4_16_ans.txt
new file mode 100755
index 000000000..9e4469d24
--- /dev/null
+++ b/728/CH4/EX4.16/Ex4_16_ans.txt
@@ -0,0 +1,4 @@
+ Characteristic impedance (in ohms) =
+
+ 163.24194
+
diff --git a/728/CH4/EX4.17/EX4_17_ans.txt b/728/CH4/EX4.17/EX4_17_ans.txt
new file mode 100755
index 000000000..e0bf44269
--- /dev/null
+++ b/728/CH4/EX4.17/EX4_17_ans.txt
@@ -0,0 +1,8 @@
+ Diameter of waveguide (in cm) =
+
+ 3.6625531
+
+ Guide wavelength (in cm) =
+
+ 8.3333333
+
diff --git a/728/CH4/EX4.17/Ex4_17.sce b/728/CH4/EX4.17/Ex4_17.sce
new file mode 100755
index 000000000..3169154d7
--- /dev/null
+++ b/728/CH4/EX4.17/Ex4_17.sce
@@ -0,0 +1,14 @@
+//Caption:Determine the diameter of waveguide & guide wavelength
+//Exa:4.17
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=6*10^9;//in Hz
+f_c=0.8*f;
+wl_c=c/f_c;
+D=1.841*wl_c/%pi;
+wl_o=c/f;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+disp(D,'Diameter of waveguide (in cm) =');
+disp(wl_g,'Guide wavelength (in cm) =');
diff --git a/728/CH4/EX4.17/Ex4_17.txt b/728/CH4/EX4.17/Ex4_17.txt
new file mode 100755
index 000000000..3169154d7
--- /dev/null
+++ b/728/CH4/EX4.17/Ex4_17.txt
@@ -0,0 +1,14 @@
+//Caption:Determine the diameter of waveguide & guide wavelength
+//Exa:4.17
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=6*10^9;//in Hz
+f_c=0.8*f;
+wl_c=c/f_c;
+D=1.841*wl_c/%pi;
+wl_o=c/f;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+disp(D,'Diameter of waveguide (in cm) =');
+disp(wl_g,'Guide wavelength (in cm) =');
diff --git a/728/CH4/EX4.18/Ex4_18.sce b/728/CH4/EX4.18/Ex4_18.sce
new file mode 100755
index 000000000..2dfa1eda2
--- /dev/null
+++ b/728/CH4/EX4.18/Ex4_18.sce
@@ -0,0 +1,20 @@
+//Caption:Show TE01 mode propagates under given conditions
+//Exa:4.18
+clc;
+clear;
+close;
+a=1.5;//in cm
+b=1;//in cm
+e_r=4;//dielectric
+c=3*10^10;//in cm/s
+wl_c=2*b;
+f_c=c/wl_c;
+f_imp=6*10^9;//impressed frequency (in Hz)
+wl_air=c/f_imp;
+//Inserting dielectric:
+wl_dielec=wl_air/sqrt(e_r);
+if(wl_dielec>wl_c)
+ disp(' TE01 can propagate');
+else
+ disp(' TE01 cannot propagate');
+end \ No newline at end of file
diff --git a/728/CH4/EX4.18/Ex4_18.txt b/728/CH4/EX4.18/Ex4_18.txt
new file mode 100755
index 000000000..8fdb9e88d
--- /dev/null
+++ b/728/CH4/EX4.18/Ex4_18.txt
@@ -0,0 +1,20 @@
+//Caption:Show TE01 mode propagates under given conditions
+//Exa:4.18
+clc;
+clear;
+close;
+a=1.5;//in cm
+b=1;//in cm
+e_r=4;//dielectric
+c=3*10^10;//in cm/s
+wl_c=2*b;
+f_c=c/wl_c;
+f_imp=6*10^9;//impressed frequency (in Hz)
+wl_air=c/f_imp;
+//Inserting dielectric:
+wl_dielec=wl_air/sqrt(e_r);
+if(wl_dielec>wl_c)
+ disp(' TE01 can propagate');
+else
+ disp(' TE01 cannot propagate');
+end
diff --git a/728/CH4/EX4.18/Ex4_18_ans.txt b/728/CH4/EX4.18/Ex4_18_ans.txt
new file mode 100755
index 000000000..ce6a79d65
--- /dev/null
+++ b/728/CH4/EX4.18/Ex4_18_ans.txt
@@ -0,0 +1,2 @@
+ TE01 can propagate
+
diff --git a/728/CH4/EX4.19/Ex4_19.sce b/728/CH4/EX4.19/Ex4_19.sce
new file mode 100755
index 000000000..f4a9dbf88
--- /dev/null
+++ b/728/CH4/EX4.19/Ex4_19.sce
@@ -0,0 +1,21 @@
+//Caption:Calculate the amount of attenuation if signal of frequency is 6GHz
+//Exa:4.19
+clc;
+clear;
+close;
+u=4*%pi*10^-7;
+e=8.85*10^-12;
+c=3*10^10;//in cm/s
+f=6*10^9;//in Hz
+a=1.5;//in cm
+b=1;//in cm
+//For TE10 mode:
+m=1;
+n=0;
+wl_c=2*a;
+f_c=c/wl_c;
+t_1=(m*%pi/a)^2;
+t_2=(n*%pi/b)^2;
+t_3=(((2*%pi*f)^2)*u*e);
+a=sqrt(t_1+t_2-t_3);//in neper/m
+disp(a*20/log(10),'Attenuation (in dB/m) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.19/Ex4_19.txt b/728/CH4/EX4.19/Ex4_19.txt
new file mode 100755
index 000000000..71689b8ee
--- /dev/null
+++ b/728/CH4/EX4.19/Ex4_19.txt
@@ -0,0 +1,21 @@
+//Caption:Calculate the amount of attenuation if signal of frequency is 6GHz
+//Exa:4.19
+clc;
+clear;
+close;
+u=4*%pi*10^-7;
+e=8.85*10^-12;
+c=3*10^10;//in cm/s
+f=6*10^9;//in Hz
+a=1.5;//in cm
+b=1;//in cm
+//For TE10 mode:
+m=1;
+n=0;
+wl_c=2*a;
+f_c=c/wl_c;
+t_1=(m*%pi/a)^2;
+t_2=(n*%pi/b)^2;
+t_3=(((2*%pi*f)^2)*u*e);
+a=sqrt(t_1+t_2-t_3);//in neper/m
+disp(a*20/log(10),'Attenuation (in dB/m) =');
diff --git a/728/CH4/EX4.19/Ex4_19_ans.txt b/728/CH4/EX4.19/Ex4_19_ans.txt
new file mode 100755
index 000000000..60c025115
--- /dev/null
+++ b/728/CH4/EX4.19/Ex4_19_ans.txt
@@ -0,0 +1,5 @@
+
+ Attenuation (in dB/m) =
+
+ 1091.8468i
+ \ No newline at end of file
diff --git a/728/CH4/EX4.2/Ex4_2.sce b/728/CH4/EX4.2/Ex4_2.sce
new file mode 100755
index 000000000..1993972cf
--- /dev/null
+++ b/728/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,25 @@
+//Caption:Calculate the attenuation, phase constants,phase velocity,relative permittivity,power loss.
+//Exa:4.2
+clc;
+clear;
+close;
+R=0.05;//in ohms
+G=0;
+l=50;//in meter
+e=2.3;//dielectric constant
+c=3*10^8;//in m/s
+L=2*(10^-7);//from Exa:4.1
+C=1.58*(10^-10);//from Exa:4.1
+P_in=480;//in watts
+f=3*10^9;//in hertz
+Z_o=sqrt(L/C);
+a=R/Z_o;//in Np/m
+b=2*%pi*f*sqrt(L*C);//in rad/m
+V_p=1/sqrt(L*C);
+e_r=(c/V_p)^2;
+P_loss=P_in*2*l;
+disp(a,'Atteneuation (in Np/m) =');
+disp(b,'Phase constant (in rad/m) =');
+disp(V_p,'Phase velocity (in m/s) =');
+disp(e_r,'Relative permittivity =');
+disp(P_loss,'Power loss (in watts) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.2/Ex4_2.txt b/728/CH4/EX4.2/Ex4_2.txt
new file mode 100755
index 000000000..72254bd36
--- /dev/null
+++ b/728/CH4/EX4.2/Ex4_2.txt
@@ -0,0 +1,25 @@
+//Caption:Calculate the attenuation, phase constants,phase velocity,relative permittivity,power loss.
+//Exa:4.2
+clc;
+clear;
+close;
+R=0.05;//in ohms
+G=0;
+l=50;//in meter
+e=2.3;//dielectric constant
+c=3*10^8;//in m/s
+L=2*(10^-7);//from Exa:4.1
+C=1.58*(10^-10);//from Exa:4.1
+P_in=480;//in watts
+f=3*10^9;//in hertz
+Z_o=sqrt(L/C);
+a=R/Z_o;//in Np/m
+b=2*%pi*f*sqrt(L*C);//in rad/m
+V_p=1/sqrt(L*C);
+e_r=(c/V_p)^2;
+P_loss=P_in*2*l;
+disp(a,'Atteneuation (in Np/m) =');
+disp(b,'Phase constant (in rad/m) =');
+disp(V_p,'Phase velocity (in m/s) =');
+disp(e_r,'Relative permittivity =');
+disp(P_loss,'Power loss (in watts) =');
diff --git a/728/CH4/EX4.2/Ex4_2_ans.txt b/728/CH4/EX4.2/Ex4_2_ans.txt
new file mode 100755
index 000000000..82c5b17e8
--- /dev/null
+++ b/728/CH4/EX4.2/Ex4_2_ans.txt
@@ -0,0 +1,19 @@
+Atteneuation (in Np/m) =
+
+ 0.0014053
+
+ Phase constant (in rad/m) =
+
+ 105.96066
+
+ Phase velocity (in m/s) =
+
+ 1.779D+08
+
+ Relative permittivity =
+
+ 2.844
+
+ Power loss (in watts) =
+
+ 48000.
diff --git a/728/CH4/EX4.20/Ex4_20.sce b/728/CH4/EX4.20/Ex4_20.sce
new file mode 100755
index 000000000..eeba92d30
--- /dev/null
+++ b/728/CH4/EX4.20/Ex4_20.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate the maximum power handling capacity
+//Exa:4.21
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+a=3;//in cm
+b=1;//in cm
+E_max=3000;//in V/cm
+wl_o=c/f;
+wl_c=2*a;//in TE10
+wl_g=ceil (wl_o/sqrt(1-(wl_o/wl_c)^2));
+P_max=(6.63*10^-4)*E_max^2*a*b*(wl_o/wl_g);
+disp(P_max/1000,'Maximum power for rectangular waveguide (in kilowatts)='); \ No newline at end of file
diff --git a/728/CH4/EX4.20/Ex4_20.txt b/728/CH4/EX4.20/Ex4_20.txt
new file mode 100755
index 000000000..eeba92d30
--- /dev/null
+++ b/728/CH4/EX4.20/Ex4_20.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate the maximum power handling capacity
+//Exa:4.21
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+a=3;//in cm
+b=1;//in cm
+E_max=3000;//in V/cm
+wl_o=c/f;
+wl_c=2*a;//in TE10
+wl_g=ceil (wl_o/sqrt(1-(wl_o/wl_c)^2));
+P_max=(6.63*10^-4)*E_max^2*a*b*(wl_o/wl_g);
+disp(P_max/1000,'Maximum power for rectangular waveguide (in kilowatts)='); \ No newline at end of file
diff --git a/728/CH4/EX4.20/Ex4_20_ans.txt b/728/CH4/EX4.20/Ex4_20_ans.txt
new file mode 100755
index 000000000..9ea8d18d2
--- /dev/null
+++ b/728/CH4/EX4.20/Ex4_20_ans.txt
@@ -0,0 +1,4 @@
+
+ Maximum power for rectangular waveguide (in kilowatts)=
+
+ 11.934 \ No newline at end of file
diff --git a/728/CH4/EX4.21/Ex4_21.sce b/728/CH4/EX4.21/Ex4_21.sce
new file mode 100755
index 000000000..3155c151b
--- /dev/null
+++ b/728/CH4/EX4.21/Ex4_21.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate the maximum power
+//Exa:4.21
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+E_max=300;//in V/cm
+d=5;
+wl_o=c/f;
+//For TE11
+wl_c=d*%pi/1.841;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+P_max=0.498*E_max^2*d^2*(wl_o/wl_g);
+disp(P_max,'Maximum power (in watts) =');
diff --git a/728/CH4/EX4.21/Ex4_21.txt b/728/CH4/EX4.21/Ex4_21.txt
new file mode 100755
index 000000000..3155c151b
--- /dev/null
+++ b/728/CH4/EX4.21/Ex4_21.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate the maximum power
+//Exa:4.21
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+E_max=300;//in V/cm
+d=5;
+wl_o=c/f;
+//For TE11
+wl_c=d*%pi/1.841;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+P_max=0.498*E_max^2*d^2*(wl_o/wl_g);
+disp(P_max,'Maximum power (in watts) =');
diff --git a/728/CH4/EX4.21/Ex4_21_ans.txt b/728/CH4/EX4.21/Ex4_21_ans.txt
new file mode 100755
index 000000000..9cc2a8ad6
--- /dev/null
+++ b/728/CH4/EX4.21/Ex4_21_ans.txt
@@ -0,0 +1,5 @@
+
+ Maximum power (in watts) =
+
+ 1031453.7
+ \ No newline at end of file
diff --git a/728/CH4/EX4.22/Ex4_22.sce b/728/CH4/EX4.22/Ex4_22.sce
new file mode 100755
index 000000000..4f50012da
--- /dev/null
+++ b/728/CH4/EX4.22/Ex4_22.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate the peak value of electric feild occuring in the waveguide
+//Exa:4.22
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=30*10^9;//inHz
+a=1;//in cm
+b=1;
+P_max=746;//in watts
+wl_o=c/f;
+wl_c=2*a;
+Z=120*%pi/sqrt(1-(wl_o/wl_c)^2);
+E_max=sqrt(P_max*4*Z/(a*b/10000));
+disp(E_max/1000,'Peak value of electric field (in kV/m) =');
diff --git a/728/CH4/EX4.22/Ex4_22.txt b/728/CH4/EX4.22/Ex4_22.txt
new file mode 100755
index 000000000..4f50012da
--- /dev/null
+++ b/728/CH4/EX4.22/Ex4_22.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate the peak value of electric feild occuring in the waveguide
+//Exa:4.22
+clc;
+clear;
+close;
+c=3*10^10;//in cm/s
+f=30*10^9;//inHz
+a=1;//in cm
+b=1;
+P_max=746;//in watts
+wl_o=c/f;
+wl_c=2*a;
+Z=120*%pi/sqrt(1-(wl_o/wl_c)^2);
+E_max=sqrt(P_max*4*Z/(a*b/10000));
+disp(E_max/1000,'Peak value of electric field (in kV/m) =');
diff --git a/728/CH4/EX4.22/Ex4_22_ans.txt b/728/CH4/EX4.22/Ex4_22_ans.txt
new file mode 100755
index 000000000..faaf679e6
--- /dev/null
+++ b/728/CH4/EX4.22/Ex4_22_ans.txt
@@ -0,0 +1,5 @@
+
+ Peak value of electric field (in kV/m) =
+
+ 113.97239
+ \ No newline at end of file
diff --git a/728/CH4/EX4.23/Ex4_23.sce b/728/CH4/EX4.23/Ex4_23.sce
new file mode 100755
index 000000000..4114660b2
--- /dev/null
+++ b/728/CH4/EX4.23/Ex4_23.sce
@@ -0,0 +1,13 @@
+//Caption:Calculate the breakdown power of air filled rectangular waveguide for dominant mode at 9.375 GHz.
+//Exa:4.23
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=2.3;//in cm
+b=1;//in cm
+f=9.375*10^9;//in Hz
+wl_o=c/f;
+P_bd_TE11=597*2.3*1*{1-{wl_o/(2*a)}^2}^0.5;
+disp(P_bd_TE11,'Breakdown power for dominant mode (in kW) =');
diff --git a/728/CH4/EX4.23/Ex4_23.txt b/728/CH4/EX4.23/Ex4_23.txt
new file mode 100755
index 000000000..4114660b2
--- /dev/null
+++ b/728/CH4/EX4.23/Ex4_23.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate the breakdown power of air filled rectangular waveguide for dominant mode at 9.375 GHz.
+//Exa:4.23
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=2.3;//in cm
+b=1;//in cm
+f=9.375*10^9;//in Hz
+wl_o=c/f;
+P_bd_TE11=597*2.3*1*{1-{wl_o/(2*a)}^2}^0.5;
+disp(P_bd_TE11,'Breakdown power for dominant mode (in kW) =');
diff --git a/728/CH4/EX4.23/Ex4_23_ans.txt b/728/CH4/EX4.23/Ex4_23_ans.txt
new file mode 100755
index 000000000..c9e104875
--- /dev/null
+++ b/728/CH4/EX4.23/Ex4_23_ans.txt
@@ -0,0 +1,4 @@
+Breakdown power for dominant mode (in kW) =
+
+ 986.40589
+ \ No newline at end of file
diff --git a/728/CH4/EX4.24/Ex4_24.sce b/728/CH4/EX4.24/Ex4_24.sce
new file mode 100755
index 000000000..102fc9259
--- /dev/null
+++ b/728/CH4/EX4.24/Ex4_24.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate the breakdown power of circular waveguide
+//Exa:4.24
+clc;
+clear;
+close;
+//Given:
+d=5;//in cm
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+//Dominant mode is TE11:
+wl_o=c/f;
+wl_c=%pi*d/1.841;
+f_c=c/wl_c;
+P_bd_TE11=1790*(d/2)^2*[1-{f_c/f}^2]^0.5;
+disp(P_bd_TE11/1000,'Breakdown power (in kW) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.24/Ex4_24.txt b/728/CH4/EX4.24/Ex4_24.txt
new file mode 100755
index 000000000..102fc9259
--- /dev/null
+++ b/728/CH4/EX4.24/Ex4_24.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate the breakdown power of circular waveguide
+//Exa:4.24
+clc;
+clear;
+close;
+//Given:
+d=5;//in cm
+c=3*10^10;//in cm/s
+f=9*10^9;//inHz
+//Dominant mode is TE11:
+wl_o=c/f;
+wl_c=%pi*d/1.841;
+f_c=c/wl_c;
+P_bd_TE11=1790*(d/2)^2*[1-{f_c/f}^2]^0.5;
+disp(P_bd_TE11/1000,'Breakdown power (in kW) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.24/Ex4_24_ans.txt b/728/CH4/EX4.24/Ex4_24_ans.txt
new file mode 100755
index 000000000..04341cde4
--- /dev/null
+++ b/728/CH4/EX4.24/Ex4_24_ans.txt
@@ -0,0 +1,5 @@
+
+ Breakdown power (in kW) =
+
+ 10.298427
+ \ No newline at end of file
diff --git a/728/CH4/EX4.3/EX4_3_ans.txt b/728/CH4/EX4.3/EX4_3_ans.txt
new file mode 100755
index 000000000..c69a000b6
--- /dev/null
+++ b/728/CH4/EX4.3/EX4_3_ans.txt
@@ -0,0 +1,3 @@
+ Breakdown Power (in kW) =
+
+ 17560.256
diff --git a/728/CH4/EX4.3/Ex4_3.sce b/728/CH4/EX4.3/Ex4_3.sce
new file mode 100755
index 000000000..6c2514717
--- /dev/null
+++ b/728/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,12 @@
+//Caption:Calculate the breakdown power of air filled coaxial cable at 9.375 GHz.
+//Exa:4.3
+clc;
+clear;
+close;
+//Given:
+a=2.42;//in cm
+x=2.3;//x=(b/a)
+P_bd=3600*a^2*log(x);//in kilowatts
+disp(P_bd,'Breakdown Power (in kW) =');
+
+//answer in book is wrongly written as 398 kW. \ No newline at end of file
diff --git a/728/CH4/EX4.3/Ex4_3.txt b/728/CH4/EX4.3/Ex4_3.txt
new file mode 100755
index 000000000..b28873e45
--- /dev/null
+++ b/728/CH4/EX4.3/Ex4_3.txt
@@ -0,0 +1,12 @@
+//Caption:Calculate the breakdown power of air filled coaxial cable at 9.375 GHz.
+//Exa:4.3
+clc;
+clear;
+close;
+//Given:
+a=2.42;//in cm
+x=2.3;//x=(b/a)
+P_bd=3600*a^2*log(x);//in kilowatts
+disp(P_bd,'Breakdown Power (in kW) =');
+
+//answer in book is wrongly written as 398 kW.
diff --git a/728/CH4/EX4.4/Ex4_4.sce b/728/CH4/EX4.4/Ex4_4.sce
new file mode 100755
index 000000000..93768e4df
--- /dev/null
+++ b/728/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,13 @@
+//Caption:Calculate charcteristic impedance & velocity of propagation.
+//Exa:4.4
+clc;
+clear;
+close;
+b=0.3175;//in cm
+d=0.0539;//in cm
+c=3*10^8;//in m/s
+e_r=2.32;
+Z_o=60*log(4*b/(%pi*d))/sqrt(e_r);//in ohms
+V_p=c/sqrt(e_r);//in m/s
+disp(Z_o,'Charcteristic impedance (in ohms) =');
+disp(V_p,'Velocity of propagation (in m/s) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.4/Ex4_4.txt b/728/CH4/EX4.4/Ex4_4.txt
new file mode 100755
index 000000000..1f3fe3240
--- /dev/null
+++ b/728/CH4/EX4.4/Ex4_4.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate charcteristic impedance & velocity of propagation.
+//Exa:4.4
+clc;
+clear;
+close;
+b=0.3175;//in cm
+d=0.0539;//in cm
+c=3*10^8;//in m/s
+e_r=2.32;
+Z_o=60*log(4*b/(%pi*d))/sqrt(e_r);//in ohms
+V_p=c/sqrt(e_r);//in m/s
+disp(Z_o,'Charcteristic impedance (in ohms) =');
+disp(V_p,'Velocity of propagation (in m/s) =');
diff --git a/728/CH4/EX4.4/Ex4_4_ans.txt b/728/CH4/EX4.4/Ex4_4_ans.txt
new file mode 100755
index 000000000..e5dead8db
--- /dev/null
+++ b/728/CH4/EX4.4/Ex4_4_ans.txt
@@ -0,0 +1,9 @@
+
+ Charcteristic impedance (in ohms) =
+
+ 79.371265
+
+ Velocity of propagation (in m/s) =
+
+ 1.970D+08
+
diff --git a/728/CH4/EX4.5/Ex4_5.sce b/728/CH4/EX4.5/Ex4_5.sce
new file mode 100755
index 000000000..0a6349b86
--- /dev/null
+++ b/728/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,25 @@
+//Caption:Calculate charcteristic impedance & effective dielectric constant & velocity of propagation
+//Exa:4.5
+clc;
+clear;
+close;
+e_r=9.7;
+c=3*10^8;//in m/s
+r_1=0.5;//when ratio: (W/h)=0.5
+r_2=5;//when ratio: (W/h)=5
+//For W/h ratio=0.5
+e_eff_1=(e_r+1)/2+((e_r-1)/2)*[1/{sqrt(1+12*(1/r_1))+0.04*(1-r_1)}];
+Z_o_1=60*log(8/r_1+r_1/4)/sqrt(e_eff_1);
+v_1=c/sqrt(e_eff_1);
+disp("For W/h=0.5 ,");
+disp(e_eff_1,'Effective dielectric constant =');
+disp(Z_o_1,'Charcteristic impedance (in ohms) =');
+disp(v_1,'Velocity of propagation (in m/s) =');
+//For W/h ratio=5
+e_eff_2=(e_r+1)/2+((e_r-1)/2)*[1/{sqrt(1+12*(1/r_2))}];
+Z_o_2=120*%pi*[1/{r_2+1.393+0.667*log(1.444+r_2)}]/sqrt(e_eff_2);
+v_2=c/sqrt(e_eff_2);
+disp("For W/h=5,");
+disp(e_eff_2,'Effective dielectric constant =');
+disp(Z_o_2,'Charcteristic impedance (in ohms) =');
+disp(v_2,'Velocity of propagation (in m/s) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.5/Ex4_5.txt b/728/CH4/EX4.5/Ex4_5.txt
new file mode 100755
index 000000000..b4c72681e
--- /dev/null
+++ b/728/CH4/EX4.5/Ex4_5.txt
@@ -0,0 +1,25 @@
+//Caption:Calculate charcteristic impedance & effective dielectric constant & velocity of propagation
+//Exa:4.5
+clc;
+clear;
+close;
+e_r=9.7;
+c=3*10^8;//in m/s
+r_1=0.5;//when ratio: (W/h)=0.5
+r_2=5;//when ratio: (W/h)=5
+//For W/h ratio=0.5
+e_eff_1=(e_r+1)/2+((e_r-1)/2)*[1/{sqrt(1+12*(1/r_1))+0.04*(1-r_1)}];
+Z_o_1=60*log(8/r_1+r_1/4)/sqrt(e_eff_1);
+v_1=c/sqrt(e_eff_1);
+disp("For W/h=0.5 ,");
+disp(e_eff_1,'Effective dielectric constant =');
+disp(Z_o_1,'Charcteristic impedance (in ohms) =');
+disp(v_1,'Velocity of propagation (in m/s) =');
+//For W/h ratio=5
+e_eff_2=(e_r+1)/2+((e_r-1)/2)*[1/{sqrt(1+12*(1/r_2))}];
+Z_o_2=120*%pi*[1/{r_2+1.393+0.667*log(1.444+r_2)}]/sqrt(e_eff_2);
+v_2=c/sqrt(e_eff_2);
+disp("For W/h=5,");
+disp(e_eff_2,'Effective dielectric constant =');
+disp(Z_o_2,'Charcteristic impedance (in ohms) =');
+disp(v_2,'Velocity of propagation (in m/s) =');
diff --git a/728/CH4/EX4.5/Ex4_5_ans.txt b/728/CH4/EX4.5/Ex4_5_ans.txt
new file mode 100755
index 000000000..1096bd3fe
--- /dev/null
+++ b/728/CH4/EX4.5/Ex4_5_ans.txt
@@ -0,0 +1,28 @@
+
+ For W/h=0.5 ,
+
+ Effective dielectric constant =
+
+ 6.2165339
+
+ Charcteristic impedance (in ohms) =
+
+ 66.908274
+
+ Velocity of propagation (in m/s) =
+
+ 1.203D+08
+
+ For W/h=5,
+
+ Effective dielectric constant =
+
+ 7.7091187
+
+ Charcteristic impedance (in ohms) =
+
+ 17.781923
+
+ Velocity of propagation (in m/s) =
+
+ 1.080D+08
diff --git a/728/CH4/EX4.6/Ex4_6.sce b/728/CH4/EX4.6/Ex4_6.sce
new file mode 100755
index 000000000..3bd2335e4
--- /dev/null
+++ b/728/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,24 @@
+//Caption:Calculate ratio of circular waveguide cross-sectional area to rectangular waveguide cross-section
+//Exa:4.6
+clc;
+clear;
+close;
+//For TE Wave propagated:
+//for Rectangular , taking (a=2b)
+r=100;//assume
+//for TE11, wavelength=2*pi*r/1.841
+//for TE10, wavelength=2a
+a=(2*%pi*r/1.841)/2;
+ar_rec_TE=(a)*(a/2);
+ar_cir_TE=%pi*r^2;
+ratio_TE=(ar_cir_TE)/(ar_rec_TE);
+disp(ratio_TE,'Ratio of Circular & Rectangular coss-section area (in TE) =');
+//For TM Wave propagated:
+//for Rectangular , taking (a=2b)
+//for TE01, wavelength=2.6155*r
+//for TE11, wavelength=4b/sqrt(5)
+b=(2.6155*r)/1.78885;
+ar_rec_TM=(b)*(b);
+ar_cir_TM=%pi*r^2;
+ratio_TM=(ar_cir_TM)/(ar_rec_TM);
+disp(ratio_TM,'Ratio of Circular & Rectangular coss-section area (in TM) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.6/Ex4_6.txt b/728/CH4/EX4.6/Ex4_6.txt
new file mode 100755
index 000000000..8ee29199d
--- /dev/null
+++ b/728/CH4/EX4.6/Ex4_6.txt
@@ -0,0 +1,24 @@
+//Caption:Calculate ratio of circular waveguide cross-sectional area to rectangular waveguide cross-section
+//Exa:4.6
+clc;
+clear;
+close;
+//For TE Wave propagated:
+//for Rectangular , taking (a=2b)
+r=100;//assume
+//for TE11, wavelength=2*pi*r/1.841
+//for TE10, wavelength=2a
+a=(2*%pi*r/1.841)/2;
+ar_rec_TE=(a)*(a/2);
+ar_cir_TE=%pi*r^2;
+ratio_TE=(ar_cir_TE)/(ar_rec_TE);
+disp(ratio_TE,'Ratio of Circular & Rectangular coss-section area (in TE) =');
+//For TM Wave propagated:
+//for Rectangular , taking (a=2b)
+//for TE01, wavelength=2.6155*r
+//for TE11, wavelength=4b/sqrt(5)
+b=(2.6155*r)/1.78885;
+ar_rec_TM=(b)*(b);
+ar_cir_TM=%pi*r^2;
+ratio_TM=(ar_cir_TM)/(ar_rec_TM);
+disp(ratio_TM,'Ratio of Circular & Rectangular coss-section area (in TM) =');
diff --git a/728/CH4/EX4.6/Ex4_6_ans.txt b/728/CH4/EX4.6/Ex4_6_ans.txt
new file mode 100755
index 000000000..88c57e7fb
--- /dev/null
+++ b/728/CH4/EX4.6/Ex4_6_ans.txt
@@ -0,0 +1,8 @@
+ Ratio of Circular & Rectangular coss-section area (in TE) =
+
+ 2.1576833
+
+ Ratio of Circular & Rectangular coss-section area (in TM) =
+
+ 1.4695632
+
diff --git a/728/CH4/EX4.7/Ex4_7.sce b/728/CH4/EX4.7/Ex4_7.sce
new file mode 100755
index 000000000..2e42445ce
--- /dev/null
+++ b/728/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,12 @@
+//Caption:Calculate breadth of rectangular waveguide
+//Exa:4.7
+clc;
+clear;
+close;
+f=9*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_g=4;//in m
+wl_o=c/f;
+wl_c=[sqrt(1-((wl_o/wl_g)^2))/wl_o]^-1;
+b=wl_c/4;
+disp(b,'Breadth of rectangular waveguide (in cm) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.7/Ex4_7.txt b/728/CH4/EX4.7/Ex4_7.txt
new file mode 100755
index 000000000..e35dca3d8
--- /dev/null
+++ b/728/CH4/EX4.7/Ex4_7.txt
@@ -0,0 +1,12 @@
+//Caption:Calculate breadth of rectangular waveguide
+//Exa:4.7
+clc;
+clear;
+close;
+f=9*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_g=4;//in m
+wl_o=c/f;
+wl_c=[sqrt(1-((wl_o/wl_g)^2))/wl_o]^-1;
+b=wl_c/4;
+disp(b,'Breadth of rectangular waveguide (in cm) =');
diff --git a/728/CH4/EX4.7/Ex4_7_ans.txt b/728/CH4/EX4.7/Ex4_7_ans.txt
new file mode 100755
index 000000000..6a17cc42a
--- /dev/null
+++ b/728/CH4/EX4.7/Ex4_7_ans.txt
@@ -0,0 +1,4 @@
+
+ Breadth of rectangular waveguide (in cm) =
+
+ 1.5075567
diff --git a/728/CH4/EX4.8/Ex4_8.sce b/728/CH4/EX4.8/Ex4_8.sce
new file mode 100755
index 000000000..c831d77c3
--- /dev/null
+++ b/728/CH4/EX4.8/Ex4_8.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate the cutoff wavelength, guide wavelength, group & phase velocities
+//Exa:4.8
+clc;
+clear;
+close;
+a=10;//in cm
+c=3*10^10;//in cm/s
+wl_c=2*a;//in cm
+f=2.5*10^9;//in Hz
+wl_o=c/f;
+wl_g=wl_o/(sqrt(1-(wl_o/wl_c)^2));//in cm
+V_p=c/(sqrt(1-(wl_o/wl_c)^2));
+V_g=c^2/V_p;
+disp(wl_c,'Cut-off wavelength (in cm) =');
+disp(wl_g,'Guide wavelength (in cm) =');
+disp(V_p,'Phase velocity (in cm/s) =');
+disp(V_g,'Group velocity (in cm/s) ='); \ No newline at end of file
diff --git a/728/CH4/EX4.8/Ex4_8.txt b/728/CH4/EX4.8/Ex4_8.txt
new file mode 100755
index 000000000..84ac1eeac
--- /dev/null
+++ b/728/CH4/EX4.8/Ex4_8.txt
@@ -0,0 +1,17 @@
+//Caption:Calculate the cutoff wavelength, guide wavelength, group & phase velocities
+//Exa:4.8
+clc;
+clear;
+close;
+a=10;//in cm
+c=3*10^10;//in cm/s
+wl_c=2*a;//in cm
+f=2.5*10^9;//in Hz
+wl_o=c/f;
+wl_g=wl_o/(sqrt(1-(wl_o/wl_c)^2));//in cm
+V_p=c/(sqrt(1-(wl_o/wl_c)^2));
+V_g=c^2/V_p;
+disp(wl_c,'Cut-off wavelength (in cm) =');
+disp(wl_g,'Guide wavelength (in cm) =');
+disp(V_p,'Phase velocity (in cm/s) =');
+disp(V_g,'Group velocity (in cm/s) =');
diff --git a/728/CH4/EX4.8/Ex4_8_ans.txt b/728/CH4/EX4.8/Ex4_8_ans.txt
new file mode 100755
index 000000000..02d1dd212
--- /dev/null
+++ b/728/CH4/EX4.8/Ex4_8_ans.txt
@@ -0,0 +1,15 @@
+ Cut-off wavelength (in cm) =
+
+ 20.
+
+ Guide wavelength (in cm) =
+
+ 15.
+
+ Phase velocity (in cm/s) =
+
+ 3.750D+10
+
+ Group velocity (in cm/s) =
+
+ 2.400D+10
diff --git a/728/CH4/EX4.9/EX4_9_ans.txt b/728/CH4/EX4.9/EX4_9_ans.txt
new file mode 100755
index 000000000..cec5ea754
--- /dev/null
+++ b/728/CH4/EX4.9/EX4_9_ans.txt
@@ -0,0 +1,15 @@
+Only TE10 mode is possible
+
+ Cut-off frequency(in Hz) =
+
+ 6.000D+09
+
+ Guide wavelength for TE10 (in cm) =
+
+ 4.869206
+
+ TM11 also propagates
+
+ Guide wavelength for TM11 (in cm) =
+
+ 4.869206
diff --git a/728/CH4/EX4.9/Ex4_9.sce b/728/CH4/EX4.9/Ex4_9.sce
new file mode 100755
index 000000000..f28c69383
--- /dev/null
+++ b/728/CH4/EX4.9/Ex4_9.sce
@@ -0,0 +1,24 @@
+//Caption:Calculate (i)-possible modes,(ii)-cut-off frequencies,(iii)-guide wavelength
+//Exa:4.9
+clc;
+clear;
+close;
+//For TE mode:
+a=2.5;//in cm
+b=1;//in cm
+f=8.6*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+wl_c_1=2*b;//for TE01
+wl_c_2=2*a;//for TE10
+disp( 'Only TE10 mode is possible');
+f_c=c/wl_c_2;
+wl_c_3=2*a*b/sqrt(a^2+b^2);//for TE11 & TM11
+wl_g_TE10=wl_o/(sqrt(1-(wl_o/wl_c_2)^2));//for TE10
+disp(f_c,'Cut-off frequency(in Hz) =');
+disp(wl_g_TE10,'Guide wavelength for TE10 (in cm) =');
+//For TM mode:
+disp('TM11 also propagates');
+wl_c_TM11=wl_c_3;
+wl_g_TM11=wl_o/(sqrt(1-(wl_o/wl_c_2)^2));//for TM11
+disp(wl_g_TM11,'Guide wavelength for TM11 (in cm) =');
diff --git a/728/CH4/EX4.9/Ex4_9.txt b/728/CH4/EX4.9/Ex4_9.txt
new file mode 100755
index 000000000..f28c69383
--- /dev/null
+++ b/728/CH4/EX4.9/Ex4_9.txt
@@ -0,0 +1,24 @@
+//Caption:Calculate (i)-possible modes,(ii)-cut-off frequencies,(iii)-guide wavelength
+//Exa:4.9
+clc;
+clear;
+close;
+//For TE mode:
+a=2.5;//in cm
+b=1;//in cm
+f=8.6*10^9;//in Hz
+c=3*10^10;//in cm/s
+wl_o=c/f;
+wl_c_1=2*b;//for TE01
+wl_c_2=2*a;//for TE10
+disp( 'Only TE10 mode is possible');
+f_c=c/wl_c_2;
+wl_c_3=2*a*b/sqrt(a^2+b^2);//for TE11 & TM11
+wl_g_TE10=wl_o/(sqrt(1-(wl_o/wl_c_2)^2));//for TE10
+disp(f_c,'Cut-off frequency(in Hz) =');
+disp(wl_g_TE10,'Guide wavelength for TE10 (in cm) =');
+//For TM mode:
+disp('TM11 also propagates');
+wl_c_TM11=wl_c_3;
+wl_g_TM11=wl_o/(sqrt(1-(wl_o/wl_c_2)^2));//for TM11
+disp(wl_g_TM11,'Guide wavelength for TM11 (in cm) =');
diff --git a/728/CH5/EX5.1/Ex5_1.sce b/728/CH5/EX5.1/Ex5_1.sce
new file mode 100755
index 000000000..fd11d4075
--- /dev/null
+++ b/728/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,12 @@
+//Caption:Determine the minimum distance between two end plates
+//Exa:5.1
+clc;
+clear;
+close;
+//Given:
+a=3;//in cm
+c=3*10^10;//in cm/s
+f=10*10^9;//in Hz
+P_01=2.405;
+d=%pi/sqrt(f^2*4*%pi^2/c^2-(P_01/a)^2);
+disp(d,'Minimum distance (in cm) ='); \ No newline at end of file
diff --git a/728/CH5/EX5.1/Ex5_1.txt b/728/CH5/EX5.1/Ex5_1.txt
new file mode 100755
index 000000000..fd11d4075
--- /dev/null
+++ b/728/CH5/EX5.1/Ex5_1.txt
@@ -0,0 +1,12 @@
+//Caption:Determine the minimum distance between two end plates
+//Exa:5.1
+clc;
+clear;
+close;
+//Given:
+a=3;//in cm
+c=3*10^10;//in cm/s
+f=10*10^9;//in Hz
+P_01=2.405;
+d=%pi/sqrt(f^2*4*%pi^2/c^2-(P_01/a)^2);
+disp(d,'Minimum distance (in cm) ='); \ No newline at end of file
diff --git a/728/CH5/EX5.1/Ex5_1_ans.txt b/728/CH5/EX5.1/Ex5_1_ans.txt
new file mode 100755
index 000000000..c633d250e
--- /dev/null
+++ b/728/CH5/EX5.1/Ex5_1_ans.txt
@@ -0,0 +1,3 @@
+ Minimum distance (in cm) =
+
+ 1.6236496 \ No newline at end of file
diff --git a/728/CH5/EX5.2/EX5_2.sce b/728/CH5/EX5.2/EX5_2.sce
new file mode 100755
index 000000000..47048b00b
--- /dev/null
+++ b/728/CH5/EX5.2/EX5_2.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate the lowest frequency of a rectangular cavity resonator
+//Exa:5.2
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=2;//in cm
+b=1;//in cm
+d=3; //in cm
+disp('Dominant mode is TE101');
+m=1;
+n=0;
+p=1;
+f=(c/2)*[(m/a)^2+(n/b)^2+(p/d)^2]^0.5;
+disp(f/10^9,'Lowest resonant frequency(in GHz) =');
diff --git a/728/CH5/EX5.2/Ex5_2.txt b/728/CH5/EX5.2/Ex5_2.txt
new file mode 100755
index 000000000..47048b00b
--- /dev/null
+++ b/728/CH5/EX5.2/Ex5_2.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate the lowest frequency of a rectangular cavity resonator
+//Exa:5.2
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=2;//in cm
+b=1;//in cm
+d=3; //in cm
+disp('Dominant mode is TE101');
+m=1;
+n=0;
+p=1;
+f=(c/2)*[(m/a)^2+(n/b)^2+(p/d)^2]^0.5;
+disp(f/10^9,'Lowest resonant frequency(in GHz) =');
diff --git a/728/CH5/EX5.2/Ex5_2_ans.txt b/728/CH5/EX5.2/Ex5_2_ans.txt
new file mode 100755
index 000000000..f54fdad4d
--- /dev/null
+++ b/728/CH5/EX5.2/Ex5_2_ans.txt
@@ -0,0 +1,6 @@
+Dominant mode is TE101
+
+ Lowest resonant frequency(in GHz) =
+
+ 9.0138782
+ \ No newline at end of file
diff --git a/728/CH5/EX5.3/Ex5_3.sce b/728/CH5/EX5.3/Ex5_3.sce
new file mode 100755
index 000000000..22e23aa84
--- /dev/null
+++ b/728/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,19 @@
+//Caption:Calculate the resonant frequency of a circular cavity resonator
+//Exa:5.3
+clc;
+clear;
+close;
+//Given:
+d=12.5;//diameter(in cm)
+c=3*10^10;//in cm/s
+l=5;//length(in cm)
+a=d/2;
+//For TM012 mode:
+n=0;
+m=1;
+p=2;
+P=2.405;
+f=(c/(2*%pi))*[(P/a)^2+(p*%pi/d)^2]^0.5;
+disp(f/10^9,'Resonant frequency (in GHz) =');
+
+//Answer in book in wrongly given as 6.27GHz \ No newline at end of file
diff --git a/728/CH5/EX5.3/Ex5_3.txt b/728/CH5/EX5.3/Ex5_3.txt
new file mode 100755
index 000000000..22e23aa84
--- /dev/null
+++ b/728/CH5/EX5.3/Ex5_3.txt
@@ -0,0 +1,19 @@
+//Caption:Calculate the resonant frequency of a circular cavity resonator
+//Exa:5.3
+clc;
+clear;
+close;
+//Given:
+d=12.5;//diameter(in cm)
+c=3*10^10;//in cm/s
+l=5;//length(in cm)
+a=d/2;
+//For TM012 mode:
+n=0;
+m=1;
+p=2;
+P=2.405;
+f=(c/(2*%pi))*[(P/a)^2+(p*%pi/d)^2]^0.5;
+disp(f/10^9,'Resonant frequency (in GHz) =');
+
+//Answer in book in wrongly given as 6.27GHz \ No newline at end of file
diff --git a/728/CH5/EX5.3/Ex5_3_ans.txt b/728/CH5/EX5.3/Ex5_3_ans.txt
new file mode 100755
index 000000000..b97abc2ee
--- /dev/null
+++ b/728/CH5/EX5.3/Ex5_3_ans.txt
@@ -0,0 +1,4 @@
+ Resonant frequency (in GHz) =
+
+ 3.022518
+ \ No newline at end of file
diff --git a/728/CH5/EX5.4/Ex5_4.sce b/728/CH5/EX5.4/Ex5_4.sce
new file mode 100755
index 000000000..c96ff8a1c
--- /dev/null
+++ b/728/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate the resonant frequency of a circular cavity resonator
+//Exa:5.4
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=3;//in cm
+b=2;//in cm
+d=4;//in cm
+//For TE101:
+m=1;
+n=0;
+p=1;
+f=(c/2)*[(m/a)^2+(n/b)^2+(p/d)^2]^0.5;
+disp(f/10^9,'Resonant frequency(in GHz) =');
diff --git a/728/CH5/EX5.4/Ex5_4.txt b/728/CH5/EX5.4/Ex5_4.txt
new file mode 100755
index 000000000..c96ff8a1c
--- /dev/null
+++ b/728/CH5/EX5.4/Ex5_4.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate the resonant frequency of a circular cavity resonator
+//Exa:5.4
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=3;//in cm
+b=2;//in cm
+d=4;//in cm
+//For TE101:
+m=1;
+n=0;
+p=1;
+f=(c/2)*[(m/a)^2+(n/b)^2+(p/d)^2]^0.5;
+disp(f/10^9,'Resonant frequency(in GHz) =');
diff --git a/728/CH5/EX5.4/Ex5_4_ans.txt b/728/CH5/EX5.4/Ex5_4_ans.txt
new file mode 100755
index 000000000..deee42067
--- /dev/null
+++ b/728/CH5/EX5.4/Ex5_4_ans.txt
@@ -0,0 +1,4 @@
+
+ Resonant frequency(in GHz) =
+
+ 6.25 \ No newline at end of file
diff --git a/728/CH6/EX6.10/Ex6_10.sce b/728/CH6/EX6.10/Ex6_10.sce
new file mode 100755
index 000000000..3c456d4e7
--- /dev/null
+++ b/728/CH6/EX6.10/Ex6_10.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.10
+clc;
+clear;
+close;
+In_loss=0.5;//insertion loss(in dB)
+C=20;//in dB
+D=35;//in dB
+Pi_Pf=10^(C/10);
+Pi=90;//in Watts
+Pf=Pi/Pi_Pf;
+Pf_Pb=10^(D/10);
+Pb=Pf/Pf_Pb;
+P_rec=(Pi-Pf-Pb);//Power received (in Watts)
+P_rec_dB=10*log(Pi/P_rec)/log(10);
+P_rec_eff=P_rec_dB-In_loss;//Effective power received (in dB)
+disp(P_rec_eff,'Effective power received (in dB)='); \ No newline at end of file
diff --git a/728/CH6/EX6.10/Ex6_10.txt b/728/CH6/EX6.10/Ex6_10.txt
new file mode 100755
index 000000000..3c456d4e7
--- /dev/null
+++ b/728/CH6/EX6.10/Ex6_10.txt
@@ -0,0 +1,17 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.10
+clc;
+clear;
+close;
+In_loss=0.5;//insertion loss(in dB)
+C=20;//in dB
+D=35;//in dB
+Pi_Pf=10^(C/10);
+Pi=90;//in Watts
+Pf=Pi/Pi_Pf;
+Pf_Pb=10^(D/10);
+Pb=Pf/Pf_Pb;
+P_rec=(Pi-Pf-Pb);//Power received (in Watts)
+P_rec_dB=10*log(Pi/P_rec)/log(10);
+P_rec_eff=P_rec_dB-In_loss;//Effective power received (in dB)
+disp(P_rec_eff,'Effective power received (in dB)='); \ No newline at end of file
diff --git a/728/CH6/EX6.10/Ex6_10_ans.txt b/728/CH6/EX6.10/Ex6_10_ans.txt
new file mode 100755
index 000000000..e524f64d7
--- /dev/null
+++ b/728/CH6/EX6.10/Ex6_10_ans.txt
@@ -0,0 +1,5 @@
+
+ Effective power received (in dB)=
+
+ - 0.4563381
+ \ No newline at end of file
diff --git a/728/CH6/EX6.11/Ex6_11.sce b/728/CH6/EX6.11/Ex6_11.sce
new file mode 100755
index 000000000..9b9bc1b54
--- /dev/null
+++ b/728/CH6/EX6.11/Ex6_11.sce
@@ -0,0 +1,13 @@
+ //Caption:Calculate (i)-directivity ,(ii)-coupling, (iii)-isolation
+//Exa:6.11
+clc;
+clear;
+close;
+S_13=0.1;
+S_14=0.05;
+C=-20*log(S_13)/log(10);
+D=20*log(S_13/S_14)/log(10);
+I=C+D;
+disp(C,'Coupling (in dB) =');
+disp(D,'Directivity (in dB)) =');
+disp(I,'Isolation (in dB) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.11/Ex6_11.txt b/728/CH6/EX6.11/Ex6_11.txt
new file mode 100755
index 000000000..9b9bc1b54
--- /dev/null
+++ b/728/CH6/EX6.11/Ex6_11.txt
@@ -0,0 +1,13 @@
+ //Caption:Calculate (i)-directivity ,(ii)-coupling, (iii)-isolation
+//Exa:6.11
+clc;
+clear;
+close;
+S_13=0.1;
+S_14=0.05;
+C=-20*log(S_13)/log(10);
+D=20*log(S_13/S_14)/log(10);
+I=C+D;
+disp(C,'Coupling (in dB) =');
+disp(D,'Directivity (in dB)) =');
+disp(I,'Isolation (in dB) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.11/Ex6_11_ans.txt b/728/CH6/EX6.11/Ex6_11_ans.txt
new file mode 100755
index 000000000..6984050e7
--- /dev/null
+++ b/728/CH6/EX6.11/Ex6_11_ans.txt
@@ -0,0 +1,12 @@
+Coupling (in dB) =
+
+ 20.
+
+ Directivity (in dB)) =
+
+ 6.0205999
+
+ Isolation (in dB) =
+
+ 26.0206
+ \ No newline at end of file
diff --git a/728/CH6/EX6.12/Ex6_12.sce b/728/CH6/EX6.12/Ex6_12.sce
new file mode 100755
index 000000000..d7a15ef1e
--- /dev/null
+++ b/728/CH6/EX6.12/Ex6_12.sce
@@ -0,0 +1,10 @@
+//Caption:Calculate the value of VSWR
+//Exa:6.12
+clc;
+clear;
+close;
+D=3.5;//distance of seperation(in cm)
+w_l=2*D;//wavelength
+d2_d1=2.5;//d2-d1(in m)
+S=w_l/(%pi*d2_d1*10^-1);
+disp(S,'VSWR ='); \ No newline at end of file
diff --git a/728/CH6/EX6.12/Ex6_12.txt b/728/CH6/EX6.12/Ex6_12.txt
new file mode 100755
index 000000000..d7a15ef1e
--- /dev/null
+++ b/728/CH6/EX6.12/Ex6_12.txt
@@ -0,0 +1,10 @@
+//Caption:Calculate the value of VSWR
+//Exa:6.12
+clc;
+clear;
+close;
+D=3.5;//distance of seperation(in cm)
+w_l=2*D;//wavelength
+d2_d1=2.5;//d2-d1(in m)
+S=w_l/(%pi*d2_d1*10^-1);
+disp(S,'VSWR ='); \ No newline at end of file
diff --git a/728/CH6/EX6.12/Ex6_12_ans.txt b/728/CH6/EX6.12/Ex6_12_ans.txt
new file mode 100755
index 000000000..8707fb0ff
--- /dev/null
+++ b/728/CH6/EX6.12/Ex6_12_ans.txt
@@ -0,0 +1,5 @@
+
+ VSWR =
+
+ 8.9126768
+ \ No newline at end of file
diff --git a/728/CH6/EX6.13/Ex6_13.sce b/728/CH6/EX6.13/Ex6_13.sce
new file mode 100755
index 000000000..6669f4a3c
--- /dev/null
+++ b/728/CH6/EX6.13/Ex6_13.sce
@@ -0,0 +1,9 @@
+//Caption:Calculate the phase shift of the component
+//Exa:6.13
+clc;
+clear;
+close;
+w_l=7.2;//wavelength (in cm)
+x=10.5-9.3;
+Phase_shift=(2*%pi*x)/(w_l);
+disp(Phase_shift*180/%pi,'Phase Shift (in degree) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.13/Ex6_13.txt b/728/CH6/EX6.13/Ex6_13.txt
new file mode 100755
index 000000000..6669f4a3c
--- /dev/null
+++ b/728/CH6/EX6.13/Ex6_13.txt
@@ -0,0 +1,9 @@
+//Caption:Calculate the phase shift of the component
+//Exa:6.13
+clc;
+clear;
+close;
+w_l=7.2;//wavelength (in cm)
+x=10.5-9.3;
+Phase_shift=(2*%pi*x)/(w_l);
+disp(Phase_shift*180/%pi,'Phase Shift (in degree) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.13/Ex6_13_ans.txt b/728/CH6/EX6.13/Ex6_13_ans.txt
new file mode 100755
index 000000000..520cc2b3b
--- /dev/null
+++ b/728/CH6/EX6.13/Ex6_13_ans.txt
@@ -0,0 +1,5 @@
+
+ Phase Shift (in degree) =
+
+ 60.
+ \ No newline at end of file
diff --git a/728/CH6/EX6.2/Ex6_2.sce b/728/CH6/EX6.2/Ex6_2.sce
new file mode 100755
index 000000000..1fe86dabd
--- /dev/null
+++ b/728/CH6/EX6.2/Ex6_2.sce
@@ -0,0 +1,13 @@
+//Caption:Find the distance that the position of port 1 should be shifted .
+//Exa:6.2
+clc;
+clear;
+close;
+Beeta=34.3;//in rad/m
+// S=[0,0.5*%e^(%i*53.13);0.5*%e^(%i*53.13),0];
+// S'=[0,0.5*%e^(%i*53.13-x);0.5*%e^(%i*53.13-x),0];
+//For S12& S21 to be real ,
+x=53.5;//in degrees
+x_rad=53.5*%pi/180;
+l=x_rad/Beeta;
+disp(l*100,'Distance (in cm)='); \ No newline at end of file
diff --git a/728/CH6/EX6.2/Ex6_2.txt b/728/CH6/EX6.2/Ex6_2.txt
new file mode 100755
index 000000000..587cdb763
--- /dev/null
+++ b/728/CH6/EX6.2/Ex6_2.txt
@@ -0,0 +1,13 @@
+//Caption:Find the distance that the position of port 1 should be shifted .
+//Exa:6.2
+clc;
+clear;
+close;
+Beeta=34.3;//in rad/m
+// S=[0,0.5*%e^(%i*53.13);0.5*%e^(%i*53.13),0];
+// S'=[0,0.5*%e^(%i*53.13-x);0.5*%e^(%i*53.13-x),0];
+//For S12& S21 to be real ,
+x=53.5;//in degrees
+x_rad=53.5*%pi/180;
+l=x_rad/Beeta;
+disp(l*100,'Distance (in cm)='); \ No newline at end of file
diff --git a/728/CH6/EX6.2/Ex6_2_ans.txt b/728/CH6/EX6.2/Ex6_2_ans.txt
new file mode 100755
index 000000000..448c41aaa
--- /dev/null
+++ b/728/CH6/EX6.2/Ex6_2_ans.txt
@@ -0,0 +1,5 @@
+
+ Distance (in cm)=
+
+ 2.7223066
+ \ No newline at end of file
diff --git a/728/CH6/EX6.3/Ex6_3.sce b/728/CH6/EX6.3/Ex6_3.sce
new file mode 100755
index 000000000..a3fd44435
--- /dev/null
+++ b/728/CH6/EX6.3/Ex6_3.sce
@@ -0,0 +1,28 @@
+//Caption:Determine the scattering parameters for 10 dB direction coupler
+//Exa:6.3
+clc;
+clear;
+close;
+D=30;//in dB
+VSWR=1;
+C=10;
+//p1_p4 = p1/p4
+p1_p4=10^(C/-10);
+S_41=sqrt(p1_p4);
+S_14=S_41;//As matched & lossless
+S_31=S_41^2/10^(D/10);
+S_11=(VSWR-1)/(VSWR+1);
+S_22=S_11;
+S_44=S_11;
+S_33=S_11;
+S_21=sqrt(1-0.1-10^-4);
+S_12=S_21;
+S_34=sqrt(1-0.1-10^-4);
+S_43=S_34;
+S_24=sqrt(1-0.1-S_34^2);
+S_42=S_24;
+S_23=S_41;
+S_32=S_23;
+S_13=S_31;
+S=[S_11,S_12,S_13,S_14;S_21,S_22,S_23,S_24;S_31,S_32,S_33,S_34;S_41,S_42,S_43,S_44];
+disp(S,'Required Scattering Parameters are'); \ No newline at end of file
diff --git a/728/CH6/EX6.3/Ex6_3.txt b/728/CH6/EX6.3/Ex6_3.txt
new file mode 100755
index 000000000..a3fd44435
--- /dev/null
+++ b/728/CH6/EX6.3/Ex6_3.txt
@@ -0,0 +1,28 @@
+//Caption:Determine the scattering parameters for 10 dB direction coupler
+//Exa:6.3
+clc;
+clear;
+close;
+D=30;//in dB
+VSWR=1;
+C=10;
+//p1_p4 = p1/p4
+p1_p4=10^(C/-10);
+S_41=sqrt(p1_p4);
+S_14=S_41;//As matched & lossless
+S_31=S_41^2/10^(D/10);
+S_11=(VSWR-1)/(VSWR+1);
+S_22=S_11;
+S_44=S_11;
+S_33=S_11;
+S_21=sqrt(1-0.1-10^-4);
+S_12=S_21;
+S_34=sqrt(1-0.1-10^-4);
+S_43=S_34;
+S_24=sqrt(1-0.1-S_34^2);
+S_42=S_24;
+S_23=S_41;
+S_32=S_23;
+S_13=S_31;
+S=[S_11,S_12,S_13,S_14;S_21,S_22,S_23,S_24;S_31,S_32,S_33,S_34;S_41,S_42,S_43,S_44];
+disp(S,'Required Scattering Parameters are'); \ No newline at end of file
diff --git a/728/CH6/EX6.3/Ex6_3_ans.txt b/728/CH6/EX6.3/Ex6_3_ans.txt
new file mode 100755
index 000000000..e4344812e
--- /dev/null
+++ b/728/CH6/EX6.3/Ex6_3_ans.txt
@@ -0,0 +1,8 @@
+
+ Required Scattering Parameters are
+
+ 0. 0.9486306 0.0001 0.3162278
+ 0.9486306 0. 0.3162278 0.01
+ 0.0001 0.3162278 0. 0.9486306
+ 0.3162278 0.01 0.9486306 0.
+ \ No newline at end of file
diff --git a/728/CH6/EX6.4/Ex6_4.sce b/728/CH6/EX6.4/Ex6_4.sce
new file mode 100755
index 000000000..ba41420e6
--- /dev/null
+++ b/728/CH6/EX6.4/Ex6_4.sce
@@ -0,0 +1,14 @@
+//Caption:Determine the powers in the remaining ports
+//Exa:6.4
+clc;
+clear;
+close;
+a_2=0;
+a_3=0;
+a_1=32;//in mW
+b_1=(a_1/2^2)+(a_2/-2)+(a_3/sqrt(2));
+b_2=(a_1/(-2)^2)+(a_2/-2)+(a_3/sqrt(2));
+b_3=(a_1/2)+(a_2/sqrt(2))+(a_3/-sqrt(2));
+disp(b_1,'Power at port1(in mW)=');
+disp(b_2,'Power at port2(in mW) =');
+disp(b_3,'Power at port3(in mW) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.4/Ex6_4.txt b/728/CH6/EX6.4/Ex6_4.txt
new file mode 100755
index 000000000..ba41420e6
--- /dev/null
+++ b/728/CH6/EX6.4/Ex6_4.txt
@@ -0,0 +1,14 @@
+//Caption:Determine the powers in the remaining ports
+//Exa:6.4
+clc;
+clear;
+close;
+a_2=0;
+a_3=0;
+a_1=32;//in mW
+b_1=(a_1/2^2)+(a_2/-2)+(a_3/sqrt(2));
+b_2=(a_1/(-2)^2)+(a_2/-2)+(a_3/sqrt(2));
+b_3=(a_1/2)+(a_2/sqrt(2))+(a_3/-sqrt(2));
+disp(b_1,'Power at port1(in mW)=');
+disp(b_2,'Power at port2(in mW) =');
+disp(b_3,'Power at port3(in mW) ='); \ No newline at end of file
diff --git a/728/CH6/EX6.4/Ex6_4_ans.txt b/728/CH6/EX6.4/Ex6_4_ans.txt
new file mode 100755
index 000000000..fb0c83db7
--- /dev/null
+++ b/728/CH6/EX6.4/Ex6_4_ans.txt
@@ -0,0 +1,13 @@
+
+ Power at port1(in mW)=
+
+ 8.
+
+ Power at port2(in mW) =
+
+ 8.
+
+ Power at port3(in mW) =
+
+ 16.
+ \ No newline at end of file
diff --git a/728/CH6/EX6.5/Ex6_5.sce b/728/CH6/EX6.5/Ex6_5.sce
new file mode 100755
index 000000000..444a76bcc
--- /dev/null
+++ b/728/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,13 @@
+//Caption:Determine the powers in the remaining ports
+//Exa:6.5
+clc;
+clear;
+close;
+b_1=20;
+b_2=20;
+p_1=abs((60-50)/(60+50));
+p_2=abs((75-50)/(75+50));
+P_1=b_1*(1-p_1^2)/2;
+P_2=b_2*(1-p_2^2)/2;
+disp(P_1,'Power in port1 (in mW) =');
+disp(P_2,'Power in port2 (in mW) =');
diff --git a/728/CH6/EX6.5/Ex6_5.txt b/728/CH6/EX6.5/Ex6_5.txt
new file mode 100755
index 000000000..444a76bcc
--- /dev/null
+++ b/728/CH6/EX6.5/Ex6_5.txt
@@ -0,0 +1,13 @@
+//Caption:Determine the powers in the remaining ports
+//Exa:6.5
+clc;
+clear;
+close;
+b_1=20;
+b_2=20;
+p_1=abs((60-50)/(60+50));
+p_2=abs((75-50)/(75+50));
+P_1=b_1*(1-p_1^2)/2;
+P_2=b_2*(1-p_2^2)/2;
+disp(P_1,'Power in port1 (in mW) =');
+disp(P_2,'Power in port2 (in mW) =');
diff --git a/728/CH6/EX6.5/Ex6_5_ans.txt b/728/CH6/EX6.5/Ex6_5_ans.txt
new file mode 100755
index 000000000..24de88468
--- /dev/null
+++ b/728/CH6/EX6.5/Ex6_5_ans.txt
@@ -0,0 +1,9 @@
+
+ Power in port1 (in mW) =
+
+ 9.9173554
+
+ Power in port2 (in mW) =
+
+ 9.6
+ \ No newline at end of file
diff --git a/728/CH6/EX6.6/Ex6_6.sce b/728/CH6/EX6.6/Ex6_6.sce
new file mode 100755
index 000000000..94e5b1c20
--- /dev/null
+++ b/728/CH6/EX6.6/Ex6_6.sce
@@ -0,0 +1,20 @@
+//Caption:Determine the powers reflected at port 3 & power divisions at other ports.
+//Exa:6.6
+clc;
+clear;
+close;
+p_1=0.5;
+p_2=0.6;
+p_4=0.8;
+b_1=0.6566;
+b_2=0.7576;
+b_3=0.6536;
+b_4=0.00797;
+a_1=p_1*b_1;
+a_2=p_2*b_2;
+a_3=1;//in Watts
+a_4=p_4*b_4;
+disp(b_1^2,'Power at port 1(in W)=');
+disp(b_2^2,'Power at port 2(in W)=');
+disp(b_3^2,'Power at port 3(in W)=');
+disp(b_4^2,'Power at port 4(in W)='); \ No newline at end of file
diff --git a/728/CH6/EX6.6/Ex6_6.txt b/728/CH6/EX6.6/Ex6_6.txt
new file mode 100755
index 000000000..94e5b1c20
--- /dev/null
+++ b/728/CH6/EX6.6/Ex6_6.txt
@@ -0,0 +1,20 @@
+//Caption:Determine the powers reflected at port 3 & power divisions at other ports.
+//Exa:6.6
+clc;
+clear;
+close;
+p_1=0.5;
+p_2=0.6;
+p_4=0.8;
+b_1=0.6566;
+b_2=0.7576;
+b_3=0.6536;
+b_4=0.00797;
+a_1=p_1*b_1;
+a_2=p_2*b_2;
+a_3=1;//in Watts
+a_4=p_4*b_4;
+disp(b_1^2,'Power at port 1(in W)=');
+disp(b_2^2,'Power at port 2(in W)=');
+disp(b_3^2,'Power at port 3(in W)=');
+disp(b_4^2,'Power at port 4(in W)='); \ No newline at end of file
diff --git a/728/CH6/EX6.6/Ex6_6_ans.txt b/728/CH6/EX6.6/Ex6_6_ans.txt
new file mode 100755
index 000000000..3e205681b
--- /dev/null
+++ b/728/CH6/EX6.6/Ex6_6_ans.txt
@@ -0,0 +1,16 @@
+ Power at port 1(in W)=
+
+ 0.4311236
+
+ Power at port 2(in W)=
+
+ 0.5739578
+
+ Power at port 3(in W)=
+
+ 0.4271930
+
+ Power at port 4(in W)=
+
+ 0.0000635
+ \ No newline at end of file
diff --git a/728/CH6/EX6.7/Ex6_7.sce b/728/CH6/EX6.7/Ex6_7.sce
new file mode 100755
index 000000000..79ab7d265
--- /dev/null
+++ b/728/CH6/EX6.7/Ex6_7.sce
@@ -0,0 +1,13 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.7
+clc;
+clear;
+close;
+In_loss=0.5;//in dB
+S_21=10^(-In_loss/20);
+Isolation=30;//in dB
+S_12=10^(-Isolation/20);
+S_11=0;
+S_22=0;
+S=[S_11,S_12;S_21,S_22];
+disp(S,'Scattering matrix ='); \ No newline at end of file
diff --git a/728/CH6/EX6.7/Ex6_7.txt b/728/CH6/EX6.7/Ex6_7.txt
new file mode 100755
index 000000000..79ab7d265
--- /dev/null
+++ b/728/CH6/EX6.7/Ex6_7.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.7
+clc;
+clear;
+close;
+In_loss=0.5;//in dB
+S_21=10^(-In_loss/20);
+Isolation=30;//in dB
+S_12=10^(-Isolation/20);
+S_11=0;
+S_22=0;
+S=[S_11,S_12;S_21,S_22];
+disp(S,'Scattering matrix ='); \ No newline at end of file
diff --git a/728/CH6/EX6.7/Ex6_7_ans.txt b/728/CH6/EX6.7/Ex6_7_ans.txt
new file mode 100755
index 000000000..e67ea55a5
--- /dev/null
+++ b/728/CH6/EX6.7/Ex6_7_ans.txt
@@ -0,0 +1,6 @@
+
+ Scattering matrix =
+
+ 0. 0.0316228
+ 0.9440609 0.
+ \ No newline at end of file
diff --git a/728/CH6/EX6.9/Ex6_9.sce b/728/CH6/EX6.9/Ex6_9.sce
new file mode 100755
index 000000000..05db74208
--- /dev/null
+++ b/728/CH6/EX6.9/Ex6_9.sce
@@ -0,0 +1,20 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.9
+clc;
+clear;
+close;
+VSWR=1;
+In_loss=0.5;//in dB
+S_21=10^(-In_loss/20);
+Isolation=20;//in dB
+S_12=10^(-Isolation/20);
+S_23=S_12;
+S_31=S_12;
+S_32=S_21;
+S_13=S_21;
+p=(VSWR-1)/(VSWR+1);
+S_11=p;
+S_22=p;
+S_33=p;
+S=[S_11,S_12,S_13;S_21,S_22,S_23;S_31,S_32,S_33];
+disp(S,'Scattering matrix ='); \ No newline at end of file
diff --git a/728/CH6/EX6.9/Ex6_9.txt b/728/CH6/EX6.9/Ex6_9.txt
new file mode 100755
index 000000000..05db74208
--- /dev/null
+++ b/728/CH6/EX6.9/Ex6_9.txt
@@ -0,0 +1,20 @@
+//Caption:Calculate the scattering matrix.
+//Exa:6.9
+clc;
+clear;
+close;
+VSWR=1;
+In_loss=0.5;//in dB
+S_21=10^(-In_loss/20);
+Isolation=20;//in dB
+S_12=10^(-Isolation/20);
+S_23=S_12;
+S_31=S_12;
+S_32=S_21;
+S_13=S_21;
+p=(VSWR-1)/(VSWR+1);
+S_11=p;
+S_22=p;
+S_33=p;
+S=[S_11,S_12,S_13;S_21,S_22,S_23;S_31,S_32,S_33];
+disp(S,'Scattering matrix ='); \ No newline at end of file
diff --git a/728/CH6/EX6.9/Ex6_9_ans.txt b/728/CH6/EX6.9/Ex6_9_ans.txt
new file mode 100755
index 000000000..d38ffdb28
--- /dev/null
+++ b/728/CH6/EX6.9/Ex6_9_ans.txt
@@ -0,0 +1,7 @@
+
+ Scattering matrix =
+
+ 0. 0.1 0.9440609
+ 0.9440609 0. 0.1
+ 0.1 0.9440609 0.
+ \ No newline at end of file
diff --git a/728/CH7/EX7.1/Ex7_1.sce b/728/CH7/EX7.1/Ex7_1.sce
new file mode 100755
index 000000000..e6da56039
--- /dev/null
+++ b/728/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate the SWR of the transmission line
+//Exa:7.1
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=4;//in cm
+b=2.5;//in cm
+f=10*10^9;//in Hz
+d=0.1;//distance between 2 minimum power points(in cm)
+//For TE10 mode:
+wl_c=2*a;
+wl_o=c/f;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+S=wl_g/(%pi*d);
+disp(S,'Voltage standing wave ratio ='); \ No newline at end of file
diff --git a/728/CH7/EX7.1/Ex7_1.txt b/728/CH7/EX7.1/Ex7_1.txt
new file mode 100755
index 000000000..e6da56039
--- /dev/null
+++ b/728/CH7/EX7.1/Ex7_1.txt
@@ -0,0 +1,17 @@
+//Caption:Calculate the SWR of the transmission line
+//Exa:7.1
+clc;
+clear;
+close;
+//Given:
+c=3*10^10;//in cm/s
+a=4;//in cm
+b=2.5;//in cm
+f=10*10^9;//in Hz
+d=0.1;//distance between 2 minimum power points(in cm)
+//For TE10 mode:
+wl_c=2*a;
+wl_o=c/f;
+wl_g=wl_o/sqrt(1-(wl_o/wl_c)^2);
+S=wl_g/(%pi*d);
+disp(S,'Voltage standing wave ratio ='); \ No newline at end of file
diff --git a/728/CH7/EX7.1/Ex7_1_ans.txt b/728/CH7/EX7.1/Ex7_1_ans.txt
new file mode 100755
index 000000000..14957e63e
--- /dev/null
+++ b/728/CH7/EX7.1/Ex7_1_ans.txt
@@ -0,0 +1,5 @@
+
+ Voltage standing wave ratio =
+
+ 10.301015
+ \ No newline at end of file
diff --git a/728/CH7/EX7.2/Ex7_2.sce b/728/CH7/EX7.2/Ex7_2.sce
new file mode 100755
index 000000000..514d3b9e6
--- /dev/null
+++ b/728/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,11 @@
+//Caption:Calculate the SWR of the main waveguide
+//Exa:7.2
+clc;
+clear;
+close;
+//Given:
+P_i=300;//in mW
+P_r=10;//in mW
+p=sqrt(P_r/P_i);
+S=(1+p)/(1-p);
+disp(S,'Voltage standing wave ratio =');
diff --git a/728/CH7/EX7.2/Ex7_2.txt b/728/CH7/EX7.2/Ex7_2.txt
new file mode 100755
index 000000000..514d3b9e6
--- /dev/null
+++ b/728/CH7/EX7.2/Ex7_2.txt
@@ -0,0 +1,11 @@
+//Caption:Calculate the SWR of the main waveguide
+//Exa:7.2
+clc;
+clear;
+close;
+//Given:
+P_i=300;//in mW
+P_r=10;//in mW
+p=sqrt(P_r/P_i);
+S=(1+p)/(1-p);
+disp(S,'Voltage standing wave ratio =');
diff --git a/728/CH7/EX7.2/Ex7_2_ans.txt b/728/CH7/EX7.2/Ex7_2_ans.txt
new file mode 100755
index 000000000..f422ebfef
--- /dev/null
+++ b/728/CH7/EX7.2/Ex7_2_ans.txt
@@ -0,0 +1,5 @@
+
+ Voltage standing wave ratio =
+
+ 1.4467052
+ \ No newline at end of file
diff --git a/728/CH7/EX7.3/Ex7_3.sce b/728/CH7/EX7.3/Ex7_3.sce
new file mode 100755
index 000000000..fd065113f
--- /dev/null
+++ b/728/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,11 @@
+//Caption:Calculate the SWR of the waveguide
+//Exa:7.3
+clc;
+clear;
+close;
+//Given:
+P_i=2.5;//in mW
+P_r=0.15;//in mW
+p=sqrt(P_r/P_i);
+S=(1+p)/(1-p);
+disp(S,'Voltage standing wave ratio =');
diff --git a/728/CH7/EX7.3/Ex7_3.txt b/728/CH7/EX7.3/Ex7_3.txt
new file mode 100755
index 000000000..fd065113f
--- /dev/null
+++ b/728/CH7/EX7.3/Ex7_3.txt
@@ -0,0 +1,11 @@
+//Caption:Calculate the SWR of the waveguide
+//Exa:7.3
+clc;
+clear;
+close;
+//Given:
+P_i=2.5;//in mW
+P_r=0.15;//in mW
+p=sqrt(P_r/P_i);
+S=(1+p)/(1-p);
+disp(S,'Voltage standing wave ratio =');
diff --git a/728/CH7/EX7.3/Ex7_3_ans.txt b/728/CH7/EX7.3/Ex7_3_ans.txt
new file mode 100755
index 000000000..22ac0931c
--- /dev/null
+++ b/728/CH7/EX7.3/Ex7_3_ans.txt
@@ -0,0 +1,4 @@
+Voltage standing wave ratio =
+
+ 1.6488276
+ \ No newline at end of file
diff --git a/728/CH7/EX7.4/Ex7_4.sce b/728/CH7/EX7.4/Ex7_4.sce
new file mode 100755
index 000000000..cf0c5872d
--- /dev/null
+++ b/728/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,13 @@
+//Caption:Calculate the value of reflected power
+//Exa:7.4
+clc;
+clear;
+close;
+//Given:
+P_i=4.5;//in mW
+S=2;//VSWR
+C=30;//in dB
+p=(S-1)/(S+1);
+P_f=P_i/(10^(C/10));
+P_r=p^2*P_i;
+disp(P_r,'Reflected power (in watts) ='); \ No newline at end of file
diff --git a/728/CH7/EX7.4/Ex7_4.txt b/728/CH7/EX7.4/Ex7_4.txt
new file mode 100755
index 000000000..cf0c5872d
--- /dev/null
+++ b/728/CH7/EX7.4/Ex7_4.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate the value of reflected power
+//Exa:7.4
+clc;
+clear;
+close;
+//Given:
+P_i=4.5;//in mW
+S=2;//VSWR
+C=30;//in dB
+p=(S-1)/(S+1);
+P_f=P_i/(10^(C/10));
+P_r=p^2*P_i;
+disp(P_r,'Reflected power (in watts) ='); \ No newline at end of file
diff --git a/728/CH7/EX7.4/Ex7_4_ans.txt b/728/CH7/EX7.4/Ex7_4_ans.txt
new file mode 100755
index 000000000..f0d2dc880
--- /dev/null
+++ b/728/CH7/EX7.4/Ex7_4_ans.txt
@@ -0,0 +1,5 @@
+
+ Reflected power (in watts) =
+
+ 0.5
+ \ No newline at end of file
diff --git a/728/CH8/EX8.1/Ex8_1.sce b/728/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..2695ffd59
--- /dev/null
+++ b/728/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,26 @@
+//Caption:Calculate (i)-dc electron velocity,(ii)-dc phase constant, (iii)-plasma frequency, (iv)-reduced plasma frequency for R=0.4, (v)-dc beam current beam density, (vi)-instantaneous beam current density
+//Exa:8.1
+clc;
+clear;
+close;
+V_o=14.5*10^3;//in volts
+I_o=1.4;//in A
+f=10*10^9;//in Hz
+p_o=10^-6;//in c/m^3
+p=10^-8;//in c/m^3
+v=10^5;//in m/s
+R=0.4;
+v_o=0.593*10^6*sqrt(V_o);
+k=2*%pi*f/v_o;
+w_p=[1.759*10^11*(10^-6/(8.854*10^-12))]^0.5;
+w_q=R*w_p;
+J_o=p_o*v_o;
+J=p*v_o+p_o*v;
+disp(v_o,'Dc electron velocity (in m/s) =');
+disp(k,'Dc phase constant (in rad/s) =');
+disp(w_p,'Plasma frequency (in rad/s) =');
+disp(w_q,'Reduced plasma frequency (in rad/s) =');
+disp(J_o,'Dc beam current density (in A/sq. m) =');
+disp(J,'Instantaneous beam current density(in A/sq. m) =');
+
+//Answer in book are wrongly written as: (Dc phase constant =1.41* 10^8 rad/sec) \ No newline at end of file
diff --git a/728/CH8/EX8.1/Ex8_1.txt b/728/CH8/EX8.1/Ex8_1.txt
new file mode 100755
index 000000000..2695ffd59
--- /dev/null
+++ b/728/CH8/EX8.1/Ex8_1.txt
@@ -0,0 +1,26 @@
+//Caption:Calculate (i)-dc electron velocity,(ii)-dc phase constant, (iii)-plasma frequency, (iv)-reduced plasma frequency for R=0.4, (v)-dc beam current beam density, (vi)-instantaneous beam current density
+//Exa:8.1
+clc;
+clear;
+close;
+V_o=14.5*10^3;//in volts
+I_o=1.4;//in A
+f=10*10^9;//in Hz
+p_o=10^-6;//in c/m^3
+p=10^-8;//in c/m^3
+v=10^5;//in m/s
+R=0.4;
+v_o=0.593*10^6*sqrt(V_o);
+k=2*%pi*f/v_o;
+w_p=[1.759*10^11*(10^-6/(8.854*10^-12))]^0.5;
+w_q=R*w_p;
+J_o=p_o*v_o;
+J=p*v_o+p_o*v;
+disp(v_o,'Dc electron velocity (in m/s) =');
+disp(k,'Dc phase constant (in rad/s) =');
+disp(w_p,'Plasma frequency (in rad/s) =');
+disp(w_q,'Reduced plasma frequency (in rad/s) =');
+disp(J_o,'Dc beam current density (in A/sq. m) =');
+disp(J,'Instantaneous beam current density(in A/sq. m) =');
+
+//Answer in book are wrongly written as: (Dc phase constant =1.41* 10^8 rad/sec) \ No newline at end of file
diff --git a/728/CH8/EX8.1/Ex8_1_ans.txt b/728/CH8/EX8.1/Ex8_1_ans.txt
new file mode 100755
index 000000000..7bc8166be
--- /dev/null
+++ b/728/CH8/EX8.1/Ex8_1_ans.txt
@@ -0,0 +1,24 @@
+ Dc electron velocity (in m/s) =
+
+ 71406656.
+
+ Dc phase constant (in rad/s) =
+
+ 879.91592
+
+ Plasma frequency (in rad/s) =
+
+ 1.409D+08
+
+ Reduced plasma frequency (in rad/s) =
+
+ 56379751.
+
+ Dc beam current density (in A/sq. m) =
+
+ 71.406656
+
+ Instantaneous beam current density(in A/sq. m) =
+
+ 0.8140666
+ \ No newline at end of file
diff --git a/728/CH8/EX8.10/Ex8_10.sce b/728/CH8/EX8.10/Ex8_10.sce
new file mode 100755
index 000000000..8539a0c54
--- /dev/null
+++ b/728/CH8/EX8.10/Ex8_10.sce
@@ -0,0 +1,23 @@
+//Caption:Calculate (i)-electron velocity,(ii)-dc transit time, (iii)-input voltage for maximum output voltage,(iv)-voltage gain in dB
+//Exa:8.10
+clc;
+clear;
+close;
+V_o=900;//in volts
+I_o=30*10^-3;//in A
+f=8*10^9;//in Hz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohm
+v_o=0.593*10^6*sqrt(V_o);
+T_o=l/v_o;
+Theeta_o=(2*%pi*f)*T_o;//Transit angles between cavities(in radian)
+Theeta_g=(2*%pi*f)*d/v_o;//Average gap transit angle (in radian)
+b=sin(Theeta_g/2)/(Theeta_g/2);
+V_in_max=V_o*3.68/(b*Theeta_o);
+//As, {J(X)/X=0.582}
+A_r=b^2*Theeta_o*0.582*R_sh/(30*10^3*1.841);
+disp(v_o,'Electron velocity (in m/s) =');
+disp(T_o,'Dc Transit Time (in sec)=');
+disp(V_in_max,'Maximum input voltage (in volts) =');
+disp(A_r,'Voltage gain (in dB) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.10/Ex8_10.txt b/728/CH8/EX8.10/Ex8_10.txt
new file mode 100755
index 000000000..8539a0c54
--- /dev/null
+++ b/728/CH8/EX8.10/Ex8_10.txt
@@ -0,0 +1,23 @@
+//Caption:Calculate (i)-electron velocity,(ii)-dc transit time, (iii)-input voltage for maximum output voltage,(iv)-voltage gain in dB
+//Exa:8.10
+clc;
+clear;
+close;
+V_o=900;//in volts
+I_o=30*10^-3;//in A
+f=8*10^9;//in Hz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohm
+v_o=0.593*10^6*sqrt(V_o);
+T_o=l/v_o;
+Theeta_o=(2*%pi*f)*T_o;//Transit angles between cavities(in radian)
+Theeta_g=(2*%pi*f)*d/v_o;//Average gap transit angle (in radian)
+b=sin(Theeta_g/2)/(Theeta_g/2);
+V_in_max=V_o*3.68/(b*Theeta_o);
+//As, {J(X)/X=0.582}
+A_r=b^2*Theeta_o*0.582*R_sh/(30*10^3*1.841);
+disp(v_o,'Electron velocity (in m/s) =');
+disp(T_o,'Dc Transit Time (in sec)=');
+disp(V_in_max,'Maximum input voltage (in volts) =');
+disp(A_r,'Voltage gain (in dB) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.10/Ex8_10_ans.txt b/728/CH8/EX8.10/Ex8_10_ans.txt
new file mode 100755
index 000000000..ce294587d
--- /dev/null
+++ b/728/CH8/EX8.10/Ex8_10_ans.txt
@@ -0,0 +1,16 @@
+Electron velocity (in m/s) =
+
+ 17790000.
+
+ Dc Transit Time (in sec)=
+
+ 2.248D-09
+
+ Maximum input voltage (in volts) =
+
+ 41.922525
+
+ Voltage gain (in dB) =
+
+ 23.277687
+ \ No newline at end of file
diff --git a/728/CH8/EX8.11/Ex8_11.sce b/728/CH8/EX8.11/Ex8_11.sce
new file mode 100755
index 000000000..8d309ed19
--- /dev/null
+++ b/728/CH8/EX8.11/Ex8_11.sce
@@ -0,0 +1,24 @@
+//Caption:Calculate (i)-dc electron velocity,(ii)-dc phase constant, (iii)-plasma frequency, (iv)-reduced plasma frequency for R=0.5, (v)-dc beam current beam density, (vi)-instantaneous beam current density
+//Exa:8.11
+clc;
+clear;
+close;
+V_o=20*10^3;//in volts
+I_o=2;//in A
+f=10*10^9;//in Hz
+p_o=10^-6;//in c/m^3
+p=10^-8;//in c/m^3
+v=10^5;//in m/s
+R=0.5;
+v_o=0.593*10^6*sqrt(V_o);
+k=2*%pi*f/v_o;
+w_p=[1.759*10^11*(10^-6/(8.854*10^-12))]^0.5;
+w_q=R*w_p;
+J_o=p_o*v_o;
+J=p*v_o-p_o*v;
+disp(v_o,'Dc electron velocity (in m/s) =');
+disp(k,'Dc phase constant (in rad/s) =');
+disp(w_p,'Plasma frequency (in rad/s) =');
+disp(w_q,'Reduced plasma frequency (in rad/s) =');
+disp(J_o,'Dc beam current density (in A/sq. m) =');
+disp(J,'Instantaneous beam current density(in A/sq. m) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.11/Ex8_11.txt b/728/CH8/EX8.11/Ex8_11.txt
new file mode 100755
index 000000000..8d309ed19
--- /dev/null
+++ b/728/CH8/EX8.11/Ex8_11.txt
@@ -0,0 +1,24 @@
+//Caption:Calculate (i)-dc electron velocity,(ii)-dc phase constant, (iii)-plasma frequency, (iv)-reduced plasma frequency for R=0.5, (v)-dc beam current beam density, (vi)-instantaneous beam current density
+//Exa:8.11
+clc;
+clear;
+close;
+V_o=20*10^3;//in volts
+I_o=2;//in A
+f=10*10^9;//in Hz
+p_o=10^-6;//in c/m^3
+p=10^-8;//in c/m^3
+v=10^5;//in m/s
+R=0.5;
+v_o=0.593*10^6*sqrt(V_o);
+k=2*%pi*f/v_o;
+w_p=[1.759*10^11*(10^-6/(8.854*10^-12))]^0.5;
+w_q=R*w_p;
+J_o=p_o*v_o;
+J=p*v_o-p_o*v;
+disp(v_o,'Dc electron velocity (in m/s) =');
+disp(k,'Dc phase constant (in rad/s) =');
+disp(w_p,'Plasma frequency (in rad/s) =');
+disp(w_q,'Reduced plasma frequency (in rad/s) =');
+disp(J_o,'Dc beam current density (in A/sq. m) =');
+disp(J,'Instantaneous beam current density(in A/sq. m) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.11/Ex8_11_ans.txt b/728/CH8/EX8.11/Ex8_11_ans.txt
new file mode 100755
index 000000000..95c0fefc5
--- /dev/null
+++ b/728/CH8/EX8.11/Ex8_11_ans.txt
@@ -0,0 +1,24 @@
+Dc electron velocity (in m/s) =
+
+ 83862864.
+
+ Dc phase constant (in rad/s) =
+
+ 749.22141
+
+ Plasma frequency (in rad/s) =
+
+ 1.409D+08
+
+ Reduced plasma frequency (in rad/s) =
+
+ 70474689.
+
+ Dc beam current density (in A/sq. m) =
+
+ 83.862864
+
+ Instantaneous beam current density(in A/sq. m) =
+
+ 0.7386286
+ \ No newline at end of file
diff --git a/728/CH8/EX8.12/Ex8_12.sce b/728/CH8/EX8.12/Ex8_12.sce
new file mode 100755
index 000000000..c1c2efa31
--- /dev/null
+++ b/728/CH8/EX8.12/Ex8_12.sce
@@ -0,0 +1,12 @@
+//Caption:Calculate the gap transit angle
+//Exa:8.12
+clc;
+clear;
+close;
+V_o=1000;//Anode voltage(in volts)
+gap=0.002;//in m
+f=5*10^9;//in Hz
+L=2.463*10^-3;//length of drift region (in m)
+u_o=5.93*10^5*sqrt(V_o);//in m/s
+Theeta_g=2*%pi*f*2*10^-3/u_o;//radians
+disp(Theeta_g,'Transit angle(in radians) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.12/Ex8_12.txt b/728/CH8/EX8.12/Ex8_12.txt
new file mode 100755
index 000000000..c1c2efa31
--- /dev/null
+++ b/728/CH8/EX8.12/Ex8_12.txt
@@ -0,0 +1,12 @@
+//Caption:Calculate the gap transit angle
+//Exa:8.12
+clc;
+clear;
+close;
+V_o=1000;//Anode voltage(in volts)
+gap=0.002;//in m
+f=5*10^9;//in Hz
+L=2.463*10^-3;//length of drift region (in m)
+u_o=5.93*10^5*sqrt(V_o);//in m/s
+Theeta_g=2*%pi*f*2*10^-3/u_o;//radians
+disp(Theeta_g,'Transit angle(in radians) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.12/Ex8_12_ans.txt b/728/CH8/EX8.12/Ex8_12_ans.txt
new file mode 100755
index 000000000..6a56e6958
--- /dev/null
+++ b/728/CH8/EX8.12/Ex8_12_ans.txt
@@ -0,0 +1,4 @@
+Transit angle(in radians) =
+
+ 3.35062
+ \ No newline at end of file
diff --git a/728/CH8/EX8.13/Ex8_13.sce b/728/CH8/EX8.13/Ex8_13.sce
new file mode 100755
index 000000000..b1098d54b
--- /dev/null
+++ b/728/CH8/EX8.13/Ex8_13.sce
@@ -0,0 +1,29 @@
+//Caption:Calculate (i)-i/p rf voltage,(ii)-voltage gain, (iii)-efficiency
+//Exa:8.13
+clc;
+clear;
+close;
+V_o=1200;//in volts
+I_o=30*10^-3;//in A
+f=10*10^9;//inHz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohms
+v_o=0.593*10^6*sqrt(V_o);
+Theeta_o=2*%pi*f*l/(20.54*10^6);
+X=1.84;//for maximum output power
+V_max=2*X*V_o/122.347;
+Theeta_g=122.347*10^-3/(4*10^-2);
+Beeta_i=sin(Theeta_g/2)/(Theeta_g/2);
+V_1_max=V_max/Beeta_i;
+J=0.58;
+Beeta_o=Beeta_i;
+I_2=2*I_o*J;
+V_2=Beeta_o*I_2*R_sh;
+A_v=V_2/V_1_max;//in dB
+eff=0.58*(V_2/V_o)*100;//in percentage
+disp(V_1_max,'Input rf voltage(in volts) =');
+disp(A_v,'Voltage gain (in dB) =');
+disp(eff,'Maximum efficiency (in percentage) =');
+
+//Answer in book is wrongly given as: A_v=24.33 dB \ No newline at end of file
diff --git a/728/CH8/EX8.13/Ex8_13.txt b/728/CH8/EX8.13/Ex8_13.txt
new file mode 100755
index 000000000..b1098d54b
--- /dev/null
+++ b/728/CH8/EX8.13/Ex8_13.txt
@@ -0,0 +1,29 @@
+//Caption:Calculate (i)-i/p rf voltage,(ii)-voltage gain, (iii)-efficiency
+//Exa:8.13
+clc;
+clear;
+close;
+V_o=1200;//in volts
+I_o=30*10^-3;//in A
+f=10*10^9;//inHz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohms
+v_o=0.593*10^6*sqrt(V_o);
+Theeta_o=2*%pi*f*l/(20.54*10^6);
+X=1.84;//for maximum output power
+V_max=2*X*V_o/122.347;
+Theeta_g=122.347*10^-3/(4*10^-2);
+Beeta_i=sin(Theeta_g/2)/(Theeta_g/2);
+V_1_max=V_max/Beeta_i;
+J=0.58;
+Beeta_o=Beeta_i;
+I_2=2*I_o*J;
+V_2=Beeta_o*I_2*R_sh;
+A_v=V_2/V_1_max;//in dB
+eff=0.58*(V_2/V_o)*100;//in percentage
+disp(V_1_max,'Input rf voltage(in volts) =');
+disp(A_v,'Voltage gain (in dB) =');
+disp(eff,'Maximum efficiency (in percentage) =');
+
+//Answer in book is wrongly given as: A_v=24.33 dB \ No newline at end of file
diff --git a/728/CH8/EX8.13/Ex8_13_ans.txt b/728/CH8/EX8.13/Ex8_13_ans.txt
new file mode 100755
index 000000000..822057a54
--- /dev/null
+++ b/728/CH8/EX8.13/Ex8_13_ans.txt
@@ -0,0 +1,12 @@
+
+ Input rf voltage(in volts) =
+
+ 55.247474
+
+ Voltage gain (in dB) =
+
+ 16.460769
+
+ Maximum efficiency (in percentage) =
+
+ 43.955103 \ No newline at end of file
diff --git a/728/CH8/EX8.14/Ex8_14.sce b/728/CH8/EX8.14/Ex8_14.sce
new file mode 100755
index 000000000..f9b921363
--- /dev/null
+++ b/728/CH8/EX8.14/Ex8_14.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate (i)-cyclotron angular frequency, (ii)-cut-off voltage, (iii)-cut-off magnetic flux
+//Exa:8.14
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+a=0.04;
+b=0.08;
+V_o=30*10^3;//in volts
+I_o=80;//in A
+B_o=0.01;//in weber/sq.m
+w=(e_m_ratio)*B_o;
+disp(w,'Cyclotron angular frequency( in rad/s) =');
+V_c={(e_m_ratio)*B_o^2*b^2*[1-(a/b)^2]^2}/8;
+disp(V_c,'Cut-off voltage (in volts) =');
+B_c=sqrt(8*V_o/e_m_ratio)/[[1-(a/b)^2]*(b)];//in weber/m^2
+disp(B_c*10^3,'Cut-off magnetic flux density (in milli weber/sq. m) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.14/Ex8_14.txt b/728/CH8/EX8.14/Ex8_14.txt
new file mode 100755
index 000000000..f9b921363
--- /dev/null
+++ b/728/CH8/EX8.14/Ex8_14.txt
@@ -0,0 +1,17 @@
+//Caption:Calculate (i)-cyclotron angular frequency, (ii)-cut-off voltage, (iii)-cut-off magnetic flux
+//Exa:8.14
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+a=0.04;
+b=0.08;
+V_o=30*10^3;//in volts
+I_o=80;//in A
+B_o=0.01;//in weber/sq.m
+w=(e_m_ratio)*B_o;
+disp(w,'Cyclotron angular frequency( in rad/s) =');
+V_c={(e_m_ratio)*B_o^2*b^2*[1-(a/b)^2]^2}/8;
+disp(V_c,'Cut-off voltage (in volts) =');
+B_c=sqrt(8*V_o/e_m_ratio)/[[1-(a/b)^2]*(b)];//in weber/m^2
+disp(B_c*10^3,'Cut-off magnetic flux density (in milli weber/sq. m) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.14/Ex8_14_ans.txt b/728/CH8/EX8.14/Ex8_14_ans.txt
new file mode 100755
index 000000000..035946828
--- /dev/null
+++ b/728/CH8/EX8.14/Ex8_14_ans.txt
@@ -0,0 +1,12 @@
+Cyclotron angular frequency( in rad/s) =
+
+ 1.759D+09
+
+ Cut-off voltage (in volts) =
+
+ 7915.5
+
+ Cut-off magnetic flux density (in milli weber/sq. m) =
+
+ 19.468005
+ \ No newline at end of file
diff --git a/728/CH8/EX8.15/Ex8_15.sce b/728/CH8/EX8.15/Ex8_15.sce
new file mode 100755
index 000000000..2a5f1cb2b
--- /dev/null
+++ b/728/CH8/EX8.15/Ex8_15.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input power ,(ii)-output power in watts, (iii)-efficiency
+//Exa:8.15
+clc;
+clear;
+close;
+n=2;
+V_o=280;//in volts
+I_o=22*10^-3;//in A
+V_i=30;//in volts
+J=1.25;//J(X')
+P_dc=V_o*I_o;
+P_ac=2*V_o*I_o*J/(2*n*%pi-%pi/2);
+eff=(P_ac/P_dc)*100;
+disp(P_dc,'Input power (in watts) =');
+disp(P_ac,'Output power (in watts) =');
+disp(eff,'Efficiency (in percent) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.15/Ex8_15.txt b/728/CH8/EX8.15/Ex8_15.txt
new file mode 100755
index 000000000..2a5f1cb2b
--- /dev/null
+++ b/728/CH8/EX8.15/Ex8_15.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input power ,(ii)-output power in watts, (iii)-efficiency
+//Exa:8.15
+clc;
+clear;
+close;
+n=2;
+V_o=280;//in volts
+I_o=22*10^-3;//in A
+V_i=30;//in volts
+J=1.25;//J(X')
+P_dc=V_o*I_o;
+P_ac=2*V_o*I_o*J/(2*n*%pi-%pi/2);
+eff=(P_ac/P_dc)*100;
+disp(P_dc,'Input power (in watts) =');
+disp(P_ac,'Output power (in watts) =');
+disp(eff,'Efficiency (in percent) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.15/Ex8_15_ans.txt b/728/CH8/EX8.15/Ex8_15_ans.txt
new file mode 100755
index 000000000..ac3ef33f1
--- /dev/null
+++ b/728/CH8/EX8.15/Ex8_15_ans.txt
@@ -0,0 +1,12 @@
+
+ Input power (in watts) =
+
+ 6.16
+
+ Output power (in watts) =
+
+ 1.4005635
+
+ Efficiency (in percent) =
+
+ 22.73642 \ No newline at end of file
diff --git a/728/CH8/EX8.16/Ex8_16.sce b/728/CH8/EX8.16/Ex8_16.sce
new file mode 100755
index 000000000..37d76d780
--- /dev/null
+++ b/728/CH8/EX8.16/Ex8_16.sce
@@ -0,0 +1,20 @@
+//Caption:Calculate (i)-repeller voltage V_r ,(ii)-beam current necessary to give gap voltage of 200V
+//Exa:8.16
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+V_o=300;//in volts
+R_sh=20*10^3;//in ohms
+f=8*10^9;//inHz
+w=2*%pi*f;
+n=2;//mode
+L=0.001;//spacing between repeller & cavity (in m)
+x=(e_m_ratio)*(2*%pi*n-%pi/2)^2/(8*w^2*L^2);
+volt_diff=sqrt(V_o/(x));
+V_r=(volt_diff)+V_o;//repeller volatge
+J=0.582;
+V_1=200;//given (in volts)
+I_o=V_1/(R_sh*2*J);
+disp(V_r,'Repeller voltage(in volts) =');
+disp(I_o*10^3,'Necessary beam current (in milliAmp.s) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.16/Ex8_16.txt b/728/CH8/EX8.16/Ex8_16.txt
new file mode 100755
index 000000000..37d76d780
--- /dev/null
+++ b/728/CH8/EX8.16/Ex8_16.txt
@@ -0,0 +1,20 @@
+//Caption:Calculate (i)-repeller voltage V_r ,(ii)-beam current necessary to give gap voltage of 200V
+//Exa:8.16
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+V_o=300;//in volts
+R_sh=20*10^3;//in ohms
+f=8*10^9;//inHz
+w=2*%pi*f;
+n=2;//mode
+L=0.001;//spacing between repeller & cavity (in m)
+x=(e_m_ratio)*(2*%pi*n-%pi/2)^2/(8*w^2*L^2);
+volt_diff=sqrt(V_o/(x));
+V_r=(volt_diff)+V_o;//repeller volatge
+J=0.582;
+V_1=200;//given (in volts)
+I_o=V_1/(R_sh*2*J);
+disp(V_r,'Repeller voltage(in volts) =');
+disp(I_o*10^3,'Necessary beam current (in milliAmp.s) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.16/Ex8_16_ans.txt b/728/CH8/EX8.16/Ex8_16_ans.txt
new file mode 100755
index 000000000..418e75500
--- /dev/null
+++ b/728/CH8/EX8.16/Ex8_16_ans.txt
@@ -0,0 +1,8 @@
+Repeller voltage(in volts) =
+
+ 833.97957
+
+ Necessary beam current (in milliAmp.s) =
+
+ 8.5910653
+ \ No newline at end of file
diff --git a/728/CH8/EX8.2/Ex8_2.sce b/728/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..39e1bdc93
--- /dev/null
+++ b/728/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input rms voltage,(ii)-output rms voltage, (iii)-power delivered to load
+//Exa:8.2
+clc;
+clear;
+close;
+A_v=15;//in dB
+P_i=5*10^-3;//in W
+R_sh_i=30000;//in ohms
+R_sh_o=40000;//in ohms
+R_l=20000;//in ohms
+V_i=sqrt(P_i*R_sh_i);
+V_o=10^((A_v/20))*12.25;
+P_out=V_o^2/R_l;
+disp(V_i,'Input rms voltage (in volts) =');
+disp(V_o,'Output rms voltage (in volts) =');
+disp(P_out,'Power delivered to load (in watts) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.2/Ex8_2.txt b/728/CH8/EX8.2/Ex8_2.txt
new file mode 100755
index 000000000..39e1bdc93
--- /dev/null
+++ b/728/CH8/EX8.2/Ex8_2.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input rms voltage,(ii)-output rms voltage, (iii)-power delivered to load
+//Exa:8.2
+clc;
+clear;
+close;
+A_v=15;//in dB
+P_i=5*10^-3;//in W
+R_sh_i=30000;//in ohms
+R_sh_o=40000;//in ohms
+R_l=20000;//in ohms
+V_i=sqrt(P_i*R_sh_i);
+V_o=10^((A_v/20))*12.25;
+P_out=V_o^2/R_l;
+disp(V_i,'Input rms voltage (in volts) =');
+disp(V_o,'Output rms voltage (in volts) =');
+disp(P_out,'Power delivered to load (in watts) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.2/Ex8_2_ans.txt b/728/CH8/EX8.2/Ex8_2_ans.txt
new file mode 100755
index 000000000..daf76b053
--- /dev/null
+++ b/728/CH8/EX8.2/Ex8_2_ans.txt
@@ -0,0 +1,12 @@
+Input rms voltage (in volts) =
+
+ 12.247449
+
+ Output rms voltage (in volts) =
+
+ 68.886812
+
+ Power delivered to load (in watts) =
+
+ 0.2372696
+ \ No newline at end of file
diff --git a/728/CH8/EX8.3/Ex8_3.sce b/728/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..bf7dd832e
--- /dev/null
+++ b/728/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input power in watts,(ii)-output power in watts, (iii)-efficiency
+//Exa:8.3
+clc;
+clear;
+close;
+n=2;
+V_o=300;//in volts
+I_o=20*10^-3;//in A
+V_i=40;//in volts
+J=1.25;//J(X')
+P_dc=V_o*I_o;
+P_ac=2*V_o*I_o*J/(2*n*%pi-%pi/2);
+eff=(P_ac/P_dc)*100;
+disp(P_dc,'Input power (in watts) =');
+disp(P_ac,'Output power (in watts) =');
+disp(eff,'Efficiency (in percent) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.3/Ex8_3.txt b/728/CH8/EX8.3/Ex8_3.txt
new file mode 100755
index 000000000..bf7dd832e
--- /dev/null
+++ b/728/CH8/EX8.3/Ex8_3.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-input power in watts,(ii)-output power in watts, (iii)-efficiency
+//Exa:8.3
+clc;
+clear;
+close;
+n=2;
+V_o=300;//in volts
+I_o=20*10^-3;//in A
+V_i=40;//in volts
+J=1.25;//J(X')
+P_dc=V_o*I_o;
+P_ac=2*V_o*I_o*J/(2*n*%pi-%pi/2);
+eff=(P_ac/P_dc)*100;
+disp(P_dc,'Input power (in watts) =');
+disp(P_ac,'Output power (in watts) =');
+disp(eff,'Efficiency (in percent) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.3/Ex8_3_ans.txt b/728/CH8/EX8.3/Ex8_3_ans.txt
new file mode 100755
index 000000000..5720a3d55
--- /dev/null
+++ b/728/CH8/EX8.3/Ex8_3_ans.txt
@@ -0,0 +1,11 @@
+Input power (in watts) =
+
+ 6.
+
+ Output power (in watts) =
+
+ 1.3641852
+
+ Efficiency (in percent) =
+
+ 22.73642 \ No newline at end of file
diff --git a/728/CH8/EX8.4/Ex8_4.sce b/728/CH8/EX8.4/Ex8_4.sce
new file mode 100755
index 000000000..e1a3590db
--- /dev/null
+++ b/728/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,23 @@
+//Caption:Calculate (i)-electron velocity,(ii)-dc transit time, (iii)-input voltage for maximum output voltage,(iv)-voltage gain in dB
+//Exa:8.4
+clc;
+clear;
+close;
+V_o=900;//in volts
+I_o=30*10^-3;//in A
+f=8*10^9;//in Hz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohm
+v_o=0.593*10^6*sqrt(V_o);
+T_o=l/v_o;
+Theeta_o=(2*%pi*f)*T_o;//Transit angles between cavities(in radian)
+Theeta_g=(2*%pi*f)*d/v_o;//Average gap transit angle (in radian)
+b=sin(Theeta_g/2)/(Theeta_g/2);
+V_in_max=V_o*3.68/(b*Theeta_o);
+//As, {J(X)/X=0.582}
+A_r=b^2*Theeta_o*0.582*R_sh/(30*10^3*1.841);
+disp(v_o,'Electron velocity (in m/s) =');
+disp(T_o,'Dc Transit Time (in sec)=');
+disp(V_in_max,'Maximum input voltage (in volts) =');
+disp(A_r,'Voltage gain (in dB) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.4/Ex8_4.txt b/728/CH8/EX8.4/Ex8_4.txt
new file mode 100755
index 000000000..e1a3590db
--- /dev/null
+++ b/728/CH8/EX8.4/Ex8_4.txt
@@ -0,0 +1,23 @@
+//Caption:Calculate (i)-electron velocity,(ii)-dc transit time, (iii)-input voltage for maximum output voltage,(iv)-voltage gain in dB
+//Exa:8.4
+clc;
+clear;
+close;
+V_o=900;//in volts
+I_o=30*10^-3;//in A
+f=8*10^9;//in Hz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohm
+v_o=0.593*10^6*sqrt(V_o);
+T_o=l/v_o;
+Theeta_o=(2*%pi*f)*T_o;//Transit angles between cavities(in radian)
+Theeta_g=(2*%pi*f)*d/v_o;//Average gap transit angle (in radian)
+b=sin(Theeta_g/2)/(Theeta_g/2);
+V_in_max=V_o*3.68/(b*Theeta_o);
+//As, {J(X)/X=0.582}
+A_r=b^2*Theeta_o*0.582*R_sh/(30*10^3*1.841);
+disp(v_o,'Electron velocity (in m/s) =');
+disp(T_o,'Dc Transit Time (in sec)=');
+disp(V_in_max,'Maximum input voltage (in volts) =');
+disp(A_r,'Voltage gain (in dB) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.4/Ex8_4_ans.txt b/728/CH8/EX8.4/Ex8_4_ans.txt
new file mode 100755
index 000000000..99b7297b9
--- /dev/null
+++ b/728/CH8/EX8.4/Ex8_4_ans.txt
@@ -0,0 +1,16 @@
+ Electron velocity (in m/s) =
+
+ 17790000.
+
+ Dc Transit Time (in sec)=
+
+ 2.248D-09
+
+ Maximum input voltage (in volts) =
+
+ 41.922525
+
+ Voltage gain (in dB) =
+
+ 23.277687
+ \ No newline at end of file
diff --git a/728/CH8/EX8.5/Ex8_5.sce b/728/CH8/EX8.5/Ex8_5.sce
new file mode 100755
index 000000000..b7fa37ac5
--- /dev/null
+++ b/728/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,28 @@
+//Caption:Calculate (i)-i/p microwave voltage,(ii)-voltage gain, (iii)-efficiency of amplifier,(iv)-beam loading conductance
+//Exa:8.5
+clc;
+clear;
+close;
+V_o=1200;//in volts
+I_o=28*10^-3;//in A
+f=8*10^9;//inHz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohms
+V_p_max=1200*3.68*0.593*10^6*sqrt(V_o)/(2*%pi*f*l);
+Theeta_g=(2*%pi*f)*d/(0.593*10^6*sqrt(V_o));//transit angle (in rad)
+beeta=sin(Theeta_g/2)/(Theeta_g/2);
+V_i_max=V_p_max/beeta;
+Beeta_o=0.768;
+J=0.582;//J(X)
+A_v=(Beeta_o)^2*97.88*J*R_sh/(1200/(28*10^-3*1.841));//calculating voltage gain
+eff=[0.58*[2*28*10^-3*J*Beeta_o*R_sh]/V_o]*100;//calculating efficiency
+G_o=23.3*10^-6;
+G_b=(G_o/2)*{Beeta_o^2-Beeta_o*cos(Theeta_g)};//beam loading conductance
+R_b=1/(G_b*1000);//beam loading resistance(in kilo ohms)
+disp(V_i_max,'Input microwave voltage(in volts) =');
+disp(A_v,'Voltage gain =');
+disp(eff,'Effeciency of amplifier (in percentage) =');
+disp(R_b,'Beam loading resistance(in kilo ohms) =');
+
+//Answer in book is wrongly given as: Voltage gain =17.034 \ No newline at end of file
diff --git a/728/CH8/EX8.5/Ex8_5.txt b/728/CH8/EX8.5/Ex8_5.txt
new file mode 100755
index 000000000..b7fa37ac5
--- /dev/null
+++ b/728/CH8/EX8.5/Ex8_5.txt
@@ -0,0 +1,28 @@
+//Caption:Calculate (i)-i/p microwave voltage,(ii)-voltage gain, (iii)-efficiency of amplifier,(iv)-beam loading conductance
+//Exa:8.5
+clc;
+clear;
+close;
+V_o=1200;//in volts
+I_o=28*10^-3;//in A
+f=8*10^9;//inHz
+d=0.001;//in m
+l=0.04;//in m
+R_sh=40*10^3;//in ohms
+V_p_max=1200*3.68*0.593*10^6*sqrt(V_o)/(2*%pi*f*l);
+Theeta_g=(2*%pi*f)*d/(0.593*10^6*sqrt(V_o));//transit angle (in rad)
+beeta=sin(Theeta_g/2)/(Theeta_g/2);
+V_i_max=V_p_max/beeta;
+Beeta_o=0.768;
+J=0.582;//J(X)
+A_v=(Beeta_o)^2*97.88*J*R_sh/(1200/(28*10^-3*1.841));//calculating voltage gain
+eff=[0.58*[2*28*10^-3*J*Beeta_o*R_sh]/V_o]*100;//calculating efficiency
+G_o=23.3*10^-6;
+G_b=(G_o/2)*{Beeta_o^2-Beeta_o*cos(Theeta_g)};//beam loading conductance
+R_b=1/(G_b*1000);//beam loading resistance(in kilo ohms)
+disp(V_i_max,'Input microwave voltage(in volts) =');
+disp(A_v,'Voltage gain =');
+disp(eff,'Effeciency of amplifier (in percentage) =');
+disp(R_b,'Beam loading resistance(in kilo ohms) =');
+
+//Answer in book is wrongly given as: Voltage gain =17.034 \ No newline at end of file
diff --git a/728/CH8/EX8.5/Ex8_5_ans.txt b/728/CH8/EX8.5/Ex8_5_ans.txt
new file mode 100755
index 000000000..71c170c00
--- /dev/null
+++ b/728/CH8/EX8.5/Ex8_5_ans.txt
@@ -0,0 +1,16 @@
+ Input microwave voltage(in volts) =
+
+ 58.70546
+
+ Voltage gain =
+
+ 57.733774
+
+ Effeciency of amplifier (in percentage) =
+
+ 48.392602
+
+ Beam loading resistance(in kilo ohms) =
+
+ 72.751565
+ \ No newline at end of file
diff --git a/728/CH8/EX8.6/Ex8_6.sce b/728/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..ce568ae31
--- /dev/null
+++ b/728/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,26 @@
+//Caption:Calculate (i)-value of repeller voltage V_r ,(ii)-beam current necessary to give gap voltage of 200V, (iii)-electronic efficiency
+//Exa:8.6
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+V_o=500;//in volts
+R_sh=20*10^3;//in ohms
+f=8*10^9;//inHz
+w=2*%pi*f;
+n=2;//mode
+L=0.001;//spacing between repeller & cavity (in m)
+x=0.023;
+volt_diff=sqrt(V_o*(x));
+V_r=volt_diff+V_o;//repeller volatge
+Beeta_o=1;//Assuming
+J=0.582;
+V_1=200;//given (in volts)
+I_o=V_1/(R_sh*2*J);
+Theeta_o=2*%pi*f*J*10^6*2*10^-3*sqrt(V_o)/(1.579*10^11*(V_r+V_o));
+X=V_1*Theeta_o/(2*V_o);//X'
+j=0.84;//J(X')
+eff=[2*j/(2*%pi*2-%pi/2)]*100;
+disp(V_r,'Repeller voltage(in volts) =');
+disp(I_o,'Necessary beam current (in Amp.s) =');
+disp(eff,'Effeciency (in percentage) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.6/Ex8_6.txt b/728/CH8/EX8.6/Ex8_6.txt
new file mode 100755
index 000000000..ce568ae31
--- /dev/null
+++ b/728/CH8/EX8.6/Ex8_6.txt
@@ -0,0 +1,26 @@
+//Caption:Calculate (i)-value of repeller voltage V_r ,(ii)-beam current necessary to give gap voltage of 200V, (iii)-electronic efficiency
+//Exa:8.6
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+V_o=500;//in volts
+R_sh=20*10^3;//in ohms
+f=8*10^9;//inHz
+w=2*%pi*f;
+n=2;//mode
+L=0.001;//spacing between repeller & cavity (in m)
+x=0.023;
+volt_diff=sqrt(V_o*(x));
+V_r=volt_diff+V_o;//repeller volatge
+Beeta_o=1;//Assuming
+J=0.582;
+V_1=200;//given (in volts)
+I_o=V_1/(R_sh*2*J);
+Theeta_o=2*%pi*f*J*10^6*2*10^-3*sqrt(V_o)/(1.579*10^11*(V_r+V_o));
+X=V_1*Theeta_o/(2*V_o);//X'
+j=0.84;//J(X')
+eff=[2*j/(2*%pi*2-%pi/2)]*100;
+disp(V_r,'Repeller voltage(in volts) =');
+disp(I_o,'Necessary beam current (in Amp.s) =');
+disp(eff,'Effeciency (in percentage) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.6/Ex8_6_ans.txt b/728/CH8/EX8.6/Ex8_6_ans.txt
new file mode 100755
index 000000000..e542c97c6
--- /dev/null
+++ b/728/CH8/EX8.6/Ex8_6_ans.txt
@@ -0,0 +1,12 @@
+Repeller voltage(in volts) =
+
+ 503.39116
+
+ Necessary beam current (in Amp.s) =
+
+ 0.0085911
+
+ Effeciency (in percentage) =
+
+ 15.278875
+ \ No newline at end of file
diff --git a/728/CH8/EX8.7/Ex8_7.sce b/728/CH8/EX8.7/Ex8_7.sce
new file mode 100755
index 000000000..25bf11dfb
--- /dev/null
+++ b/728/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-efficiency of reflex klystron ,(ii)-total output power in mW, (iii)-power delivered to load
+//Exa:8.7
+clc;
+clear;
+close;
+P_dc_in=40;//in mW
+ratio=0.278;//V_1/V_o;
+n=1;
+X=ratio*(2*n*%pi-%pi/2);
+J=2.35;
+eff=ratio*J*100;//in percentage
+P_out= 8.91*P_dc_in/100;
+P_load=3.564*80/100;
+disp(eff,'Effeciency (in percentage) =');
+disp(P_out,'Total power output (in mW) =');
+disp(P_load,'Power delivered to load (in mW) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.7/Ex8_7.txt b/728/CH8/EX8.7/Ex8_7.txt
new file mode 100755
index 000000000..25bf11dfb
--- /dev/null
+++ b/728/CH8/EX8.7/Ex8_7.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-efficiency of reflex klystron ,(ii)-total output power in mW, (iii)-power delivered to load
+//Exa:8.7
+clc;
+clear;
+close;
+P_dc_in=40;//in mW
+ratio=0.278;//V_1/V_o;
+n=1;
+X=ratio*(2*n*%pi-%pi/2);
+J=2.35;
+eff=ratio*J*100;//in percentage
+P_out= 8.91*P_dc_in/100;
+P_load=3.564*80/100;
+disp(eff,'Effeciency (in percentage) =');
+disp(P_out,'Total power output (in mW) =');
+disp(P_load,'Power delivered to load (in mW) ='); \ No newline at end of file
diff --git a/728/CH8/EX8.7/Ex8_7_ans.txt b/728/CH8/EX8.7/Ex8_7_ans.txt
new file mode 100755
index 000000000..4e47a845d
--- /dev/null
+++ b/728/CH8/EX8.7/Ex8_7_ans.txt
@@ -0,0 +1,12 @@
+ Effeciency (in percentage) =
+
+ 65.33
+
+ Total power output (in mW) =
+
+ 3.564
+
+ Power delivered to load (in mW) =
+
+ 2.8512
+ \ No newline at end of file
diff --git a/728/CH8/EX8.8/Ex8_8.sce b/728/CH8/EX8.8/Ex8_8.sce
new file mode 100755
index 000000000..47ced396b
--- /dev/null
+++ b/728/CH8/EX8.8/Ex8_8.sce
@@ -0,0 +1,20 @@
+//Caption:Calculate (i)-Hull cut-off voltage ,(ii)-cut-off magnetic flux density if beam voltage V_o is 6000V, (iii)-cyclotron frequency in GHz
+//Exa:8.8
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+R_a=0.15;//in m
+R_o=0.45;//in m
+B_o=1.2*10^-3;//in weber/m^2
+V_o={(e_m_ratio)*B_o^2*R_o^2*[1-(R_a/R_o)^2]^2}/8;
+//Given:
+V=6000;//in volts
+B_c=sqrt(8*V/e_m_ratio)/[[1-(R_a/R_o)^2]*(R_o)];//in weber/m^2
+w_c=(e_m_ratio)*B_o;
+f_c=w_c/(2*%pi);//in Hz
+disp(V_o,'Cut-off voltage (in volts) =');
+disp(B_c*10^5,'Cut-off magnetic flux density (in milli weber/sq. m) =');
+disp(f_c*10^-9,'Cyclotron frequency (in GHz) =');
+
+//Answer in book is wrongly given as: f_c=0.336Hz & V_o=50.666 kV \ No newline at end of file
diff --git a/728/CH8/EX8.8/Ex8_8.txt b/728/CH8/EX8.8/Ex8_8.txt
new file mode 100755
index 000000000..47ced396b
--- /dev/null
+++ b/728/CH8/EX8.8/Ex8_8.txt
@@ -0,0 +1,20 @@
+//Caption:Calculate (i)-Hull cut-off voltage ,(ii)-cut-off magnetic flux density if beam voltage V_o is 6000V, (iii)-cyclotron frequency in GHz
+//Exa:8.8
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+R_a=0.15;//in m
+R_o=0.45;//in m
+B_o=1.2*10^-3;//in weber/m^2
+V_o={(e_m_ratio)*B_o^2*R_o^2*[1-(R_a/R_o)^2]^2}/8;
+//Given:
+V=6000;//in volts
+B_c=sqrt(8*V/e_m_ratio)/[[1-(R_a/R_o)^2]*(R_o)];//in weber/m^2
+w_c=(e_m_ratio)*B_o;
+f_c=w_c/(2*%pi);//in Hz
+disp(V_o,'Cut-off voltage (in volts) =');
+disp(B_c*10^5,'Cut-off magnetic flux density (in milli weber/sq. m) =');
+disp(f_c*10^-9,'Cyclotron frequency (in GHz) =');
+
+//Answer in book is wrongly given as: f_c=0.336Hz & V_o=50.666 kV \ No newline at end of file
diff --git a/728/CH8/EX8.8/Ex8_8_ans.txt b/728/CH8/EX8.8/Ex8_8_ans.txt
new file mode 100755
index 000000000..88398e501
--- /dev/null
+++ b/728/CH8/EX8.8/Ex8_8_ans.txt
@@ -0,0 +1,13 @@
+
+ Cut-off voltage (in volts) =
+
+ 5065.92
+
+ Cut-off magnetic flux density (in milli weber/sq. m) =
+
+ 130.59535
+
+ Cyclotron frequency (in GHz) =
+
+ 0.0335944
+ \ No newline at end of file
diff --git a/728/CH8/EX8.9/EX8_9_ans.txt b/728/CH8/EX8.9/EX8_9_ans.txt
new file mode 100755
index 000000000..c89bf2e44
--- /dev/null
+++ b/728/CH8/EX8.9/EX8_9_ans.txt
@@ -0,0 +1,8 @@
+Axial phase velocity (in m/s) =
+
+ 9549296.6
+
+ Anode Voltage (in kV) =
+
+ 259.20712
+ \ No newline at end of file
diff --git a/728/CH8/EX8.9/Ex8_9.sce b/728/CH8/EX8.9/Ex8_9.sce
new file mode 100755
index 000000000..adb62dde3
--- /dev/null
+++ b/728/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-Axial phase velocity ,(ii)-Anode voltage at which TWT can be operated for useful gain
+//Exa:8.9
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+c=3*10^8;//in m/s
+d=0.002;//diameter(in m)
+pitch=(1/50)/100;//As,50 turns per cm (in m)
+circum=%pi*d;
+v_p=c*pitch/circum;
+V_o=v_p^2/(2*e_m_ratio);
+disp(v_p,'Axial phase velocity (in m/s) =');
+disp(V_o,'Anode Voltage (in kV) =');
+
+//Answer in book is wrongly given as V_o=25.92 V \ No newline at end of file
diff --git a/728/CH8/EX8.9/Ex8_9.txt b/728/CH8/EX8.9/Ex8_9.txt
new file mode 100755
index 000000000..adb62dde3
--- /dev/null
+++ b/728/CH8/EX8.9/Ex8_9.txt
@@ -0,0 +1,16 @@
+//Caption:Calculate (i)-Axial phase velocity ,(ii)-Anode voltage at which TWT can be operated for useful gain
+//Exa:8.9
+clc;
+clear;
+close;
+e_m_ratio=1.759*10^11;//(e/m)
+c=3*10^8;//in m/s
+d=0.002;//diameter(in m)
+pitch=(1/50)/100;//As,50 turns per cm (in m)
+circum=%pi*d;
+v_p=c*pitch/circum;
+V_o=v_p^2/(2*e_m_ratio);
+disp(v_p,'Axial phase velocity (in m/s) =');
+disp(V_o,'Anode Voltage (in kV) =');
+
+//Answer in book is wrongly given as V_o=25.92 V \ No newline at end of file
diff --git a/728/CH9/EX9.1/Ex9_1.sce b/728/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..ede2a92c1
--- /dev/null
+++ b/728/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,9 @@
+//Caption:Calculate operating frequency of IMPATT diode
+//Exa:9.1
+clc;
+clear;
+close;
+v_d=10^7*10^-2;//drift velocity(in m/s)
+L=2*10^-6;//drift length(in m)
+f=v_d/(2*L);//in Hz
+disp(f/10^9,'Operating Frequency (in GHz) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.1/Ex9_1.txt b/728/CH9/EX9.1/Ex9_1.txt
new file mode 100755
index 000000000..ede2a92c1
--- /dev/null
+++ b/728/CH9/EX9.1/Ex9_1.txt
@@ -0,0 +1,9 @@
+//Caption:Calculate operating frequency of IMPATT diode
+//Exa:9.1
+clc;
+clear;
+close;
+v_d=10^7*10^-2;//drift velocity(in m/s)
+L=2*10^-6;//drift length(in m)
+f=v_d/(2*L);//in Hz
+disp(f/10^9,'Operating Frequency (in GHz) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.1/Ex9_1_ans.txt b/728/CH9/EX9.1/Ex9_1_ans.txt
new file mode 100755
index 000000000..20ac8c613
--- /dev/null
+++ b/728/CH9/EX9.1/Ex9_1_ans.txt
@@ -0,0 +1,5 @@
+
+ Operating Frequency (in GHz) =
+
+ 25.
+ \ No newline at end of file
diff --git a/728/CH9/EX9.10/Ex9_10.sce b/728/CH9/EX9.10/Ex9_10.sce
new file mode 100755
index 000000000..29f96e369
--- /dev/null
+++ b/728/CH9/EX9.10/Ex9_10.sce
@@ -0,0 +1,11 @@
+//Caption:Calculate (i)-drift time of carrier ,(ii)-operating frequency of diode
+//Exa:9.10
+clc;
+clear;
+close;
+L=2*10^-6;//drift length (in m)
+v_d=10^5;//in cm/s
+drift_time=L/v_d;
+f=1/(2*drift_time);
+disp(drift_time,'Drift time (in sec) =');
+disp(f*10^-9,'Operating Frequency (in GHz)='); \ No newline at end of file
diff --git a/728/CH9/EX9.10/Ex9_10.txt b/728/CH9/EX9.10/Ex9_10.txt
new file mode 100755
index 000000000..29f96e369
--- /dev/null
+++ b/728/CH9/EX9.10/Ex9_10.txt
@@ -0,0 +1,11 @@
+//Caption:Calculate (i)-drift time of carrier ,(ii)-operating frequency of diode
+//Exa:9.10
+clc;
+clear;
+close;
+L=2*10^-6;//drift length (in m)
+v_d=10^5;//in cm/s
+drift_time=L/v_d;
+f=1/(2*drift_time);
+disp(drift_time,'Drift time (in sec) =');
+disp(f*10^-9,'Operating Frequency (in GHz)='); \ No newline at end of file
diff --git a/728/CH9/EX9.10/Ex9_10_ans.txt b/728/CH9/EX9.10/Ex9_10_ans.txt
new file mode 100755
index 000000000..6704d01cd
--- /dev/null
+++ b/728/CH9/EX9.10/Ex9_10_ans.txt
@@ -0,0 +1,9 @@
+
+ Drift time (in sec) =
+
+ 2.000D-11
+
+ Operating Frequency (in GHz)=
+
+ 25.
+ \ No newline at end of file
diff --git a/728/CH9/EX9.11/Ex9_11.sce b/728/CH9/EX9.11/Ex9_11.sce
new file mode 100755
index 000000000..690020d67
--- /dev/null
+++ b/728/CH9/EX9.11/Ex9_11.sce
@@ -0,0 +1,14 @@
+//Caption:Calculate (i)-breakdown voltage ,(ii)-breakdown electric field.
+//Exa:9.11
+clc;
+clear;
+close;
+E_r=11.8;
+E_o=8.85*10^-12;
+N=3*10^21;//in per cubic meter
+L=6.2*10^-6;//in meter
+q=1.6*10^-19;//in coulombs
+V_bd=q*N*L^2/(E_o*E_r);
+E_bd=V_bd/L;
+disp(V_bd,'Breakdown voltage (in volts) =');
+disp(E_bd,'Breakdown electric field (in V/m) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.11/Ex9_11.txt b/728/CH9/EX9.11/Ex9_11.txt
new file mode 100755
index 000000000..690020d67
--- /dev/null
+++ b/728/CH9/EX9.11/Ex9_11.txt
@@ -0,0 +1,14 @@
+//Caption:Calculate (i)-breakdown voltage ,(ii)-breakdown electric field.
+//Exa:9.11
+clc;
+clear;
+close;
+E_r=11.8;
+E_o=8.85*10^-12;
+N=3*10^21;//in per cubic meter
+L=6.2*10^-6;//in meter
+q=1.6*10^-19;//in coulombs
+V_bd=q*N*L^2/(E_o*E_r);
+E_bd=V_bd/L;
+disp(V_bd,'Breakdown voltage (in volts) =');
+disp(E_bd,'Breakdown electric field (in V/m) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.11/Ex9_11_ans.txt b/728/CH9/EX9.11/Ex9_11_ans.txt
new file mode 100755
index 000000000..4fe1ffb5a
--- /dev/null
+++ b/728/CH9/EX9.11/Ex9_11_ans.txt
@@ -0,0 +1,8 @@
+ Breakdown voltage (in volts) =
+
+ 176.68486
+
+ Breakdown electric field (in V/m) =
+
+ 28497558.
+ \ No newline at end of file
diff --git a/728/CH9/EX9.12/Ex9_12.sce b/728/CH9/EX9.12/Ex9_12.sce
new file mode 100755
index 000000000..4c97fc7e5
--- /dev/null
+++ b/728/CH9/EX9.12/Ex9_12.sce
@@ -0,0 +1,18 @@
+ //Caption:Calculate (i)-maximum power gain in dBs ,(ii)-noise figure F in dBs, (iii)-bandwidth for r=0.2
+//Exa:9.12
+clc;
+clear;
+close;
+ratio=8;
+r=0.2;
+r_Q=8;
+T_d=300;//in Kelvin
+T_o=300;//in Kelvin
+X=8;
+G=(ratio)*X/(1+sqrt(1+X))^2;
+G_in_dB=(10*log(G))/log(10);//gain
+disp(G_in_dB,'Maximum Gain (in dB)=');
+F=[10*log(1+(2*T_d/T_o)*[(1/(r_Q))+(1/(r_Q)^2)])]/log(10);//noise figure
+disp(F,'Noise figure (in dB) =');
+B_W=2*r*sqrt(ratio);//bandwidth
+disp(B_W,'bandwidth ='); \ No newline at end of file
diff --git a/728/CH9/EX9.12/Ex9_12.txt b/728/CH9/EX9.12/Ex9_12.txt
new file mode 100755
index 000000000..4c97fc7e5
--- /dev/null
+++ b/728/CH9/EX9.12/Ex9_12.txt
@@ -0,0 +1,18 @@
+ //Caption:Calculate (i)-maximum power gain in dBs ,(ii)-noise figure F in dBs, (iii)-bandwidth for r=0.2
+//Exa:9.12
+clc;
+clear;
+close;
+ratio=8;
+r=0.2;
+r_Q=8;
+T_d=300;//in Kelvin
+T_o=300;//in Kelvin
+X=8;
+G=(ratio)*X/(1+sqrt(1+X))^2;
+G_in_dB=(10*log(G))/log(10);//gain
+disp(G_in_dB,'Maximum Gain (in dB)=');
+F=[10*log(1+(2*T_d/T_o)*[(1/(r_Q))+(1/(r_Q)^2)])]/log(10);//noise figure
+disp(F,'Noise figure (in dB) =');
+B_W=2*r*sqrt(ratio);//bandwidth
+disp(B_W,'bandwidth ='); \ No newline at end of file
diff --git a/728/CH9/EX9.12/Ex9_12_ans.txt b/728/CH9/EX9.12/Ex9_12_ans.txt
new file mode 100755
index 000000000..9aac91011
--- /dev/null
+++ b/728/CH9/EX9.12/Ex9_12_ans.txt
@@ -0,0 +1,11 @@
+Maximum Gain (in dB)=
+
+ 6.0205999
+
+ Noise figure (in dB) =
+
+ 1.0763388
+
+ bandwidth =
+
+ 1.1313708 \ No newline at end of file
diff --git a/728/CH9/EX9.13/Ex9_13.sce b/728/CH9/EX9.13/Ex9_13.sce
new file mode 100755
index 000000000..b7d8dca6d
--- /dev/null
+++ b/728/CH9/EX9.13/Ex9_13.sce
@@ -0,0 +1,32 @@
+//Caption:Calculate (i)-equivalent noise resistance, (ii)-gain, (iii)-noise figure, (iv)-bandwidth
+//Exa:9.13
+clc;
+clear;
+close;
+f_s=2*10^9;//in Hz
+f_p=12*10^9;//in Hz
+f_i=10*10^9;//in Hz
+f_d=5*10^9;//in Hz
+R_i=1*10^3;//in ohm
+R_g=1*10^3;//in ohm
+R_T_s=1*10^3;//in ohm
+R_T_i=1*10^3;//in ohm
+T_d=300;//in Kelvin
+T_o=300;//in Kelvin
+w_s=2*%pi*f_s;
+w_i=2*%pi*f_i;
+r=0.35;
+r_Q=10;
+r_d=300;//in ohm
+C=0.01*10^-12;//in Farad
+R=r^2/(w_s*w_i*C^2*R_T_i);
+a=R/R_T_s;
+g=((4*f_i*R_g*R_i*a)/(f_s*R_T_s*R_T_i*(1-a)^2));//gain
+Gain=[10*log(g)]/log(10);//gain in dB
+f={1+((2*T_d)/T_o)*[(1/r_Q)+(1/r_Q^2)]};//noise figure
+F=[10*log(f)]/log(10);//noise figure in dB
+B_W=(r/2)*sqrt(f_d/(f_s*Gain));//bandwidth
+disp(a,'Equivalent noise resistance (in ohm) =');
+disp(Gain,'Gain (in dB) =');
+disp(F,'Noise figure (in dB) =');
+disp(B_W,'Bandwidth ='); \ No newline at end of file
diff --git a/728/CH9/EX9.13/Ex9_13.txt b/728/CH9/EX9.13/Ex9_13.txt
new file mode 100755
index 000000000..b7d8dca6d
--- /dev/null
+++ b/728/CH9/EX9.13/Ex9_13.txt
@@ -0,0 +1,32 @@
+//Caption:Calculate (i)-equivalent noise resistance, (ii)-gain, (iii)-noise figure, (iv)-bandwidth
+//Exa:9.13
+clc;
+clear;
+close;
+f_s=2*10^9;//in Hz
+f_p=12*10^9;//in Hz
+f_i=10*10^9;//in Hz
+f_d=5*10^9;//in Hz
+R_i=1*10^3;//in ohm
+R_g=1*10^3;//in ohm
+R_T_s=1*10^3;//in ohm
+R_T_i=1*10^3;//in ohm
+T_d=300;//in Kelvin
+T_o=300;//in Kelvin
+w_s=2*%pi*f_s;
+w_i=2*%pi*f_i;
+r=0.35;
+r_Q=10;
+r_d=300;//in ohm
+C=0.01*10^-12;//in Farad
+R=r^2/(w_s*w_i*C^2*R_T_i);
+a=R/R_T_s;
+g=((4*f_i*R_g*R_i*a)/(f_s*R_T_s*R_T_i*(1-a)^2));//gain
+Gain=[10*log(g)]/log(10);//gain in dB
+f={1+((2*T_d)/T_o)*[(1/r_Q)+(1/r_Q^2)]};//noise figure
+F=[10*log(f)]/log(10);//noise figure in dB
+B_W=(r/2)*sqrt(f_d/(f_s*Gain));//bandwidth
+disp(a,'Equivalent noise resistance (in ohm) =');
+disp(Gain,'Gain (in dB) =');
+disp(F,'Noise figure (in dB) =');
+disp(B_W,'Bandwidth ='); \ No newline at end of file
diff --git a/728/CH9/EX9.13/Ex9_13_ans.txt b/728/CH9/EX9.13/Ex9_13_ans.txt
new file mode 100755
index 000000000..fb2d03536
--- /dev/null
+++ b/728/CH9/EX9.13/Ex9_13_ans.txt
@@ -0,0 +1,16 @@
+
+ Equivalent noise resistance (in ohm) =
+
+ 1.5514806
+
+ Gain (in dB) =
+
+ 20.087158
+
+ Noise figure (in dB) =
+
+ 0.8635983
+
+ Bandwidth =
+
+ 0.0617375 \ No newline at end of file
diff --git a/728/CH9/EX9.2/Ex9_2.sce b/728/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..dda2350c7
--- /dev/null
+++ b/728/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,10 @@
+//Caption:Determine threshold electric field
+//Exa:9_2
+clc;
+clear;
+close;
+f=10*10^9;//in Hz
+L=75*10^-6;//Device length (in m)
+V=25;//Voltage pulse amplified (in volts)
+E_th=V/L;
+disp(E_th,'Threshold Electric field (in kV/cm) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.2/Ex9_2.txt b/728/CH9/EX9.2/Ex9_2.txt
new file mode 100755
index 000000000..dda2350c7
--- /dev/null
+++ b/728/CH9/EX9.2/Ex9_2.txt
@@ -0,0 +1,10 @@
+//Caption:Determine threshold electric field
+//Exa:9_2
+clc;
+clear;
+close;
+f=10*10^9;//in Hz
+L=75*10^-6;//Device length (in m)
+V=25;//Voltage pulse amplified (in volts)
+E_th=V/L;
+disp(E_th,'Threshold Electric field (in kV/cm) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.2/Ex9_2_ans.txt b/728/CH9/EX9.2/Ex9_2_ans.txt
new file mode 100755
index 000000000..ef873ecf1
--- /dev/null
+++ b/728/CH9/EX9.2/Ex9_2_ans.txt
@@ -0,0 +1,4 @@
+Threshold Electric field (in kV/cm) =
+
+ 333333.33
+ \ No newline at end of file
diff --git a/728/CH9/EX9.3/Ex9_3.sce b/728/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..89c25e195
--- /dev/null
+++ b/728/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,13 @@
+//Caption:Calculate (i)-power gain in dB ,(ii)-power gain if it is USB converter.
+//Exa:9.3
+clc;
+clear;
+close;
+f_s=2*10^9;//in Hz
+f_p=12*10^9;//in Hz
+R_i=16;
+R_s=1000;
+A_p=10*log((f_p-f_s)/f_s);
+A_p_usb=10*log((f_p+f_s)/f_s);
+disp(log(10),'Power gain (in dB) =');
+disp(A_p_usb,'Power gain as USB converter (in dB) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.3/Ex9_3.txt b/728/CH9/EX9.3/Ex9_3.txt
new file mode 100755
index 000000000..89c25e195
--- /dev/null
+++ b/728/CH9/EX9.3/Ex9_3.txt
@@ -0,0 +1,13 @@
+//Caption:Calculate (i)-power gain in dB ,(ii)-power gain if it is USB converter.
+//Exa:9.3
+clc;
+clear;
+close;
+f_s=2*10^9;//in Hz
+f_p=12*10^9;//in Hz
+R_i=16;
+R_s=1000;
+A_p=10*log((f_p-f_s)/f_s);
+A_p_usb=10*log((f_p+f_s)/f_s);
+disp(log(10),'Power gain (in dB) =');
+disp(A_p_usb,'Power gain as USB converter (in dB) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.3/Ex9_3_ans.txt b/728/CH9/EX9.3/Ex9_3_ans.txt
new file mode 100755
index 000000000..42e9bba62
--- /dev/null
+++ b/728/CH9/EX9.3/Ex9_3_ans.txt
@@ -0,0 +1,8 @@
+ Power gain (in dB) =
+
+ 2.3025851
+
+ Power gain as USB converter (in dB) =
+
+ 19.459101
+ \ No newline at end of file
diff --git a/728/CH9/EX9.4/Ex9_4.sce b/728/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..70fa56276
--- /dev/null
+++ b/728/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,17 @@
+//Caption:Calculate (i)-critical voltage ,(ii)-breakdown voltage, (iii)-breakdown electric field
+//Exa:9.4
+clc;
+clear;
+close;
+E_s=12.5;
+E_o=8.85*10^-12;
+E=E_o*E_s;
+N=3.2*10^22;//per cubic meter
+L=8*10^-6;//in m
+q=1.6*10^-19;//in coulombs
+V_c=q*N*L^2/(2*E);
+V_bd=2*V_c;
+E_bd=V_bd/L;
+disp(V_c/10^3,'Critical voltage(in kV) =');
+disp(V_bd/10^3,'Breakdown Voltage (in kV) =');
+disp(E_bd,'Breakdown Electric field (in V/cm) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.4/Ex9_4.txt b/728/CH9/EX9.4/Ex9_4.txt
new file mode 100755
index 000000000..70fa56276
--- /dev/null
+++ b/728/CH9/EX9.4/Ex9_4.txt
@@ -0,0 +1,17 @@
+//Caption:Calculate (i)-critical voltage ,(ii)-breakdown voltage, (iii)-breakdown electric field
+//Exa:9.4
+clc;
+clear;
+close;
+E_s=12.5;
+E_o=8.85*10^-12;
+E=E_o*E_s;
+N=3.2*10^22;//per cubic meter
+L=8*10^-6;//in m
+q=1.6*10^-19;//in coulombs
+V_c=q*N*L^2/(2*E);
+V_bd=2*V_c;
+E_bd=V_bd/L;
+disp(V_c/10^3,'Critical voltage(in kV) =');
+disp(V_bd/10^3,'Breakdown Voltage (in kV) =');
+disp(E_bd,'Breakdown Electric field (in V/cm) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.4/Ex9_4_ans.txt b/728/CH9/EX9.4/Ex9_4_ans.txt
new file mode 100755
index 000000000..5fab198df
--- /dev/null
+++ b/728/CH9/EX9.4/Ex9_4_ans.txt
@@ -0,0 +1,11 @@
+ Critical voltage(in kV) =
+
+ 1.4810395
+
+ Breakdown Voltage (in kV) =
+
+ 2.9620791
+
+ Breakdown Electric field (in V/cm) =
+
+ 3.703D+08 \ No newline at end of file
diff --git a/728/CH9/EX9.5/Ex9_5.sce b/728/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..5616bbe89
--- /dev/null
+++ b/728/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,10 @@
+//Caption:Calculate (i)-power gain in dB ,(ii)-power gain if it is USB converter.
+//Exa:9.5
+clc;
+clear;
+close;
+N_a=2.5*10^16;//per cubic cm
+J=33;//in kA/cm^2
+q=1.6*10^-19;
+V_z=J/(q*N_a);//Avalanche zone velocity (in cm/s)
+disp(V_z,'Avalanche zone velocity (in cm/s) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.5/Ex9_5.txt b/728/CH9/EX9.5/Ex9_5.txt
new file mode 100755
index 000000000..edc4504f9
--- /dev/null
+++ b/728/CH9/EX9.5/Ex9_5.txt
@@ -0,0 +1,10 @@
+//Caption:Calculate (i)-power gain in dB ,(ii)-power gain if it is USB converter.
+//Exa:9.5
+clc;
+clear;
+close;
+N_a-2.5*10^16;//per cubic cm
+J=33;//in kA/cm^2
+q=1.6*10^-19;
+V_z=J/(q*N_a);//Avalanche zone velocity (in cm/s)
+disp(V_z,'Avalanche zone velocity (in cm/s) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.5/Ex9_5_ans.txt b/728/CH9/EX9.5/Ex9_5_ans.txt
new file mode 100755
index 000000000..713d9cc1b
--- /dev/null
+++ b/728/CH9/EX9.5/Ex9_5_ans.txt
@@ -0,0 +1,4 @@
+ Avalanche zone velocity (in cm/s) =
+
+ 8250.
+ \ No newline at end of file
diff --git a/728/CH9/EX9.6/Ex9_6.sce b/728/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..c4d4b69e7
--- /dev/null
+++ b/728/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,9 @@
+//Caption:Calculate the power gain in dB
+//Exa:9.6
+clc;
+clear;
+close;
+R_neg=25;//in ohm
+R_load=50;//in ohm
+G={[- abs(R_neg)-R_load]/[- abs(R_neg)+R_load]}^2;
+disp(G,'Power gain ='); \ No newline at end of file
diff --git a/728/CH9/EX9.6/Ex9_6.txt b/728/CH9/EX9.6/Ex9_6.txt
new file mode 100755
index 000000000..c4d4b69e7
--- /dev/null
+++ b/728/CH9/EX9.6/Ex9_6.txt
@@ -0,0 +1,9 @@
+//Caption:Calculate the power gain in dB
+//Exa:9.6
+clc;
+clear;
+close;
+R_neg=25;//in ohm
+R_load=50;//in ohm
+G={[- abs(R_neg)-R_load]/[- abs(R_neg)+R_load]}^2;
+disp(G,'Power gain ='); \ No newline at end of file
diff --git a/728/CH9/EX9.6/Ex9_6_ans.txt b/728/CH9/EX9.6/Ex9_6_ans.txt
new file mode 100755
index 000000000..8609bb195
--- /dev/null
+++ b/728/CH9/EX9.6/Ex9_6_ans.txt
@@ -0,0 +1,5 @@
+
+ Power gain =
+
+ 9.
+ \ No newline at end of file
diff --git a/728/CH9/EX9.7/Ex9_7.sce b/728/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..d7209aec6
--- /dev/null
+++ b/728/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,9 @@
+//Caption:Calculate the minimum voltage needed to GUNN effect
+//Exa:9.7
+clc;
+clear;
+close;
+volt_grad=3.3*10^3;//voltage gradient
+L=5*10^-4;//in drift length
+V_min=volt_grad*L;//in volts
+disp(V_min,'Minimum voltage needed (in Volts) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.7/Ex9_7.txt b/728/CH9/EX9.7/Ex9_7.txt
new file mode 100755
index 000000000..d7209aec6
--- /dev/null
+++ b/728/CH9/EX9.7/Ex9_7.txt
@@ -0,0 +1,9 @@
+//Caption:Calculate the minimum voltage needed to GUNN effect
+//Exa:9.7
+clc;
+clear;
+close;
+volt_grad=3.3*10^3;//voltage gradient
+L=5*10^-4;//in drift length
+V_min=volt_grad*L;//in volts
+disp(V_min,'Minimum voltage needed (in Volts) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.7/Ex9_7_ans.txt b/728/CH9/EX9.7/Ex9_7_ans.txt
new file mode 100755
index 000000000..6cff07d63
--- /dev/null
+++ b/728/CH9/EX9.7/Ex9_7_ans.txt
@@ -0,0 +1,4 @@
+ Minimum voltage needed (in Volts) =
+
+ 1.65
+ \ No newline at end of file
diff --git a/728/CH9/EX9.8/EX9_8.txt b/728/CH9/EX9.8/EX9_8.txt
new file mode 100755
index 000000000..9fa5dfdf2
--- /dev/null
+++ b/728/CH9/EX9.8/EX9_8.txt
@@ -0,0 +1,12 @@
+//Caption:Calculate the rational frequency & critical velocity of diode.
+//Exa:9.8
+clc;
+clear;
+close;
+v_d=2*10^7;//in cm/s
+L=20*10^-4;//in cm
+f=v_d/L;
+disp(f*10^-9,'Natural frequency (in GHz) =');
+critical_field=3.3*10^3;
+V=L*critical_field;
+disp(V,'Critical voltage (in volts) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.8/Ex9_8.sce b/728/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..9fa5dfdf2
--- /dev/null
+++ b/728/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,12 @@
+//Caption:Calculate the rational frequency & critical velocity of diode.
+//Exa:9.8
+clc;
+clear;
+close;
+v_d=2*10^7;//in cm/s
+L=20*10^-4;//in cm
+f=v_d/L;
+disp(f*10^-9,'Natural frequency (in GHz) =');
+critical_field=3.3*10^3;
+V=L*critical_field;
+disp(V,'Critical voltage (in volts) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.8/Ex9_8_ans.txt b/728/CH9/EX9.8/Ex9_8_ans.txt
new file mode 100755
index 000000000..cd57710e1
--- /dev/null
+++ b/728/CH9/EX9.8/Ex9_8_ans.txt
@@ -0,0 +1,9 @@
+
+ Natural frequency (in GHz) =
+
+ 10.
+
+ Critical voltage (in volts) =
+
+ 6.6
+ \ No newline at end of file
diff --git a/728/CH9/EX9.9/Ex9_9.sce b/728/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..bb4f13e16
--- /dev/null
+++ b/728/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,15 @@
+//Caption:Calculate the resonant frequency & efficiency.
+//Exa:9.9
+clc;
+clear;
+close;
+L_p=0.5*10^-9;//in H
+C_j=0.5*10^-12;//in F
+V_bd=100;//breakdown voltage (in volts)
+I_bias=100*10^-3;//bias current(in A)
+I_rf_peak=0.8;
+R_l=2;
+f=1/(2*%pi*sqrt(L_p*C_j));
+eff={(0.5*I_rf_peak^2*R_l)/(V_bd*I_bias)}*100;
+disp(f*10^-9,'Resonant frequency (in GHz) =');
+disp(eff,'Efficiency (in percentage) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.9/Ex9_9.txt b/728/CH9/EX9.9/Ex9_9.txt
new file mode 100755
index 000000000..bb4f13e16
--- /dev/null
+++ b/728/CH9/EX9.9/Ex9_9.txt
@@ -0,0 +1,15 @@
+//Caption:Calculate the resonant frequency & efficiency.
+//Exa:9.9
+clc;
+clear;
+close;
+L_p=0.5*10^-9;//in H
+C_j=0.5*10^-12;//in F
+V_bd=100;//breakdown voltage (in volts)
+I_bias=100*10^-3;//bias current(in A)
+I_rf_peak=0.8;
+R_l=2;
+f=1/(2*%pi*sqrt(L_p*C_j));
+eff={(0.5*I_rf_peak^2*R_l)/(V_bd*I_bias)}*100;
+disp(f*10^-9,'Resonant frequency (in GHz) =');
+disp(eff,'Efficiency (in percentage) ='); \ No newline at end of file
diff --git a/728/CH9/EX9.9/Ex9_9_ans.txt b/728/CH9/EX9.9/Ex9_9_ans.txt
new file mode 100755
index 000000000..7bafbc460
--- /dev/null
+++ b/728/CH9/EX9.9/Ex9_9_ans.txt
@@ -0,0 +1,8 @@
+ Resonant frequency (in GHz) =
+
+ 10.065842
+
+ Efficiency (in percentage) =
+
+ 6.4
+ \ No newline at end of file