diff options
author | prashantsinalkar | 2017-10-10 12:38:01 +0530 |
---|---|---|
committer | prashantsinalkar | 2017-10-10 12:38:01 +0530 |
commit | f35ea80659b6a49d1bb2ce1d7d002583f3f40947 (patch) | |
tree | eb72842d800ac1233e9d890e020eac5fd41b0b1b /60 | |
parent | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (diff) | |
download | Scilab-TBC-Uploads-f35ea80659b6a49d1bb2ce1d7d002583f3f40947.tar.gz Scilab-TBC-Uploads-f35ea80659b6a49d1bb2ce1d7d002583f3f40947.tar.bz2 Scilab-TBC-Uploads-f35ea80659b6a49d1bb2ce1d7d002583f3f40947.zip |
updated the code
Diffstat (limited to '60')
-rwxr-xr-x | 60/CH3/EX3.10/ex_10.sce | 2 | ||||
-rwxr-xr-x | 60/CH3/EX3.2.a/ex_2_a.sce | 50 | ||||
-rwxr-xr-x | 60/CH4/EX4.21/ex_21.sce | 4 | ||||
-rwxr-xr-x | 60/CH6/EX6.12/ex_12.sce | 14 |
4 files changed, 29 insertions, 41 deletions
diff --git a/60/CH3/EX3.10/ex_10.sce b/60/CH3/EX3.10/ex_10.sce index 855fd40df..42ccfa1cc 100755 --- a/60/CH3/EX3.10/ex_10.sce +++ b/60/CH3/EX3.10/ex_10.sce @@ -1,6 +1,6 @@ //example(3.10) - c=[-3 1 0 1]] + c=[-3 1 0 1] p3=poly(c,'x','coeff') roots(p3) //here diff --git a/60/CH3/EX3.2.a/ex_2_a.sce b/60/CH3/EX3.2.a/ex_2_a.sce index 684d1b91e..44d062385 100755 --- a/60/CH3/EX3.2.a/ex_2_a.sce +++ b/60/CH3/EX3.2.a/ex_2_a.sce @@ -1,30 +1,22 @@ -//example(3.2a) - - -//bisection method - -deff('[y]=f(x)','y=x^3-x-1') -bisection(1.0,1.5,f) - - -//regula falsi method - -deff('[y]=f(x)','y=x^3-x-1') -regularfalsi(1.0,1.5,f) - -//secant method - -deff('[y]=f(x)','y=x^3-x-1') -secant(1.0,1.5,f) - - +// // Example(3.2a) +// +//// finding roots using bisection method +// +// +// deff('[y]=f(x)','y=x^3-x-1') +// bisection(0,10,f); +// +// +////regula falsi method +// +//deff('[y]=f(x)','y=x^3-x-1') +//regularfalsi(0,10,f); +// +// +// //newton rapson method - - -x=(0.5+1)/2 -deff('[y]=f(x)','y=x^3-x-1') -deff('[y]=g(x)','y=1-0.2*cos(x)') - -newton(x,f,g) - - +// +//x=(0+10)/2 +//deff('[y]=f(x)','x^3-x-1') +//deff ('[y]=g(x)','y=2*x^2 -1') +//x=newton(x,f,g);
\ No newline at end of file diff --git a/60/CH4/EX4.21/ex_21.sce b/60/CH4/EX4.21/ex_21.sce index aee0e9737..a1542a30f 100755 --- a/60/CH4/EX4.21/ex_21.sce +++ b/60/CH4/EX4.21/ex_21.sce @@ -4,5 +4,5 @@ B=[1 2 0;2 1 0;0 0 -1] I=[1 0 0;0 1 0;0 0 1] //here we are taking lamda=a //det(B-a)*I=0 is characteristic equation to get lamda -deff('[y]=p(a)','p(a)=det(B-a)*I') -p(a)=0
\ No newline at end of file +deff('[y]=pa','pa=det(B-a)*I') +pa=0
\ No newline at end of file diff --git a/60/CH6/EX6.12/ex_12.sce b/60/CH6/EX6.12/ex_12.sce index 7eaa16436..5020647cb 100755 --- a/60/CH6/EX6.12/ex_12.sce +++ b/60/CH6/EX6.12/ex_12.sce @@ -1,15 +1,12 @@ - - - -/Example 6.12 +//Example 6.12 //Least squares approximation deff('[y]=f(x)','y=10-2*x+((x^2)/10)') -//we seek the polynomial of degree <= 2 which minimizes +//we seek the polynomial of degree <= 2 which minimizes //sum(n=1 to 9)[fn-p(xn)]^2 //we are dealing with scalar product with w(x)=1 -P0(x)=1 +P0x=1 //hence s0=0; B=0; @@ -36,7 +33,7 @@ x=poly(0,"x") y2=((x-B0)^2)-0.1166667 //similarly calculate s2 s2=0.05973332 -//p*(x)=(d0*)*P0(x)+(d1*)*P1(x)+(d2*)*P2(x) +//p*(x)=(d0*)*P0x+(d1*)*P1(x)+(d2*)*P2(x) //d0*=d0,d1*=d1,d2*=d2 are least squares approximation //d0*=d0=sigma(n=1 to 6)[fn/6] where fn=f(xn) @@ -49,5 +46,4 @@ p=d0+d1*(x-B0)+d2*{[(x-B0)^2]-C1} //c1=c1* ,c2=c2*,c3=c3* c1=9.99998 c2=-1.9999998 -c3=0.0999999 -
\ No newline at end of file +c3=0.0999999
\ No newline at end of file |