求助，逼近积分

splelvis

u''(x)=f(x),
* z  t4 q) O- ~9 i边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

  Q' x7 E( G* z' ^$ K) nmaltab code:
; H" t2 ~7 n9 w7 B
. i% E' w! p$ |; U6 D" wclear all
a=0;
b=1;
2 a- H9 v. N2 S9 tn=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
( j4 O+ T9 I# i# Ux=x_with_boundary(2:n+1);
& b" J8 [" U" I. Z* _A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
rhs=f4(x);
3 E5 r) c; l: z8 bsol=A\rhs;
' [0 ?5 O$ Q1 }* `: R3 T* Isol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);

) J8 s: s$ O+ n$ N3 u! x
: W) J6 {& o! n* K) L另开一文档，f4,
function y=f4(x)
; _2 C2 e  a/ K% I& o" p: S& ^y=ones(length(x),1);
# W# k+ y% F/ v( i0 h" e; |) z* Y8 y
出来后是抛物线。

8 [' R" O9 k  K; Q" B
现在问题是，
u''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]
; x4 V# t" @4 d* w# v; e8 m8 O
+ M% s, d2 e# K9 T4 H请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
; k' ?8 I7 J4 A7 f# ]+ z2 {& L还是有什么资料不明确？