function [y0] = solve(w,x0)
n = max(size(w));
syms y x u0 u1 tmp;
y = w(1)*x^(0);
for i = 2:n
    y = y + w(i)*x^(i-1);
end
y = simplify(y);
ezplot(y);
y
x = x0;
y0 = subs(y);
y0 = double(y0);
end