function​ nonint

clear ​all
clc
h = 0.1;
sc = 15;
n = sc/h;
h1(1) = 4.5;
h2(1) = 4;
t(1)= 0;
for​ i = 1:n
t(i+1) = t(i)+h;
p1=h*nonint1(t(i),h1(i));
p2=h*nonint1(t(i)+(h/2),h1(i)+(p1/2));
p3=h*nonint1(t(i)+(h/2),h1(i)+(p2/2));
p4=h*nonint1(t(i)+h,h1(i));
h1(i+1)=h1(i)+((p1+(2*p2)+(2*p3)+p4)/6);
m1=h*nonint2(t(i),h1(i),h2(i));
m2=h*nonint2(t(i)+(h/2),h1(i)+(p1/2),h2(i)+(m1/2));
m3=h*nonint2(t(i)+(h/2),h1(i)+(p2/2),h2(i)+(m2/2));
m4=h*nonint2(t(i)+h,h1(i)+p3,h2(i)+m3);
h2(i+1)=h2(i)+((m1+(2*m2)+(2*m3)+m4)/6);
end
disp(​' Time,min Height 1,ft Height 2,ft '​)
disp([t;h1;h2]')
subplot(2,1,1);plot(t,h1);
title(​'Variation of the liquid levels w.r.t time in the two tanks non
interacting system using Runge Kutta method'​)
xlabel(​'Time,min'​)
ylabel(​'Height 1,ft'​)
subplot(2,1,2);plot(t,h2);
xlabel(​'Time,min'​)
ylabel(​'Height 2,ft'​)
function​ adot = nonint1(t,h1)
q = 5;
r1 = 1;
q1 = (h1)/r1;
dh1dt = q-q1;
adot =[dh1dt];
function​ cdot = nonint2(t,h1,h2)
r1 = 1;
q1 = (h1)/r1;
r2 = 1;
q2 = (h2/r2);
dh2dt = q1-q2;
cdot = [dh2dt];
Time,min Height 1,ft Height 2,ft

0 4.5000 4.0000

0.1000 4.5484 4.0499

0.2000 4.5921 4.0995

0.3000 4.6315 4.1483

0.4000 4.6672 4.1960

0.5000 4.6994 4.2424

0.6000 4.7285 4.2873

0.7000 4.7547 4.3305

0.8000 4.7785 4.3720

0.9000 4.7999 4.4118

1.0000 4.8193 4.4496

1.1000 4.8367 4.4856

1.2000 4.8525 4.5198

1.3000 4.8668 4.5522

1.4000 4.8797 4.5827

1.5000 4.8913 4.6116

1.6000 4.9018 4.6387

1.7000 4.9113 4.6642

1.8000 4.9199 4.6881

1.9000 4.9277 4.7106

2.0000 4.9347 4.7316

2.1000 4.9410 4.7512

2.2000 4.9467 4.7695

2.3000 4.9519 4.7866

2.4000 4.9565 4.8026

2.5000 4.9607 4.8174

2.6000 4.9645 4.8313

2.8000 4.9711 4.8560

2.9000 4.9739 4.8671

3.0000 4.9764 4.8774

3.1000 4.9787 4.8869

3.2000 4.9807 4.8958

3.3000 4.9826 4.9039

3.4000 4.9843 4.9115

3.5000 4.9858 4.9185

3.6000 4.9872 4.9250

3.7000 4.9884 4.9310

3.8000 4.9895 4.9365

3.9000 4.9905 4.9416

4.0000 4.9915 4.9463

4.1000 4.9923 4.9506

4.2000 4.9930 4.9546

4.3000 4.9937 4.9583

4.4000 4.9943 4.9617

4.5000 4.9949 4.9648

4.6000 4.9954 4.9677

4.7000 4.9958 4.9704

4.8000 4.9962 4.9728

4.9000 4.9966 4.9751

5.0000 4.9969 4.9771

5.1000 4.9972 4.9790

5.2000 4.9975 4.9808

5.3000 4.9977 4.9824

5.4000 4.9979 4.9838

5.6000 4.9983 4.9864

5.7000 4.9985 4.9876

5.8000 4.9986 4.9886

5.9000 4.9988 4.9896

6.0000 4.9989 4.9905

6.1000 4.9990 4.9913

6.2000 4.9991 4.9920

6.3000 4.9992 4.9927

6.4000 4.9993 4.9933

6.5000 4.9993 4.9939

6.6000 4.9994 4.9944

6.7000 4.9995 4.9949

6.8000 4.9995 4.9953

6.9000 4.9996 4.9957

7.0000 4.9996 4.9961

7.1000 4.9996 4.9964

7.2000 4.9997 4.9967

7.3000 4.9997 4.9970

7.4000 4.9997 4.9973

7.5000 4.9998 4.9975

7.6000 4.9998 4.9977

7.7000 4.9998 4.9979

7.8000 4.9998 4.9981

7.9000 4.9998 4.9983

8.0000 4.9999 4.9984

8.1000 4.9999 4.9985

8.2000 4.9999 4.9987

8.4000 4.9999 4.9989

8.5000 4.9999 4.9990

8.6000 4.9999 4.9991

8.7000 4.9999 4.9992

8.8000 4.9999 4.9992

8.9000 4.9999 4.9993

9.0000 4.9999 4.9994

9.1000 5.0000 4.9994

9.2000 5.0000 4.9995

9.3000 5.0000 4.9995

9.4000 5.0000 4.9996

9.5000 5.0000 4.9996

9.6000 5.0000 4.9996

9.7000 5.0000 4.9997

9.8000 5.0000 4.9997

9.9000 5.0000 4.9997

10.0000 5.0000 4.9997

10.1000 5.0000 4.9998

10.2000 5.0000 4.9998

10.3000 5.0000 4.9998

10.4000 5.0000 4.9998

10.5000 5.0000 4.9998

10.6000 5.0000 4.9999

10.7000 5.0000 4.9999

10.8000 5.0000 4.9999

10.9000 5.0000 4.9999

11.0000 5.0000 4.9999

11.2000 5.0000 4.9999

11.3000 5.0000 4.9999

11.4000 5.0000 4.9999

11.5000 5.0000 4.9999

11.6000 5.0000 4.9999

11.7000 5.0000 4.9999

11.8000 5.0000 5.0000

11.9000 5.0000 5.0000

12.0000 5.0000 5.0000

12.1000 5.0000 5.0000

12.2000 5.0000 5.0000

12.3000 5.0000 5.0000

12.4000 5.0000 5.0000

12.5000 5.0000 5.0000

12.6000 5.0000 5.0000

12.7000 5.0000 5.0000

12.8000 5.0000 5.0000

12.9000 5.0000 5.0000

13.0000 5.0000 5.0000

13.1000 5.0000 5.0000

13.2000 5.0000 5.0000

13.3000 5.0000 5.0000

13.4000 5.0000 5.0000

13.5000 5.0000 5.0000

13.6000 5.0000 5.0000

13.7000 5.0000 5.0000

13.8000 5.0000 5.0000

14.0000 5.0000 5.0000

14.1000 5.0000 5.0000

14.2000 5.0000 5.0000

14.3000 5.0000 5.0000

14.4000 5.0000 5.0000

14.5000 5.0000 5.0000

14.6000 5.0000 5.0000

14.7000 5.0000 5.0000

14.8000 5.0000 5.0000

14.9000 5.0000 5.0000

15.0000 5.0000 5.0000

>>