求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
" w8 ~4 I" }; d  f: {
maltab code:
# W) o2 `- ]( B, i- j# \
4 A1 {3 Y& `$ @clear all
' @$ |1 d' S/ }0 ]& a' Ga=0;
b=1;
  W9 c# z4 l2 p" r) t) x; Yn=10;
- g5 g/ P- M1 M  D8 j) Sh=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
* ~' O1 I+ e, |. w; s9 ox=x_with_boundary(2:n+1);
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
rhs=f4(x);
9 S; r. I; {. p; L' A! gsol=A\rhs;
) e' ?8 ~- q# N6 Q  O( |! Hsol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
; D, p, ^5 L5 ]" E9 U
8 Q8 t" g8 y( f) ~
4 B2 B6 \! i8 s& [7 @: r$ w另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);
2 h1 Z6 z1 w% r) ~3 v+ {
) g2 V  c$ s0 \1 t7 ~7 \9 p出来后是抛物线。

6 T  z! Y8 B7 o! s) R$ J
8 O1 W" z+ ^2 @& w1 C9 F# q5 J现在问题是，
9 `6 X; B0 w' Uu''(x)=sin(2pix);
u(-1)=0;
) Z; f: B7 l0 N8 `. @& Ku(1)=0;
! {' z! @% x, b) f范围[-1;1]

! N9 U) }7 m, H) B请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
+ d0 b* J7 {! U- @2 Z% [还是有什么资料不明确？