BSample.m ( File view )

  • By 宝马带斗 2013-11-27
  • View(s):186
  • Download(s):6
  • Point(s): 3
			function f0 = BSample(a,b,n,y,y_1,y_N,x0)
f0 = 0.0;
h = (b-a)/n;
c = zeros(n+3,1);
b = zeros(n+1,1);

for i=0:n-1
    if(a+i*h<=x0) && (a+i*h+h>=x0)
        index = i;
end                              %找到x0所在区间

A = diag(4*ones(n+1,1));
I = eye(n+1,n+1);
AL = [I(2:n+1,:);zeros(1,n+1)];
AU = [zeros(1,n+1);I(1:n,:)];
A = A+AL+AU;                      %形成系数矩阵
for i=2:n
    b(i,1) = 6*y(i);
b(1) = 6*y(1)+2*h*y_1;
b(n+1) = 6*y(n+1)-2*h*y_N;
d = followup(A,b);          %用追赶法求出系数
c(2:n+2) = d;
c(1) = c(2) - 2*h*y_1;      %c(-1)
c(n+3) = c(3)+2*h*y_N;      %c(n+1)

x1 = (a+index*h-
(Please download the complete source code to view)
Expand> <Close

Want complete source code? Download it here

Point(s): 3

0 lines left, continue to read
Sponsored links

File list

Tips: You can preview the content of files by clicking file names^_^
Name Size Date
Atken.m681.00 B24-11-06|22:27
BSample.m988.00 B09-04-08|18:14
DCS.m531.00 B25-11-06|00:51
DH.m1.08 kB09-04-08|18:17
DL.m731.00 B09-04-08|18:15
DTL.m1.66 kB09-04-08|18:16
FCZ.m476.00 B09-04-08|18:15
Gauss.m2.16 kB24-11-06|23:56
Hermite.m661.00 B24-11-06|10:49
Language.m804.00 B24-11-06|22:29
Neville.m616.00 B25-11-06|01:03
Newton.m602.00 B24-11-06|22:39
Newtonback.m725.00 B24-11-06|23:41
Newtonforward.m734.00 B24-11-06|23:25
SecSample.m934.00 B09-04-08|18:12
SubHermite.m844.00 B09-04-08|18:11
ThrSample1.m1.32 kB09-04-08|18:12
ThrSample2.m1.43 kB09-04-08|18:13
ThrSample3.m1.54 kB09-04-08|18:13
Sponsored links

BSample.m (9.25 kB)

Need 3 point
Your Point(s)

Your Point isn't enough.

Get point immediately by PayPal

More(Debit card / Credit card / PayPal Credit / Online Banking)

Submit your source codes. Get more point


Don't have an account? Register now
Need any help?
Mail to:


CodeForge Chinese Version
CodeForge English Version

Where are you going?

^_^"Oops ...

Sorry!This guy is mysterious, its blog hasn't been opened, try another, please!

Warm tip!

CodeForge to FavoriteFavorite by Ctrl+D