function [A,b,c]=varp(n,w)
A=[ones(1,n);1:n];
b=[1;w];
c=((1:n).^2)';
