:ڶVolterraֺ

ʽ:
[z,er,n]=IntegralEquVolterra2(x,K,F,a,n0,n1)
[z,er,n]=IntegralEquVolterra2(x,K,F,a,n0)
[z,er,n]=IntegralEquVolterra2(x,K,F,a)
[z,er,n]=IntegralEquVolterra2(x,K,F)

x:ֵָ,˱Ǿ
K:ű̵ַں˺,ıΪx,y,yʾֱ
F:ű洢ұֻxƵĺ
a:,ĬΪ0
n0:ø˹ֵĳʼ,ĬΪ2
n1:ø˹ֵ󸴻,ĬΪ80

z:صxֵӦļ
er:ص
n:ؼ㵽ֵʱĸ˹

ע:Ĳxеֵ붼ڵa߶Сڵa.ֵaߵֵ,ִα.

ԭ:

1µ

  z(x) = Int[K(x,y)*z(y) , y , a , x] + F(x)

  Int[K(x,y)*z(y) , y , a , b]ʾK(x,y)*z(y)ıy[ax]ж.

2Day3ɭʽ,Ȼ̲3ɭʽ4ɭʽιʽϵķʽк,вǰμΪжϼֹ.

ο:ִӦѧֲ-ֵ[M].廪ѧ,,2007:784-788 

:

kf = "x*y/5";
ff = "x";
x=[0,0.1,0.2,0.3,0.6,0.9,1.2,1.5,1.8,1.9,2.0];

[y,er,n]=IntegralEquVolterra2(x,kf,ff)//سõ½,ʵ y=x*exp(x^3/15)
y =
[ 0.00000000000000    0.10000666688889    0.20010669512925    0.30054048629174    0.60870250767669    0.94482031194692    1.34651740164395    1.87848407430932    2.65537706395036    3.00152518019025    3.40920973236942 ]
er =
[ 1.0159937300E-16 ]
n =
[ 1600.00000000000 ]