From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1332/CH21/EX21.2/21_2.sce | 36 ++++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100755 1332/CH21/EX21.2/21_2.sce (limited to '1332/CH21/EX21.2/21_2.sce') diff --git a/1332/CH21/EX21.2/21_2.sce b/1332/CH21/EX21.2/21_2.sce new file mode 100755 index 000000000..d9ba4eb82 --- /dev/null +++ b/1332/CH21/EX21.2/21_2.sce @@ -0,0 +1,36 @@ +//Example 21.2 +//Lagrange Interpolation in Parallel Computing +//Page no. 723 +clc;close;clear; + +xi=[-1,0,2,5]; +yi=[9,5,3,15]; +s=["x=1","n=4","Data:","(-1,9)","(0,5)","(2,3)","(5,15)"] +for i=1:4 + printf('\tProcessor\t') +end +printf('\n') +for i=1:4 + printf('\t N%i\t\t',i) +end +printf('\n') +for i=1:7 + for j=1:4 + printf(' %s\t\t',s(i)) + end + printf('\n') +end + +x=1;T=0; +for k=0:3 + p=yi(k+1) + for j=0:3 + if(j~=k) + p=p*((x-xi(j+1))/(xi(k+1)-xi(j+1))) + end + end + T=T+p; + printf('\nT(%i) = %g',k+1,p) +end +printf('\n\nT = %g',T) + -- cgit