《一维椭圆方程数值计算》由会员分享,可在线阅读,更多相关《一维椭圆方程数值计算(4页珍藏版)》请在金锄头文库上搜索。
1、 1I 6 9 20140806201 K u du2 dx2 f x x a b u a u b 0 p a b 0 1 uT k C u i 1 2u i u i 1 h2f x i T d Gauss L 0 1 N h 1 N xi ih 0 i N x0 0 xN 1 u00 xi u0 xi x 2 u0 xi x 2 x u xi x u xi u xi u xi x x 2 u xi x 2u xi u xi x x 2 1 x h k u00 xi u xi 1 2u xi u xi 1 h2 1 i N 1 n u xi 1 2u xi u xi 1 u00 xi h2 f
2、 xi h2 z 1 21 1 21 1 2 u x1 u x2 u xN 1 f x1 h2 f x2 h2 f x3 h2 f xN 1 h2 2 u x0 u xN 0 eu x x 1 x f x u00 x 2 gauss 9 function x Gauss s1 A b x Gauss s1 A b A x b n length b A A b for k 1 n 1 A k 1 n k 1 n 1 A k 1 n k 1 n 1 A k 1 n k A k k A k k 1 n 1 A k 1 n k zeros n k 1 end x zeros n 1 x n A n n 1 A n n for k n 1 1 1 x k A k n 1 A k k 1 n x k 1 n A k k end N 500 A zeros N 1 N 1 fori 1 N 2 A i i 2 A i i 1 1 A i 1 i 1 end b ones N 1 1 forj 1 N 1 b j 1 2 N 2 end 2 u1 Gauss s1 A b x1 1 N 1 N 1 1 N plot x1 u1 g hold on u x1 x1 2 plot x1 u r legend u xi u x 1 u x 3 2 u x 100 ddw pd J 4