求助，逼近积分

splelvis

u''(x)=f(x),
: _1 ?2 w. R! V- t" C) R边界boundary conditions u(a)=0,u(b)=0.
5 F* g4 I+ [* P/ ?1 K(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

1 f$ f( t( e" J5 |7 h1 C1 \maltab code:
0 N1 h2 k' v2 ]/ q
clear all
' ?+ J0 @3 j- n3 t6 D  P" fa=0;
b=1;
; C  [; H0 f2 g/ @3 A( |  Q  Z9 q; Ln=10;
h=(b-a)/(n+1);
# m5 T) {) ]' X4 L) sx_with_boundary=linspace(a,b,n+2)';
x=x_with_boundary(2:n+1);
( M/ V# G2 P+ F1 pA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
+ ^% v* n" H4 Q. ~- ]rhs=f4(x);
; d, F, s8 e; X3 k( U9 jsol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
$ O! @+ z1 r0 Y  h' M8 @plot(x_with_boundary,sol_with_boundary_conds);


( ~7 R! N- ]2 {5 r另开一文档，f4,
* j0 W# x9 B9 a$ v$ l  Hfunction y=f4(x)
8 h3 V; b2 {- v, t1 j( _y=ones(length(x),1);

出来后是抛物线。


! K* e( g" y+ t4 O# e# m现在问题是，
- A3 Z* A8 i' s- q4 ^4 m( iu''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]
- y4 W$ L  w( N, h
5 _) [6 }; F- g: m1 Y请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
) h3 G) g! Z* y  T7 r还是有什么资料不明确？