KASUS PERSAMAAN DIFERENSIAL
PENYELESAIAN MENGGUNAKAN MANUAL TERLEBIH DARULU MENJADI PERSAMAAN UMUM
clear all
clc
%program penyelesaian persamaan diferensial
ya=input(' masukan nilai y awal (0) = ');
xa=input(' masukan nilai x awal
= ');
xt=input(' masukan nilai x akhir
= ');
N=input (' masukan banyak iterasi
= ');
dx=(xt-xa)/N;
X=zeros(N,1);
Y=X;
Xi=X;
Y(1)=ya+xa*((ya)^(1/2))*dx;
Xi(1)=xa+dx;
for s=2:1:N;
Xi(s)=(s)*dx+xa;
Y(s)=Y(s-1)+(Xi(s-1))*((Y(s-1))^(1/2))*dx;
end
d=1:1:N;
dd=d';
disp (' no. Xi+1 Yi+1')
disp ([ dd Xi Y ])
Komentar
Posting Komentar