Versão do Scilab: 5.5.2 clear clc function [orderValues] = getOrder(lastOrder, vectorX, orderNumber) for i = 1 : size(lastOrder, 2) - 1 // size 1 = quantas linhas, size 2 = quantas colunas orderValues(1,i) = (lastOrder(1,i+1) - lastOrder(1,i)) / (vectorX(1,orderNumber+i-1) - vectorX(1,i)) end endfunction function [f, polyF] = newtonInterpol(vectorY, vectorX) lastOrder = vectorY polyF = poly([vectorY(1,1)],"x", "coeff") for i = 2 : size(vectorX, 2) // size 1 = quantas linhas, size 2 = quantas colunas lastOrder = getOrder(lastOrder, vectorX, i) polyA = poly([-vectorX(1,1) 1],"x","coeff") for j = 2 : i-1 polyA = polyA * poly([-vectorX(1,j) 1],"x","coeff") end polyA = polyA * lastOrder(1,1) ...
Comentários
Postar um comentário
Obrigado pelo contato, retornarei o mais breve possível!