function [] = jus ()
t = -2:0.01:2;
y = xpi(t, -1/2, 3);
y2 = xpi(t, 0, 1/2);

y3 = y - y2;

w1=xpi(t,-1,2) - xpi(t,1,2);
w2=xpi(t,-3/2,1)-xpi(t,0,2)+xpi(t,3/2,1);
w3=xpi(t,-3/2,1)-xpi(t,-1/2,1)+xpi(t,1/2,1)-xpi(t,3/2,1);

w1_calka = calka(w1.*w1,0.001),
w2_calka= calka(w2.*w2, 0.001),
w3_calka= calka(w3.*w3, 0.001),

a1 = calka(y3.*w1,0.001) ./ w1_calka,
a2 = calka(y3.*w2,0.001) ./ w2_calka,
a3 = calka(y3.*w3,0.001) ./ w3_calka,
a = [a1,a2,a3],

xa1=calka(a1.*xpi(t,-1,2),0.001) + calka(a2.*xpi(t,-3/2,1),0.001) + calka(a3.* xpi(t, -3/2 , 1), 0.001);
xa2 = calka(a1.* xpi(t, -1, 2), 0.001) + calka(a2.* xpi(t,-1/2, 1),0.001).*(-1) + calka(a3.*xpi(t, -1/2, 1),0.001).*(-1);
xa3 = calka(a1.*xpi(t, 1, 2),0.001).*(-1) + calka(a2.*xpi(t, 1/2, 1),0.001).*(-1) + calka(a3.*xpi(t, 1/2,1),0.001);
xa4 = calka(a1.*xpi(t, 1,2),0.001).*(-1) + calka(a2.* xpi(t, 3/2, 1),0.001) + calka(a3.*xpi(t,3/2,1),0.001).*(-1);

xa = [xa1,xa2,xa3,xa4],


grid on;
figure(1)
 plot(t,y3,'r')
 figure(2)
 plot(t,w1)
 figure(3)
 plot(t,w2)
 figure(4)
 plot(t,w3)
 figure(5)

 xapp=  xa1.* xpi(t, -3/2,1) + xa2.*xpi(t, -1/2,1) + xa3.*xpi(t, 1/2 ,1) + xa4.*xpi(t,3/2,1);

 plot(t, xapp)
 grid on;