<literal>max(m(j)) <= collpnt <= 7</literal>.
</para>
<para>
- Если <literal>ipar(2)=0</literal>, то <literal>collpnt</literal> установлен равным <literal>max( max(m(j))+1, 5-max(m(j)) )</literal>.
+ Если <literal>ipar(2)=0</literal>, то <literal>collpnt</literal> установлен равным <literal>max( max(m(j))+1, 5-max(m(j)))</literal>.
</para>
</listitem>
</varlistentry>
<term>subint: ipar(3)</term>
<listitem>
<para>
- Задаёт количество подынтервалов в исходной сетке. Если
- <literal>ipar(3) = 0</literal>, то <literal>bvode</literal> произвольным образом устанавливается <literal>subint = 5</literal>.
+ Задаёт количество подынтервалов в исходной сетке. Если <literal>ipar(3) = 0</literal>, то <literal>bvode</literal> произвольным образом устанавливается <literal>subint = 5</literal>.
</para>
</listitem>
</varlistentry>
<para>
массив размерности <literal>ntol=ipar(4)</literal>.
<literal>ltol(j) = l</literal> определяет, что <literal>j</literal>-тый допуск в массиве
- <literal> tol</literal> управляет ошибкой в <literal>l</literal>-том элементе
+ <literal>tol</literal> управляет ошибкой в <literal>l</literal>-том элементе
<inlinemediaobject>
<imageobject>
<imagedata>
</mml:math>
</imagedata>
</imageobject>
- </inlinemediaobject>
- .
+ </inlinemediaobject>.
Также требуется, чтобы:
</para>
<para>
</mml:math>
</imagedata>
</imageobject>
- </inlinemediaobject>
- . Таким образом код пытается удовлетворить
+ </inlinemediaobject>.
+ Таким образом код пытается удовлетворить
<inlinemediaobject>
<imageobject>
<imagedata>
<para>
массив размером <literal>nfxpnt=ipar(11)</literal>. Он содержит точки отличные от
<literal>x_low</literal> и <literal>x_up</literal>, которые нужно включить во все
- сетки. Код требует, чтобы все точки дополнительных условий, отличные от
+ сетки. Код требует, чтобы все точки дополнительных условий, отличные от
<literal>x_low</literal> и <literal>x_up</literal> (см. описание <literal>zeta</literal> )
были включены в качестве фиксированных точек в <literal>fixpnt</literal>.
</para>
end
endfunction
-// - Функция, которая вычисляет начальное предположение, здесь не используется
+// - Функция, которая вычисляет начальное предположение, здесь не используемое
function [zu,mpar]=guess(x)
zu=0;
mpar=0;
fixpnt=[ ];//Все граничные условия находятся в x_low и x_up
-// nonlin collpnt n ntol ndimf ndimi iprint iread iguess rstart nfxpnt
-+ipar=[1 collpnt 10 4 nmax*nsizef nmax*nsizei -1 0 0 0 0 ]
+// nonlin collpnt n ntol ndimf ndimi iprint iread iguess rstart nfxpnt
+ipar=[0 0 1 2 2000 200 1 0 0 0 0 ]
ltol=[1,3];//установка контроля допуска на zu(1) и zu(3)
tol=[1.e-11,1.e-11];//установка значений допуска для контроля допуска
fixpnt=[ ];// все граничные условия размещены в x_low и x_up
collpnt=4;
nsizef=4+3*M+(5+collpnt*N)*(collpnt*N+M)+(2*M-2)*2*M ;
-nsizei=3 + collpnt*N+M;;
+nsizei=3 + collpnt*N+M;
nmax=200;
// nonlin collpnt n ntol ndimf ndimi iprint iread iguess rstart nfxpnt
-ipar=[1 k 10 4 nmax*nsizef nmax*nsizei -1 0 0 0 0 ]
+ipar=[1 collpnt 10 4 nmax*nsizef nmax*nsizei -1 0 0 0 0 ]
ltol=1:4;//установка контроля допусков на zu(1), zu(2), zu(3) и zu(4)
-tol=[1.e-5,1.e-5,1.e-5,1.e-5];//установка значений для контроля допусков
+tol=[1.e-5,1.e-5,1.e-5,1.e-5];//установка значений допусков для этих четырёх значений контроля
xpoints=x_low:0.01:x_up;
-// - Функция, которая вычисляет начальное предположение, здесь не используется
+// - Функция, которая вычисляет начальное предположение, здесь не используемое
function [zu,dmval]=guess2(x,gam),
cons=gam*x*(1-x^2/2)
dcons=gam*(1-3*x^2/2)
N=101;
x=linspace(a,b,N)';
-//Имеем s(n)-(n+1/2)*pi -> 0 для n, стремящимся к бесконечности.
-la0=input('n-ное собственное значение: n= ?');la0=(%pi/2+la0*%pi)^2;
+//Имеем s(n)-(n+1/2)*pi -> 0 для n, стремящемся к бесконечности.
+la0=evstr(x_dialog('n-ное собственное значение: n= ?','10'));
+la0=(%pi/2+la0*%pi)^2;
z=bvodeS(x,m,n,a,b,fsub,gsub,zeta,ystart=list(ystart,la0));
// Тот же вызов без вывода на экран
clf()
plot(x,[z(1,:)' z(2,:)'])
xtitle(['Начальное значение = '+string(la0);'Собственное значение = '+string(z(3,1))],'x',' ')
-legend(['y(x)';'y''(x)'])
+legend(['y(x)';'y''(x)']);
]]></programlisting>
<scilab:image localized="true">
function rhs=fsub(x,z)