1 files changed, 25 insertions, 16 deletions
diff --git a/851/CH5/EX5.13/Figure5_13.sce b/851/CH5/EX5.13/Figure5_13.sce
index 00b537f8a..61a4c6180 100755
--- a/851/CH5/EX5.13/Figure5_13.sce
+++ b/851/CH5/EX5.13/Figure5_13.sce
@@ -1,16 +1,25 @@
-//clear//
-//Caption:u-Law companding
-//Figure5.13(a)Mulaw companding Nonlinear Quantization
-//Plotting mulaw characteristics for different 
-//Values of mu
-clc;
-x = 0:0.01:1; //Normalized input
-mu = [0,5,255];//different values of mu
-for i = 1:length(mu)
-  [Cx(i,:),Xmax(i)] =  mulaw(x,mu(i));
-end
-plot2d(x/Xmax(1),Cx(1,:),2)
-plot2d(x/Xmax(2),Cx(2,:),4)
-plot2d(x/Xmax(3),Cx(3,:),6)
-xtitle('Compression Law: u-Law companding','Normalized Input |x|','Normalized Output |c(x)|');
-legend(['u =0'],['u=5'],['u=255'])
+//clear//
+//Caption:u-Law companding
+//Figure5.13(a)Mulaw companding Nonlinear Quantization
+//Plotting mulaw characteristics for different 
+//Values of mu
+clc;
+ [Cx,Xmax] =  mulaw(x,mu)
+  Xmax  = max(abs(x));
+  if(log(1+mu)~=0)
+    Cx = (log(1+mu*abs(x/Xmax))./log(1+mu));
+  else
+    Cx = x/Xmax;
+  end
+Cx = Cx/Xmax; //normalization of output vector
+endfunction
+x = 0:0.01:1; //Normalized input
+mu = [0,5,255];//different values of mu
+for i = 1:length(mu)
+  [Cx(i,:),Xmax(i)] =  mulaw(x,mu(i));
+end
+plot2d(x/Xmax(1),Cx(1,:),2)
+plot2d(x/Xmax(2),Cx(2,:),4)
+plot2d(x/Xmax(3),Cx(3,:),6)
+xtitle('Compression Law: u-Law companding','Normalized Input |x|','Normalized Output |c(x)|');
+legend(['u =0'],['u=5'],['u=255'])
+\ No newline at end of file