From baf7ff4b6da3b5670426991e1ed38f95a03d2440 Mon Sep 17 00:00:00 2001 From: Stanislav KROTER Date: Tue, 26 Mar 2013 22:33:01 +0600 Subject: [PATCH] Revision help pages for bvode (en_US, ru_RU). Fixed some typos and changed the example "An eigenvalue problem". This example asked about n-th eigenvalue but did not wait the answer and, as a result, the graphic figure and the sample image were different. Change-Id: I4bbecf02e061a0907df459cc79517f3c81a66af4 --- .../differential_equations/help/en_US/bvode.xml | 32 ++++++++--------- .../differential_equations/help/ru_RU/bvode.xml | 37 ++++++++++---------- 2 files changed, 33 insertions(+), 36 deletions(-) diff --git a/scilab/modules/differential_equations/help/en_US/bvode.xml b/scilab/modules/differential_equations/help/en_US/bvode.xml index ae230cd..d9b5d35 100644 --- a/scilab/modules/differential_equations/help/en_US/bvode.xml +++ b/scilab/modules/differential_equations/help/en_US/bvode.xml @@ -56,7 +56,7 @@ N a scalar with integer value, number of differential equations - (N <= 20). + (N <= 20). @@ -213,7 +213,7 @@ = 0 - causes bvode to generate a uniform initial + causes bvode to generate a uniform initial mesh. @@ -228,12 +228,9 @@ = 1 - if the initial mesh is provided by the user. - it is defined in fspace as follows: the mesh + if the initial mesh is provided by the user it is defined in fspace as follows: the mesh will occupy fspace(1), ..., fspace(n+1). The user needs to supply only - the interior mesh points fspace(j) = x(j), - j = 2, ..., n. - + the interior mesh points fspace(j) = x(j),j = 2, ..., n. @@ -286,8 +283,8 @@ = 3 if a former mesh and approximate solution - coefficients are provided by the user in fspace , and the - new mesh is to be taken twice as coarse; i.e.,every + coefficients are provided by the user in fspace, and the + new mesh is to be taken twice as coarse; i.e. every second point from the former mesh. @@ -360,7 +357,7 @@ an array of dimension ntol=ipar(4). ltol(j) = l specifies that the j-th tolerance in - the tol array controls the error in the l-th + the tol array controls the error in the l-th component of @@ -375,8 +372,8 @@ - - . It is also required that: + . + It is also required that: 1 <= ltol(1) < ltol(2) < ... < ltol(ntol) @@ -532,7 +529,7 @@ x_up, which are to be included in every mesh. The code requires that all side condition points other than x_low and x_up (see - description of zeta ) be included as fixed points in + description of zeta) be included as fixed points in fixpnt. @@ -1431,7 +1428,7 @@ function [dmval,zu]=fsub(x,parameters) The problem 1 describes a - uniformy loaded beam of variable stifness, simply supported at both + uniformly loaded beam of variable stiffness, simply supported at both end. It may be defined as follow : @@ -2146,7 +2143,7 @@ xt=sqrt(2*(gam-1)/gam) fixpnt=[ ];//All boundary conditions are located at x_low and 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 collpnt 10 4 nmax*nsizef nmax*nsizei -1 0 0 0 0 ] @@ -2222,7 +2219,8 @@ N=101; x=linspace(a,b,N)'; // We have s(n)-(n+1/2)*pi -> 0 for n to infinity. -la0=input('n-th eigenvalue: n= ?');la0=(%pi/2+la0*%pi)^2; +la0=evstr(x_dialog('n-th eigenvalue: n= ?','10')); +la0=(%pi/2+la0*%pi)^2; z=bvodeS(x,m,n,a,b,fsub,gsub,zeta,ystart=list(ystart,la0)); // The same call without any display @@ -2233,7 +2231,7 @@ z=bvodeS(x,m,n,a,b,fsub,gsub,zeta,ystart=list(ystart,la0),iprint=-1); clf() plot(x,[z(1,:)' z(2,:)']) xtitle(['Startvalue = '+string(la0);'Eigenvalue = '+string(z(3,1))],'x',' ') -legend(['y(x)';'y''(x)']) +legend(['y(x)';'y''(x)']); ]]> function rhs=fsub(x,z) diff --git a/scilab/modules/differential_equations/help/ru_RU/bvode.xml b/scilab/modules/differential_equations/help/ru_RU/bvode.xml index 1cb3d40..1322544 100644 --- a/scilab/modules/differential_equations/help/ru_RU/bvode.xml +++ b/scilab/modules/differential_equations/help/ru_RU/bvode.xml @@ -124,7 +124,7 @@ max(m(j)) <= collpnt <= 7. - Если ipar(2)=0, то collpnt установлен равным max( max(m(j))+1, 5-max(m(j)) ). + Если ipar(2)=0, то collpnt установлен равным max( max(m(j))+1, 5-max(m(j))). @@ -132,8 +132,7 @@ subint: ipar(3) - Задаёт количество подынтервалов в исходной сетке. Если - ipar(3) = 0, то bvode произвольным образом устанавливается subint = 5. + Задаёт количество подынтервалов в исходной сетке. Если ipar(3) = 0, то bvode произвольным образом устанавливается subint = 5. @@ -360,7 +359,7 @@ массив размерности ntol=ipar(4). ltol(j) = l определяет, что j-тый допуск в массиве - tol управляет ошибкой в l-том элементе + tol управляет ошибкой в l-том элементе @@ -374,8 +373,7 @@ - - . + . Также требуется, чтобы: @@ -407,8 +405,8 @@ - - . Таким образом код пытается удовлетворить + . + Таким образом код пытается удовлетворить @@ -527,7 +525,7 @@ массив размером nfxpnt=ipar(11). Он содержит точки отличные от x_low и x_up, которые нужно включить во все - сетки. Код требует, чтобы все точки дополнительных условий, отличные от + сетки. Код требует, чтобы все точки дополнительных условий, отличные от x_low и x_up (см. описание zeta ) были включены в качестве фиксированных точек в fixpnt. @@ -1795,7 +1793,7 @@ function dg=dgsub(i,z) end endfunction -// - Функция, которая вычисляет начальное предположение, здесь не используется +// - Функция, которая вычисляет начальное предположение, здесь не используемое function [zu,mpar]=guess(x) zu=0; mpar=0; @@ -1812,8 +1810,8 @@ endfunction 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];//установка значений допуска для контроля допуска @@ -2129,16 +2127,16 @@ xt=sqrt(2*(gam-1)/gam) 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) @@ -2204,8 +2202,9 @@ zeta=[a a b]; 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)); // Тот же вызов без вывода на экран @@ -2216,7 +2215,7 @@ z=bvodeS(x,m,n,a,b,fsub,gsub,zeta,ystart=list(ystart,la0),iprint=-1); 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)']); ]]> function rhs=fsub(x,z) -- 1.7.9.5