% ADAMS-FOURTH ORDER PREDICTOR-CORRECTOR ALGORITHM 5.4
%
% To approximate the solution of the initial value problem
% y' = f(t,y), a <= t <= b, y(a) = alpha,
% at N+1 equally spaced points in the interval [a,b].
%
% INPUT: endpoints a,b; initial condition alpha; integer N.
%
% OUTPUT: approximation w to y at the (N+1) values of t.
syms('F', 'OK', 'A', 'B', 'ALPHA', 'N', 'FLAG', 'NAME', 'OUP');
syms('H', 'T', 'W', 'I', 'K1', 'K2', 'K3', 'K4', 'T0', 'W0', 'J');
syms('t','y', 's','Part1','Part2');
TRUE = 1;
FALSE = 0;
T = zeros(1,4);
W = zeros(1,4);
fprintf(1,'This is Adams-Bashforth Predictor Corrector Method\n');
fprintf(1,'Input the function F(t,y) in terms of t and y\n');
fprintf(1,'For example: y-t^2+1 \n');
s = input(' ','s');
F = inline(s,'t','y');
OK = FALSE;
while OK == FALSE
fprintf(1,'Input left and right endpoints on separate lines.\n');
A = input(' ');
B = input(' ');
if A >= B
fprintf(1,'Left endpoint must be less than right endpoint\n');
else
OK = TRUE;
end;
end;
fprintf(1,'Input the initial condition\n');
ALPHA = input(' ');
OK = FALSE;
while OK == FALSE
fprintf(1,'Input an integer > 3 for the number of subintervals\n');
N = input(' ');
if N <= 3
fprintf(1,'Number must be at least 4.\n');
else
OK = TRUE;
end;
end;
if OK == TRUE
fprintf(1,'Choice of output method:\n');
fprintf(1,'1. Output to screen\n');
fprintf(1,'2. Output to text file\n');
fprintf(1,'Please enter 1 or 2\n');
FLAG = input(' ');
if FLAG == 2
fprintf(1,'Input the file name in the form - drive:\\name.ext\n');
fprintf(1,'For example A:\\OUTPUT.DTA\n');
NAME = input(' ','s');
docsity.com
