http://freakshare.com/files/98izd2lm/Spline.m.htmlla funcion de spline

function [valor,a,b,c,d] = Spline(vy,vx,xinter)
    a=vy;
    c=0;b=0;d=0;
    h=0;A=eye(length(vx));B=zeros([length(vx)],1);
    for i=1:(length(vx)-1)
        h(i)=vx(i+1)-vx(i);
    end
    for i=2:length(h)
        A(i,i)=2*(h(i)+h(i-1));
        A(i,i-1)=h(i-1);
        A(i,i+1)=h(i);
        B(i,1)=((3/h(i))*(a(i+1)-a(i)))-((3/h(i-1))*(a(i)-a(i-1)));
    end
    c=inv(A)*B;
    c=c';
    for(i=1:length(h))
        b(i)=((a(i+1)-a(i))/h(i))-((h(i)/3)*(c(i+1)+2*c(i)));
    end
    for(i=1:length(h))
        d(i)=(c(i+1)-c(i))/(3*h(i));
    end
    for(i=1:length(d)-1)
        if(vx(i)<=xinter&xinter<=vx(i+1))
            valor=a(i)+b(i)*(xinter-vx(i))+c(i)*(xinter-vx(i))^2+d(i)*(xinter-vx(i))^3
        end
    end
    disp(valor)
end

Aca un ejemplo de la vida sobre sistemas mal condicionados http://folk.uib.no/ssu029/Pdf_file/Deuflhard74.pdf si no lo pueden bajar mañana se los paso :)

SISTEMAS MAL CONDICIONADOS

SISTEMAS MAL CONDICIONADOS http://www.cibernetia.com/tesis_es/CIENCIAS_DE_LA_TIERRA_Y_DEL_ESPACIO/GEODESIA/GEODESIA_DE_SATELITES/1 Nos habla acerca de la aplicacion de sistemas de ecuaciones para la observacion de satelites, los sistemas mal condicionados con frecuencia en múltiples campos relacionados con las Ciencias de la Tierra. Tambien pueden escuchar el resumen de la tesis en la pagina

Hola chicos este es el link de descarga del libro:   Basta de Historias-Andres Oppenheimer
http://www.mediafire.com/?eafkt6a5yc071gv

Funcion solo para dos variables
function [val,ea]=error2(f,x1,ex,y1,ey)
syms x y
dy=diff(f,y);
dx=diff(f,x);
dx=inline(dx);
dy=inline(dy)
ea=ex*dx(x1,y1)+ey*dy(x1)
f=inline(f);
val=f(x1,y1)
ea/val
end
ejemplo en el matlab
error2('x^2+4*x*y',2,0.2,3,0.1)
pruebenlo

Esta es la funcion que nos sirve para la propagacion de errores para dos varialbles:

function propagacion2var(f,x0,y0,ex,ey)
    syms x
    syms y;
    y1=subs(f,{x,y},{x0,y0});
    f1=diff(f,x);
    f2=diff(f,y);
    valF1=subs(f1,{x,y},{x0,y0});
    valF2=subs(f2,{x,y},{x0,y0});
    ea=abs(valF1)*ex+abs(valF2)*ey;
    er=ea/y1;
    disp('valor de y:');
    disp(y1);
    disp('error:');
    disp(ea)
    disp('er:')
    disp(er);
end


y ejecutenlo de la siguiente forma:
>> propagacion2var('x^2+4*x*y',5,3,0.2,0.1)

y el resultado sera:
valor de y:
    85

error:
32/5

er:
32/425

el resultado nos devolvera en fracciones pero si lo dividen les dara el mismo resultado