%% 7. {(-2,-15),(-1,-8),(0,-3)}À» º¸°£ÇÏ´Â 2Â÷ ´ÙÇ×½ÄÀ» ¸¸µå½Ã¿À. ±×¸®°í p2(x)ÀÇ ±Ù Áß¿¡¼ x = 0ÀÇ ¿ìÃøÀ¸·Î ÇâÇÑ Ã¹ ¹ø° ±ÙÀ» ãÀ¸½Ã¿À. ´ÙÇ×½Ä p2(x)°¡ ÃÖ´ë°ªÀ» °®´Â°¡? % x0=-2; x1=-1; x2=0; % y0=-15; y1=-8; y2=-3;
% L2,0(x) = (x-x1)*(x-x2)/(x0-x1)*(x0-x2); % L2,1(x) = (x-x0)*(x-x2)/(x1-x0)*(x1-x2); % L2,2(x) = (x-x0)*(x-x1)/(x2-x0)*(x2-x1);
% p2(x) = L2,0(x)*(-15) + L2,1(x)*(-8) + L2,2(x)*(-3); % p2(x) = (-20/3)*x^2 + (37/3)*x - 3;
f=@(x) (-20/3)*x.^2 + (37/3)*x - 3; h = 1; xj = -2:h:0; yj = f(xj); xv = -2:0.01:0; fv=f(xv); yv = interp_dd(xj,yj,xv); plot (xv,fv) title('f(x) = (-20/3)*x^2 + (37/3)*x - 3') hold on plot(xj,yj,'ro')
%%%%%%%%%%%½ÇÇà°á°ú%%%%%%%%%%%%%%%%% ------ Divided Differences of f -------- i x_i f(x_i) D_i ---------------------------------------- 0 -2.00 -54.3333 -5.4333e+001 1 -1.00 -22.0000 3.2333e+001 2 0.00 -3.0000 -6.6667e+000 ----------------------------------------
>> solve('(-20/3)*x^2 + (37/3)*x - 3=0') ans = 37/40-1/40*649^(1/2) 37/40+1/40*649^(1/2)
ÃÖ´ë°ª = 37/40 = 0.9250
|
|