求助，逼近积分

splelvis

u''(x)=f(x),
$ }* s* @& W" ?+ G边界boundary conditions u(a)=0,u(b)=0.
: P5 I/ l! a  R  M- C(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

! o3 u4 Y6 I# lmaltab code:
! b( E1 {7 @- P1 l; {* s6 X4 z
( A" k& D( x; @% g9 P; q8 Bclear all
4 O  J6 S/ B0 K+ ja=0;
b=1;
n=10;
h=(b-a)/(n+1);
) g' Z8 k$ H6 J9 jx_with_boundary=linspace(a,b,n+2)';
# L7 t5 p" s0 q- Jx=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);
) I+ e9 Z0 H0 e* Ysol=A\rhs;
3 g8 F8 G: x! S% t( g1 ysol_with_boundary_conds=[0;sol;0];
) _9 k! K9 z4 j# D, E; Kplot(x_with_boundary,sol_with_boundary_conds);


0 G! P4 u7 |  p! @& y3 v另开一文档，f4,
function y=f4(x)
: E! T* W- q" R; C9 C5 q; ry=ones(length(x),1);

2 R* N# V8 u7 X出来后是抛物线。


5 p8 H7 X$ G* P7 Z  c现在问题是，
u''(x)=sin(2pix);
0 z' W% J$ y3 qu(-1)=0;
5 j; j* j0 {" Uu(1)=0;
范围[-1;1]
+ U0 p2 r. X* w, L6 o3 j' u. G
请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
8 a3 S9 M! |1 U5 g: E# Z0 J还是有什么资料不明确？