CONTOH PENYELESAIAN KASUS 3 MENGGUNAKAN METODE EULER

KASUS PERSAMAAN 

LISTING PROGRAMNYA ADALAH

 

clear all

clc

%program penyelesaian persamaan diferensial

Q=10;

In=0.0005;

n=0.025;

g=9.81;

B=5;

disp (' ')

ya=input('kedalaman air diujung hilir = ');

xa=0;

xt=input(' masukan nilai x akhir     = ');

N=input ('masukan banyak iterasi     = ');

dx=(xt-xa)/N

X=zeros(N,1);

A=X;

P=X;

R=X;

T=X;

Y=X;

A(1)=B*ya;

P(1)=B+2*ya;

R(1)=A(1)/P(1);

T(1)=B;

Y(1)=ya+((In-((n)^2*(Q)^2)/((A(1))^2)*(R(1))^(4/3))/(1-((Q)^2*T(1))/(g*(A(1))^3)))*dx;

for s=2:1:N;

A(s)=B*Y(s-1);

P(s)=B+2*Y(s-1);

R(s)=A(s)/P(s);

T(s)=B;

Y(s)=Y(s-1)+((In-((n)^2*(Q)^2)/(((A(s-1))^2)*(R(s-1))^(4/3)))/(1-((Q)^2*T(s-1))/(g*(A(s-1))^3)))*dx;

end

d=1:1:N;

dd=d';

disp (' ')

disp (' ')

disp (' ')

disp ('       no.         Yi+1')

disp ([ dd  Y ])