Grzegorz Stolarczyk

Grzegorz Stolarczyk Specjalista ds.
Elektrycznych

Temat: MatDziad

To jest dzięki Alanowi :)


function [w] = splinee (x,f)
lx=length(x);
lf=length(f);
h=zeros(lx-1,1)
h=[x(2)-x(1),x(3)-x(2),x(4)-x(3)]
A=zeros(2,2)
for i=1:2
for j=1:2
if (i==j)
A(i,j)=2*(h(i)+h(i+1))
end
if(i==j+1)
A(i,j)=h(i)
end
if(i==j-1)
A(i,j)=h(i+1)
end
end
end
b=zeros(2,1)
for i=1:2
b(i)=(3/h(i+1))*(f(3)-f(2))-(3/h(i))*(f(2)-f(1))
end
S=zeros(3);
for j=1:3
S(j)=f(j)+b(j)*(x-x(j))

Jakby ktoś jeszcze miał inne programy to wrzucajcie...
Pewnie przyda się Newton i Lagrnaż :)