数值分析3

主函数：
 

clc;
clear all;
close all;

format long;

Aufgabe3_1()
Aufgabe3_2()


function Aufgabe3_1()
    point1 = [0.0; 0.0];
    point2 = [0.577; -0.577];
    f_values = [0.0; 1.0; 3.0; 1.0];%注意间隔开的是；不是，

    val1 = linquadref(point1(1),point1(2));
    val2 = linquadref(point2(1),point2(2));

    f_L1 = val1' * f_values;
    f_L2 = val2' * f_values;

    disp(f_L1);
    disp(f_L2);
end

function Aufgabe3_2()
    point1 = [0.0; 0.0];
    point2 = [0.577; -0.577];
    f_values = [0.0; 1.0; 3.0; 1.0];%注意间隔开的是；不是，

    derviv1 = linquadderivref(point1(1),point1(2));
    derviv2 = linquadderivref(point2(1),point2(2));

    df_L1 = derviv1' * f_values;
    df_L2 = derviv2' * f_values;

    disp(df_L1);
    disp(df_L2);
end

需要用到的函数：二维拉格朗日值和导数
 

function val = linquadref(xi,eta)
    
    N1 = 1/4*(1-xi).*(1-eta);
    N2 = 1/4*(1+xi).*(1-eta);
    N3 = 1/4*(1+xi).*(1+eta);
    N4 = 1/4*(1-xi).*(1+eta);
    val = [N1; N2; N3; N4];
end

function deriv = linquadderivref(xi,eta)
    deriv=[-1/4*(1-eta),-1/4*(1-xi)
                1/4*(1-eta), -1/4*(1+xi)
                1/4*(1+eta),1/4*(1+xi)
                -1/4*(1+eta),1/4*(1-xi)];
end