From a23ba49056a36aff221d0c23e228b981222c3576 Mon Sep 17 00:00:00 2001
From: Yuan Xu <yuanxu@student.ethz.ch>
Date: Wed, 24 May 2023 17:35:46 +0200
Subject: add submission zip

---
 Submission files/chebycenter.m | 20 ++++++++++++++++++++
 1 file changed, 20 insertions(+)
 create mode 100644 Submission files/chebycenter.m

(limited to 'Submission files/chebycenter.m')

diff --git a/Submission files/chebycenter.m b/Submission files/chebycenter.m
new file mode 100644
index 0000000..cbdce1c
--- /dev/null
+++ b/Submission files/chebycenter.m	
@@ -0,0 +1,20 @@
+function [c,r] = chebycenter(A,b)
+%CHEBYCENTER Compute Chebyshev center of polytope Ax <= b.
+%  The Chebyshev center of a polytope is the center of the largest
+%  hypersphere enclosed by the polytope. 
+%  Requires optimization toolbox.
+
+[n,p] = size(A);
+an = sqrt(sum(A.^2,2));
+A1 = zeros(n,p+1);
+A1(:,1:p) = A;
+A1(:,p+1) = an;
+f = zeros(p+1,1);
+f(p+1) = -1;
+
+options = optimset;
+options = optimset(options,'Display', 'off');
+c = linprog(f,A1,b,[],[],[],[],[],options);
+r = c(p+1);
+c = c(1:p);
+end
-- 
cgit v1.2.1