Komputasi untuk Sains dan Teknik - Universitas Indonesia

8.3. PDP PARABOLIK 137sampai 1000 kali, maka indeks j bergerak dari 1 sampai 1000. Dengan bantuan script Matlab,proses perhitungan menjadi sangat singkat.8.3.2.1 Script Forward-DifferenceScript matlab Forward-Difference untuk menyelesaikan contoh masalah ini, dimana h = 0, 1 dank = 0, 00051 clear all2 clc34 n=9;5 alpha=1.0;6 k=0.0005;7 h=0.1;8 lambda=(alpha^2)*k/(h^2);910 % Kondisi awal11 for i=1:n12 suhu(i)=sin(pi*i*0.1);13 end1415 %Mengcopy kondisi awal ke w16 for i=1:n17 w0(i,1)=suhu(i);18 end1920 A=[ (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=1000;31 for k=1:iterasi32 disp(’perkalian matriks’)33 %======================================34 for i=1:n35 w(i,1)=0.0;36 end3738 for i=1:n39 for j=1:n40 w(i,1)=w(i,1)+A(i,j)*w0(j,1);41 end42 end43 %====================================44 w45 w0=w;46 end