close all;
clear all;
clc;

%initialize test problem
c = [-300; -500; 0; 0; 0];
A = [1 2 1 0 0;
     1 1 0 1 0;
     0 3 0 0 1];
b = [170; 150; 180]; 

%(arbitrary but feasible) initial values
x      = [10 10 10 10 10]';
lambda = [10 10 10]';
s      = [10 10 10 10 10]';

%specify algorithm paramters
tol = 1e-8;
sigma = 0.5;
alpha = 0.75;

%solve linear program
[x,lambda,s,fx,nit] = interior_point(A,b,c,x,lambda,s,tol,sigma,alpha)