Question

CAN YOU PLEASE CHECK WHATS WRONG?

function[F]=func(t,h)
global Atank Ahole kOne kTwo g rho
%F=zeros(2,1)
%NL=y(1)
%NG=y(2)
F=(kOne+kTwo*cos(2*pi*t/24)-rho*Ahole*sqrt(2*g*h))/(rho*Atank);

close all
clear all
%clear all
%close all
global Atank Ahole kOne kTwo g rho

Atank = 3.13;% Atank value
Ahole =0.06; %Ahole value
kOne = 300; % K1 value
kTwo = 200; % K2 value
g = 9.81; % Gravity
rho = 1000; % rho value

t0=0; tf=150; %range of time

%Initial value

h0_2=h(1);
h0_1=h(2);
h0_0=h(3);

n=input('Enter the number of steps:');
dt=(tf-t0)/n; %step size
t=[t0:dt:tf]; %vector of independent variabll (The time)

h0=3 %Initial Height
h=[h0_2; h0_1; h0_0; zeros(n-2,1)];

%Solution (Using Third-order Adam-Bashford Method)

for i=3:n
h(i+1)=h(i)+(dt/12)*(23*func(t(i),h(i))-16*func(t(i-1),h(i-1))+5*func(t(i-2),h(i-2)));
end
plot(t,h,'r')
title('Plot Employing Adams-Bashforth Method');
xlabel('time')
ylabel('Level of Water')

Answer 1

Hey,

Note: Brother function file are not meant to be run directly they are either run through the command window or some script and should be saved with the same name mentioned below. In case of any queries brother just comment in box I would be very happy to assist all your queries.  

========================================

Save file as func.m

========================================

function[F]=func(t,h)
global Atank Ahole kOne kTwo g rho
%F=zeros(2,1)
%NL=y(1)
%NG=y(2)
F=(kOne+kTwo*cos(2*pi*t/24)-rho*Ahole*sqrt(2*g*h))/(rho*Atank);

=======================================

Executable file

=======================================

clc;

close all
clear all
%clear all
%close all
global Atank Ahole kOne kTwo g rho

Atank = 3.13;% Atank value
Ahole =0.06; %Ahole value
kOne = 300; % K1 value
kTwo = 200; % K2 value
g = 9.81; % Gravity
rho = 1000; % rho value

t0=0; tf=150; %range of time

%Initial value

h0_2=0.01;
h0_1=0.02;
h0_0=0.03;

n=input('Enter the number of steps:');
dt=(tf-t0)/n; %step size
t=[t0:dt:tf]; %vector of independent variabll (The time)

h0=3 %Initial Height
h=[h0_2; h0_1; h0_0; zeros(n-2,1)];

%Solution (Using Third-order Adam-Bashford Method)

for i=3:n
h(i+1)=h(i)+(dt/12)*(23*func(t(i),h(i))-16*func(t(i-1),h(i-1))+5*func(t(i-2),h(i-2)));
end
plot(t,h,'r')
title('Plot Employing Adams-Bashforth Method');
xlabel('time')
ylabel('Level of Water')

Kindly revert for any queries

Thanks.