% Aufgabe 10.3 Darstellung der Daubechies Skalierungsfunktion D4
function y=aufgabe103(anz)
% clf reset
% clc
% clear all;
% close all;
% echo off
% anz Anzahl der berechneten Punkte ist 3*2^J +1

c0=(1+sqrt(3))/4; c1=(3+sqrt(3))/4; c2=(3-sqrt(3))/4; c3=(1-sqrt(3))/4;
v0= [0; 2+sqrt(3); -1];
v1=[2+sqrt(3); -1;0];
A0=[c0 0 0 ; c2 c1 c0 ; 0 c3 c2];
A1=[c1 c0 0; c3 c2 c1 ; 0 0 c3];

anz
N=2^anz;
phi=zeros(3,N);
phi(:,1)
phi(:,1)=v0;
phi(:,N/2+1)= A0*v1; 
for j=2:anz 
  for k=1:2^(j-2)
   phi(:,N*(2*k-1)/2^j+1)=A0*phi(:,N*(2*k-1)*2/2^j+1);  phi(:,N*(2*k-1)/2^j+N/2+1)=A1*phi(:,2*N*(2*k-1)/2^j+1); 
  end;
end;

z=0:1/N:3;
z1=[phi(1,:), phi(2,:), phi(3,:), 0];
figure
   plot(z,z1);
    set(gca,'FontSize',15);

