f=@(x) exp(x)./(1+x.^2);

>>I=integral(f,0,1)

I =

    1.2707

>> t=@(x) sin(x);
>> [A,b]=sistem(1,t,5)

A =

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000
         0    0.2000    0.4000    0.6000    0.8000    1.0000
         0    0.0400    0.1600    0.3600    0.6400    1.0000
         0    0.0080    0.0640    0.2160    0.5120    1.0000
         0    0.0016    0.0256    0.1296    0.4096    1.0000
         0    0.0003    0.0102    0.0778    0.3277    1.0000


b =

    0.4597
    0.3012
    0.2232
    0.1771
    0.1467
    0.1251


>> koeficijenti(1,t,5)
koeficijenti su:
    0.0051    0.0270    0.1151    0.0531    0.2077    0.0517