Answer:
clc
clear all
close all
format long
f=@(n) sum((1:n).^4);
n=1:6;
A=[n(:).^5 n(:).^4 n(:).^3 n(:).^2 n(:) ones(numel(n),1)];
b=[f(1);f(2);f(3);f(4);f(5);f(6)];
C=A\b;
for i=5:-1:0
fprintf('c%d = %f\\',i,C(5-i+1))
end
8.1m questions
10.7m answers