%% Harj. 1, teht 12.   H1T12R, Polynomien summa
format compact
% Kaytetaan edellisen tehtavan dataa yleisess"a kehittelyss"a
p=[1 0 2 -1]
q=[2 0 0 0 0 3]
np=length(p);
nq=length(q);
m=max(np,nq);
S=zeros(2,m);
S(1,m-np+1:end)=p;
S(2,m-nq+1:end)=q;
summa=sum(S)

%% Sitten muutetaan edellinen ajatus funktioksi polysum
%   Kirjoitetaan tiedostoon polysum.m seuraava (Rivit 18-30):
%%
%{
function [ summa ] = polysum(p,q)
% polysum  laskee argumenttina annettujen 
% polynomivektorien summan
%  Kutsu:
%  S=polysum(p,q);
np=length(p);
nq=length(q);
m=max(np,nq);
S=zeros(2,m);
S(1,m-np+1:end)=p;
S(2,m-nq+1:end)=q;
summa=sum(S);
end
%}

%%
help polysum
p=[1 0 -2]   %x^2 -2
q=[2 2 -1 0] %2*x^3 + 2*x^2 -x
r=polysum(p,q)
