% Jatkuvan tehtävän ominaisfkt: 
% v(k,x)=sin(beta(k)*x), missä beta(k)=k*pi
% Tapa 2: ulkotulo (vrt. diskreetti.

clear all; close all
n=4;  % Muuttele tätä!

x=linspace(0,1);
v=zeros(length(x),n);
k=1:n;
v=sin(pi*k'*x);
plot(x',v')


%%%%%%%%
%Diskreetin teht om. fkt.

%%%%%%%%%%%%%%%%%%%%%%
hold on

X=linspace(0,1,n+2);
%X=X(2:end-1);   % (Annetaan nyt päätepisteidenkin olla)
K=1:n;
V=sin(pi*K'*X);
plot(X,V,'o')
title(['n= ',num2str(n)]) ;grid on; shg
