diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /260/CH4/EX4.21 | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '260/CH4/EX4.21')
| -rw-r--r-- | 260/CH4/EX4.21/4_21.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/260/CH4/EX4.21/4_21.sce b/260/CH4/EX4.21/4_21.sce new file mode 100644 index 000000000..8f9ea9d5c --- /dev/null +++ b/260/CH4/EX4.21/4_21.sce @@ -0,0 +1,28 @@ +//Eg-4.21
+//pg-202
+
+clear
+clc
+
+
+//After complex root substitution a+ib
+
+deff('[z]=u(a,b)','z=a^2-b^2-2*a+2');
+deff('[z]=v(a,b)','z=2*a*b-2*b');
+deff('[z]=ua(a)','z=2*a-2');
+deff('[z]=ub(b)','z=-2*b');
+ain=0.5;
+bin=.5;
+
+for i=1:6
+ U=feval(ain,bin,u);
+ V=feval(ain,bin,v);
+ UA=feval(ain,ua);
+ UB=feval(bin,ub);
+ anew=ain+(V*UB-U*UA)/(UA^2+UB^2);
+ bnew=bin+(V*UA-U*UB)/(UA^2+UB^2);
+ ain=anew;
+ bin=bnew;
+end
+
+printf('\nThe complex root is %f + i(%f) \n',ain,bin)
|