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 --- 26/CH1/EX1.2.34/1_2_34.sce | 48 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 48 insertions(+) create mode 100755 26/CH1/EX1.2.34/1_2_34.sce (limited to '26/CH1/EX1.2.34/1_2_34.sce') diff --git a/26/CH1/EX1.2.34/1_2_34.sce b/26/CH1/EX1.2.34/1_2_34.sce new file mode 100755 index 000000000..cd7923dc5 --- /dev/null +++ b/26/CH1/EX1.2.34/1_2_34.sce @@ -0,0 +1,48 @@ +disp('the augmented matrix is:') +a=[1 0 0 0 0 0 0;1 2 4 8 16 32 2.9;1 4 16 64 256 1024 14.8;1 6 36 216 1296 7776 39.6;1 8 64 512 4096 32768 74.3;1 10 10^2 10^3 10^4 10^5 119]; +disp(a) +disp('performing row transformations') +for k=2:6 + a(k,:)=a(k,:)-a(1,:) +end +disp(a) +j=2; +for k=3:6 + a(k,:)=a(k,:)-j*a(2,:) + j=j+1; +end +disp(a) +j=[0 0 0 3 6 10]; +for k=4:6 + a(k,:)=a(k,:)-j(k)*a(3,:) +end +disp(a) +a(5,:)=a(5,:)-4*a(4,:) +a(6,:)=a(6,:)-10*a(4,:) +disp(a) +a(6,:)=a(6,:)-5*a(5,:) +disp(a) +a(6,:)=a(6,:)/a(6,6) +disp(a) +j=[0 32 960 4800 7680] +for k=1:5 + a(k,:)=a(k,:)-j(k)*a(6,:) +end +disp(a) +a(5,:)=a(5,:)/a(5,5) +j=[0 16 224 576] +for k=2:4 + a(k,:)=a(k,:)-j(k)*a(5,:) +end +a(4,:)=a(4,:)/48 +a(2,:)=a(2,:)-8*a(4,:) +a(3,:)=a(3,:)-48*a(4,:) +a(3,:)=a(3,:)/8 +a(2,:)=a(2,:)-4*a(3,:) +a(2,:)=a(2,:)/2 +disp(a) +v=[a(1,7) a(2,7) a(3,7) a(4,7) a(5,7) a(6,7)] +p=poly(v,"t","coeff") +disp('p(t)=') +disp(p) +disp('p(7.5)=64.6 hundred lb') \ No newline at end of file -- cgit