clc
    funcprot(0)
    printf("Лабораторная работа №15\n")
    printf(" Задание 1\n")
    function w =f(x,y)
    w=(y/x)+3/((x^2)+x)
    endfunction

    format('v',6)
    x0=3; y0=2;
    x=[3:10];
    z=ode(y0,x0,x,f);
    plot(x,z,'-*'); xgrid()
    mprintf(" x = %2.0f %2.0f %2.0f %2.0f %2.0f %2.0f %2.0f %2.0f\n",x)
    mprintf(" z = %2.0f %2.2f %2.2f %2.1f %2.2f %2.2f %2.2f %2.1f",z)
    a=gca();
    a.x_location = "origin";
    a.y_location = "origin";

    scf
    printf("\n **************************************************\n")
    printf(" Задание 2\n")
    function w = f(x,y)
    w = (y+(x^2)*cos(x))/2
    endfunction

    format('v',5)
    x0=2.5; y0=2.3;
    x=[2.5:0.5:9];
    z=ode(y0,x0,x,f);
    plot(x,z,'-*'); xgrid()
    mprintf(" x = %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f\n",x)
    mprintf(" z = %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f %2.2f",z)

clc
    printf("Лабраторная работа №14\n")
    printf(" Задание 1\n\n")
    x = 1:0.01:5
    y = 2 * sin(x)^2
    Jspl = intsplin(x,y)
    printf(' Jсплайн = %7.5f\n', Jspl)
    Jtr = inttrap(x,y)
    printf(' Jтрап = %7.5f\n', Jtr)

    printf("\n *****************************\n")
    printf(" Задание 2\n")
    Jtr=inttrap([0 1 3 4 6 9 10], [1 3 2 3 5 3 4])
    printf('\n Jтрап = %5.2f\n', Jtr)

    printf("\n *****************************\n")
    printf(" Задание 3\n")
    x = -10:0.01:8
    function y = m(x)
    y = x^2 + sin(x)
    endfunction

    Jintg = intg(-10, 8, m)
    printf('\n Jintg = %7.5f\n', Jintg)

    Jkv = integrate('x^2 + sin(x)', 'x', -10, 8)
    printf(' Jквадратур = %7.5f\n', Jkv)

    y = x^2 + sin(x)
    Jtr = inttrap(x,y)
    printf(' Jтрап = %7.5f\n', Jtr)

    printf("\n *****************************\n")
    printf(" Задание 4\n")
    x = [1 3 4 6 9 10]
    y = log(x)
    Jspl = intsplin(x,y)
    printf('\n Jсплайн = %7.5f\n', Jspl)
    Jtr = inttrap(x,y)
    printf(' Jтрап = %7.5f\n', Jtr)