% Käytetään valmiita työkaluja:

aste = 4;
alapoly = polyfit([1 2 3 4], [0 0 0 0], aste);
oikeapoly = polyfit([1 2 3 4], [20 40 60 80], aste);
ylapoly = polyfit([1 2 3 4], [50 20 10 0], aste);
vasenpoly = polyfit([1 2 3 4], [20 40 60 80], aste);

% Lasketaan ratkaisupolynomien arvot pisteissä X

X = 0:0.5:5;
alaY = polyval(alapoly, X);
oikeaY = polyval(oikeapoly, X);
ylaY = polyval(ylapoly, X);
vasenY = polyval(vasenpoly, X);

pinta = lapsolve1(alaY, oikeaY, ylaY, vasenY);

surf(pinta)

% lapsolve1.m on lainattu malliratkaisuista...
