Komputasi untuk Sains dan Teknik - Universitas Indonesia

8.3. PDP PARABOLIK 14117 w0(i,1)=suhu(i);18 end1920 AA=[ (1+2*lambda) -lambda 0 0 0 0 0 0 0;21 -lambda (1+2*lambda) -lambda 0 0 0 0 0 0;22 0 -lambda (1+2*lambda) -lambda 0 0 0 0 0 ;23 0 0 -lambda (1+2*lambda) -lambda 0 0 0 0;24 0 0 0 -lambda (1+2*lambda) -lambda 0 0 0;25 0 0 0 0 -lambda (1+2*lambda) -lambda 0 0;26 0 0 0 0 0 -lambda (1+2*lambda) -lambda 0 ;27 0 0 0 0 0 0 -lambda (1+2*lambda) -lambda ;28 0 0 0 0 0 0 0 -lambda (1+2*lambda) ];2930 iterasi=50;31 for i=1:iterasi32 %&&&&&& Proses Eliminasi Gauss &&&&&&&&&&&&&&&&&&&&&&&&&33 A=AA; %Matriks Backward Difference dicopy supaya fix3435 for i=1:n36 A(i,n+1)=w0(i,1);37 end3839 %---------Proses Triangularisasi-----------40 for j=1:(n-1)4142 %----mulai proses pivot---43 if (A(j,j)==0)44 for p=1:n+145 u=A(j,p);46 v=A(j+1,p);47 A(j+1,p)=u;48 A(j,p)=v;49 end50 end51 %----akhir proses pivot---52 jj=j+1;53 for i=jj:n54 m=A(i,j)/A(j,j);55 for k=1:(n+1)56 A(i,k)=A(i,k)-(m*A(j,k));57 end58 end59 end60 %-------------------------------------------6162 %------Proses Substitusi mundur-------------63 w(n,1)=A(n,n+1)/A(n,n);6465 for i=n-1:-1:166 S=0;67 for j=n:-1:i+168 S=S+A(i,j)*w(j,1);69 end70 w(i,1)=(A(i,n+1)-S)/A(i,i);71 end72 %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&73 w0=w;74 end75 w