function [fungsi1,fungsi2,fungsi3]=Soal1b(fs1,a,fs2,b,fs3,c,n)
delta1=(a(2)-a(1))/n;
delta2=(b(2)-b(1))/n;
delta3=(c(2)-c(1))/n;
for i=1:1:n
    x1(i)=a(1)+(i-1)*delta1;
    x2(i)=b(1)+(i-1)*delta2;
    x3(i)=c(1)+(i-1)*delta3;
    nx1=a(1)+(i-1)*delta1;
    nx2=b(1)+(i-1)*delta2;
    nx3=c(1)+(i-1)*delta3;
    fungsi1(i)=fs1(nx1);
    fungsi2(i)=fs2(nx2);
    fungsi3(i)=fs3(nx3);
end;
plot(x1,fungsi1,x2,fungsi2,x3,fungsi3);
return;