MWR

//podobszary

function alfa2=a2(x)

    alfa2=x^2-x+2

endfunction

function prawa=p(x)

    prawa=-x

endfunction

x0=0

x1=1

A=integrate('a2(x)','x',x0,x1)

B=integrate('p(x)','x',x0,x1)

roz=linsolve(A,B)

x=0:0.1:x1

y=1-x+roz*(x^2-x)

plot2d(x,y,4)

/////M.naj.kwadr

function alfa2=a2(x)

    alfa2=x^2-x+2

endfunction

function prawa=p(x)

    prawa=-x

endfunction

function waga=w(x) //pochodna po a2 z u

    waga=x^2-x+2

endfunction

x0=0

x1=1

A=integrate('a2(x)*w(x)','x',x0,x1)

B=integrate('p(x)*w(x)','x',x0,x1)

roz=linsolve(A,B)

x=0:0.1:x1

y=1-x+roz*(x^2-x)

plot2d(x,y,5)

//gallerkin pewny

//Residuum R(x)=a1(x^2-x+2)+a2(x^4+-2x^3+13x^2-12x+2)

c

// gallerkin ?????

function fi3=f3(x)

    fi3=x^2

endfunction

function alfa2=a2(x)

    alfa2=x^2-x+2

endfunction

function prawa=p(x)

    prawa=-x

endfunction

x0=0

x1=1

A=integrate('a2(x)*f3(x)','x',x0,x1)

B=integrate('p(x)*f3(x)','x',x0,x1)

roz=linsolve(A,B)

x=0:0.1:x1

y=1-x+roz*(x^2-x)

plot2d(x,y,6)