ritz_method.sci
clc clear // d2y/dx2 + 1000x^2 = 0 // y(0) = 0, y(1) = 0, function r = res(a,x) // n1 = (x-xa)*(xb-xa)*x; // yp = a1*n1 ddn1 = 2.0; r = a*ddn1 +1000*x.^2; endfunction function r = integr(a) x = linspace(0,1.0); disp(x) r = inttrap(x, res(a,x)); endfunction a = fsolve(0, integr); disp(a) x = linspace(0,1,20); yp = a(1)*x.*(x-1) plot(x,yp,'b')