From: Stanislav KROTER Date: Sat, 30 Mar 2013 14:22:20 +0000 (+0600) Subject: Revision of help pages for dae function (en_US, ru_RU). X-Git-Tag: 5.5.0-beta-1~867 X-Git-Url: http://gitweb.scilab.org/?p=scilab.git;a=commitdiff_plain;h=c26850e5bdcfc4c78af2225acd1c1d18ce63881c Revision of help pages for dae function (en_US, ru_RU). Inserted a checking for a C compiler. Change-Id: I98493151a32cf485aec5c11d6ba7cd4e74b48c49 --- diff --git a/scilab/modules/differential_equations/help/en_US/dae.xml b/scilab/modules/differential_equations/help/en_US/dae.xml index c6107af..d74f714 100644 --- a/scilab/modules/differential_equations/help/en_US/dae.xml +++ b/scilab/modules/differential_equations/help/en_US/dae.xml @@ -453,13 +453,14 @@ x(t0)=x0 and xdot(t0)=xdot0 Examples %t + [btn] = messagebox(["You need a C compiler for this example.";"Execution of this example is canceled."], "Software problem", 'info'); + return +end + //-1- create the C codes in TMPDIR - Vanderpol equation, implicit form code=['#include ' 'void res22(double *t,double *y,double *yd,double *res,int *ires,double *rpar,int *ipar)' @@ -493,11 +501,14 @@ code=['#include ' ' ' 'void gr22(int *neq, double *t, double *y, int *ng, double *groot, double *rpar, int *ipar)' '{ groot[0] = y[0];}'] +previous_dir = pwd(); cd TMPDIR; mputl(code, 't22.c') + //-2- compile and load them -ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c',TMPDIR+'/Makefile',TMPDIR+'/t22loader.sce'); +ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c','Makefile','t22loader.sce'); exec('t22loader.sce') + //-3- run rtol=[1.d-6;1.d-6]; atol=[1.d-6;1.d-4]; @@ -505,6 +516,7 @@ t0=0;y0=[2;0]; y0d=[0;-2]; t=[20:20:200]; ng=1; + //simple simulation t=0:0.003:300; yy=dae([y0,y0d],t0,t,atol,rtol,'res22','jac22'); @@ -522,6 +534,7 @@ y0d1=yy(3:4,qq); [yy,nn,hotd]=dae("root",[y01,y0d1],t01,300,atol,rtol,'res22','jac22',ng,'gr22',hotd); plot(yy(1,1),yy(2,1),'r+') xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))) +cd(previous_dir); ]]> ' @@ -537,9 +550,10 @@ code=['#include ' ' ' 'void gr22(int *neq, double *t, double *y, int *ng, double *groot, double *rpar, int *ipar)' '{ groot[0] = y[0];}'] +previous_dir = pwd() cd TMPDIR; mputl(code, 't22.c') -ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c',TMPDIR+'/Makefile',TMPDIR+'/t22loader.sce'); +ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c','Makefile','t22loader.sce'); exec('t22loader.sce') rtol=[1.d-6;1.d-6]; atol=[1.d-6;1.d-4]; @@ -559,7 +573,8 @@ y01=yy(2:3,qq); y0d1=yy(3:4,qq); [yy,nn,hotd]=dae("root",[y01,y0d1],t01,300,atol,rtol,'res22','jac22',ng,'gr22',hotd); plot(yy(1,1),yy(2,1),'r+') -xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))) +xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))); +cd(previous_dir); ]]> diff --git a/scilab/modules/differential_equations/help/ru_RU/dae.xml b/scilab/modules/differential_equations/help/ru_RU/dae.xml index 766e4d4..cd96136 100644 --- a/scilab/modules/differential_equations/help/ru_RU/dae.xml +++ b/scilab/modules/differential_equations/help/ru_RU/dae.xml @@ -35,9 +35,9 @@ Ð²ÐµÐºÑÐ¾Ñ-ÑÑÐ¾Ð»Ð±ÐµÑ. ÐÐ½ Ð¼Ð¾Ð¶ÐµÑ Ð±ÑÑÑ ÑÐ°Ð²ÐµÐ½ x0 Ð¸Ð»Ð¸ - [x0;xdot0], Ð³Ð´Ðµ x0 -- ÑÑÐ¾ + [x0;xdot0], Ð³Ð´Ðµ x0 - ÑÑÐ¾ Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ Ð² Ð½Ð°ÑÐ°Ð»ÑÐ½ÑÐ¹ Ð¼Ð¾Ð¼ÐµÐ½Ñ Ð²ÑÐµÐ¼ÐµÐ½Ð¸ t0, - Ð° xdot0 -- Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ Ð¿ÑÐ¾Ð¸Ð·Ð²Ð¾Ð´Ð½Ð¾Ð¹ Ð¸ÑÑÐ¾Ð´Ð½Ð¾Ð³Ð¾ + Ð° xdot0 - Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ Ð¿ÑÐ¾Ð¸Ð·Ð²Ð¾Ð´Ð½Ð¾Ð¹ Ð¸ÑÑÐ¾Ð´Ð½Ð¾Ð³Ð¾ ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ Ð¸Ð»Ð¸ ÐµÐ³Ð¾ Ð¾ÑÐµÐ½ÐºÐ° (ÑÐ¼. Ð½Ð¸Ð¶Ðµ). @@ -178,35 +178,35 @@ - t -- ÑÐµÐºÑÑÐµÐµ Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ + t - ÑÐµÐºÑÑÐµÐµ Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ Ð²ÑÐµÐ¼ÐµÐ½Ð¸; - x -- Ð¼Ð°ÑÑÐ¸Ð² ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ; + x - Ð¼Ð°ÑÑÐ¸Ð² ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ; - xdot -- Ð¼Ð°ÑÑÐ¸Ð² Ð¿ÑÐ¾Ð¸Ð·Ð²Ð¾Ð´Ð½ÑÑ + xdot - Ð¼Ð°ÑÑÐ¸Ð² Ð¿ÑÐ¾Ð¸Ð·Ð²Ð¾Ð´Ð½ÑÑ ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ; - res -- Ð¼Ð°ÑÑÐ¸Ð² Ð¾ÑÑÐ°ÑÐºÐ¾Ð²; + res - Ð¼Ð°ÑÑÐ¸Ð² Ð¾ÑÑÐ°ÑÐºÐ¾Ð²; - ires -- Ð¸Ð½Ð´Ð¸ÐºÐ°ÑÐ¾Ñ + ires - Ð¸Ð½Ð´Ð¸ÐºÐ°ÑÐ¾Ñ Ð²ÑÐ¿Ð¾Ð»Ð½ÐµÐ½Ð¸Ñ; - rpar -- Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ»Ð¾ÑÐ¸ÑÐ»ÐµÐ½Ð½ÑÑ + rpar - Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ»Ð¾ÑÐ¸ÑÐ»ÐµÐ½Ð½ÑÑ Ð·Ð½Ð°ÑÐµÐ½Ð¸Ð¹ Ð¿Ð°ÑÐ°Ð¼ÐµÑÑÐ¾Ð² Ñ Ð¿Ð»Ð°Ð²Ð°ÑÑÐµÐ¹ Ð·Ð°Ð¿ÑÑÐ¾Ð¹, ÐºÐ¾ÑÐ¾ÑÑÐµ Ð½ÑÐ¶Ð½Ñ, Ð½Ð¾ Ð½Ðµ Ð¼Ð¾Ð³ÑÑ Ð±ÑÑÑ ÑÑÑÐ°Ð½Ð¾Ð²Ð»ÐµÐ½Ñ Ñ Ð¿Ð¾Ð¼Ð¾ÑÑÑ ÑÑÐ½ÐºÑÐ¸Ð¸ dae. @@ -214,7 +214,7 @@ - ipar -- Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ»Ð¾ÑÐ¸ÑÐ»ÐµÐ½Ð½ÑÑ + ipar - Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ»Ð¾ÑÐ¸ÑÐ»ÐµÐ½Ð½ÑÑ Ð·Ð½Ð°ÑÐµÐ½Ð¸Ð¹ Ð¿Ð°ÑÐ°Ð¼ÐµÑÑÐ¾Ð² Ñ Ð¿Ð»Ð°Ð²Ð°ÑÑÐµÐ¹ Ð·Ð°Ð¿ÑÑÐ¾Ð¹, ÐºÐ¾ÑÐ¾ÑÑÐµ Ð½ÑÐ¶Ð½Ñ, Ð½Ð¾ Ð½Ðµ Ð¼Ð¾Ð³ÑÑ Ð±ÑÑÑ ÑÑÑÐ°Ð½Ð¾Ð²Ð»ÐµÐ½Ñ Ñ Ð¿Ð¾Ð¼Ð¾ÑÑÑ ÑÑÐ½ÐºÑÐ¸Ð¸ dae. @@ -242,7 +242,7 @@ r=jac(t,x,xdot,cj), Ð° jac Ð´Ð¾Ð»Ð¶Ð½Ð° Ð²Ð¾Ð·Ð²ÑÐ°ÑÐ°ÑÑ r=dg(t,x,xdot)/dy+cj*dg(t,x,xdot)/dxdot, - Ð³Ð´Ðµ cj -- Ð²ÐµÑÐµÑÑÐ²ÐµÐ½Ð½ÑÐ¹ ÑÐºÐ°Ð»ÑÑ. + Ð³Ð´Ðµ cj - Ð²ÐµÑÐµÑÑÐ²ÐµÐ½Ð½ÑÐ¹ ÑÐºÐ°Ð»ÑÑ. @@ -309,7 +309,7 @@ Ð³Ð´Ðµ t, x, xdot, ires, rpar, ipar Ð¸Ð¼ÐµÑÑ Ð°Ð½Ð°Ð»Ð¾Ð³Ð¸ÑÐ½Ð¾Ðµ Ð¾Ð¿ÑÐµÐ´ÐµÐ»ÐµÐ½Ð¸Ðµ, ÑÑÐ¾ Ð¸ Ð²ÑÑÐµ, - r -- Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ·ÑÐ»ÑÑÐ°ÑÐ¾Ð² + r - Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ·ÑÐ»ÑÑÐ°ÑÐ¾Ð² @@ -393,9 +393,9 @@ Ð³Ð´Ðµ t, x, xdot, ires, rpar, ipar Ð¸Ð¼ÐµÑÑ Ð°Ð½Ð°Ð»Ð¾Ð³Ð¸ÑÐ½Ð¾Ðµ Ð¾Ð¿ÑÐµÐ´ÐµÐ»ÐµÐ½Ð¸Ðµ, ÑÑÐ¾ Ð¸ Ð²ÑÑÐµ, - ng -- ÐºÐ¾Ð»Ð¸ÑÐµÑÑÐ²Ð¾ Ð¿Ð¾Ð²ÐµÑÑÐ½Ð¾ÑÑÐµÐ¹, - nx -- ÑÐ°Ð·Ð¼ÐµÑÐ½Ð¾ÑÑÑ ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ Ð¸ - r -- Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ·ÑÐ»ÑÑÐ°ÑÐ¾Ð². + ng - ÐºÐ¾Ð»Ð¸ÑÐµÑÑÐ²Ð¾ Ð¿Ð¾Ð²ÐµÑÑÐ½Ð¾ÑÑÐµÐ¹, + nx - ÑÐ°Ð·Ð¼ÐµÑÐ½Ð¾ÑÑÑ ÑÐ¾ÑÑÐ¾ÑÐ½Ð¸Ñ Ð¸ + r - Ð¼Ð°ÑÑÐ¸Ð² ÑÐµÐ·ÑÐ»ÑÑÐ°ÑÐ¾Ð². @@ -407,8 +407,8 @@ Ð²ÐµÐºÑÐ¾Ñ Ñ Ð´Ð²ÑÐ¼Ñ ÑÐ»ÐµÐ¼ÐµÐ½ÑÐ°Ð¼Ð¸ [times num], Ð³Ð´Ðµ - times -- Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ Ð¼Ð¾Ð¼ÐµÐ½ÑÐ° Ð²ÑÐµÐ¼ÐµÐ½Ð¸ Ð¿ÐµÑÐµÑÐµÑÐµÐ½Ð¸Ñ - Ð¿Ð¾Ð²ÐµÑÑÐ½Ð¾ÑÑÐ¸, num -- ÑÐ¸ÑÐ»Ð¾ Ð¿ÐµÑÐµÑÐµÑÑÐ½Ð½ÑÑ + times - Ð·Ð½Ð°ÑÐµÐ½Ð¸Ðµ Ð¼Ð¾Ð¼ÐµÐ½ÑÐ° Ð²ÑÐµÐ¼ÐµÐ½Ð¸ Ð¿ÐµÑÐµÑÐµÑÐµÐ½Ð¸Ñ + Ð¿Ð¾Ð²ÐµÑÑÐ½Ð¾ÑÑÐ¸, num - ÑÐ¸ÑÐ»Ð¾ Ð¿ÐµÑÐµÑÐµÑÑÐ½Ð½ÑÑ Ð¿Ð¾Ð²ÐµÑÑÐ½Ð¾ÑÑÐµÐ¹ @@ -433,8 +433,8 @@ ,ÑÐ¾ ÐºÐ°Ð¶Ð´ÑÐ¹ ÑÑÐ¾Ð»Ð±ÐµÑ ÑÐ²Ð»ÑÐµÑÑÑ Ð²ÐµÐºÑÐ¾ÑÐ¾Ð¼ Ð²Ð¸Ð´Ð° [t;x(t);xdot(t)], - Ð³Ð´Ðµ t -- Ð¸Ð½Ð´ÐµÐºÑ Ð²ÑÐµÐ¼ÐµÐ½Ð¸ Ð´Ð»Ñ ÐºÐ¾ÑÐ¾ÑÐ¾Ð³Ð¾ Ð²ÑÑÐ¸ÑÐ»ÐµÐ½Ð¾ - ÑÐµÑÐµÐ½Ð¸Ðµ. Ð Ð¿ÑÐ¾ÑÐ¸Ð²Ð½Ð¾Ð¼ ÑÐ»ÑÑÐ°Ðµ y -- Ð²ÐµÐºÑÐ¾Ñ Ð²Ð¸Ð´Ð° + Ð³Ð´Ðµ t - Ð¸Ð½Ð´ÐµÐºÑ Ð²ÑÐµÐ¼ÐµÐ½Ð¸ Ð´Ð»Ñ ÐºÐ¾ÑÐ¾ÑÐ¾Ð³Ð¾ Ð²ÑÑÐ¸ÑÐ»ÐµÐ½Ð¾ + ÑÐµÑÐµÐ½Ð¸Ðµ. Ð Ð¿ÑÐ¾ÑÐ¸Ð²Ð½Ð¾Ð¼ ÑÐ»ÑÑÐ°Ðµ y - Ð²ÐµÐºÑÐ¾Ñ Ð²Ð¸Ð´Ð° [x(t);xdot(t)]. @@ -472,71 +472,133 @@ ÐÑÐ¸Ð¼ÐµÑÑ - - //ÐÑÐ¸Ð¼ÐµÑ Ñ ÐºÐ¾Ð´Ð¾Ð¼ Scilab - function [r,ires]=chemres(t,y,yd) - r(1) = -0.04*y(1) + 1d4*y(2)*y(3) - yd(1); - r(2) = 0.04*y(1) - 1d4*y(2)*y(3) - 3d7*y(2)*y(2) - yd(2); - r(3) = y(1) + y(2) + y(3)-1; - ires = 0; - endfunction - function pd=chemjac(x,y,yd,cj) - pd=[-0.04-cj , 1d4*y(3) , 1d4*y(2); - 0.04 ,-1d4*y(3)-2*3d7*y(2)-cj ,-1d4*y(2); - 1 , 1 , 1 ] - endfunction + +%DAEOPTIONS=list([],1,[],[],[],0,0); // Ð¿ÑÐ¾ÑÐ¸Ñ Ð²ÐµÑÐ½ÑÑÑ ÑÐµÑÐºÑ ÑÐ¾ÑÐµÐº dae +y=dae([x0,xd0],0,4d10,chemres); // Ð±ÐµÐ· ÑÐºÐ¾Ð±Ð¸Ð°Ð½Ð° +y=dae([x0,xd0],0,4d10,chemres,chemjac); // Ñ ÑÐºÐ¾Ð±Ð¸Ð°Ð½Ð¾Ð¼ + ]]> + %t + [btn] = messagebox(["ÐÐ»Ñ ÑÑÐ¾Ð³Ð¾ Ð¿ÑÐ¸Ð¼ÐµÑÐ° ÐÐ°Ð¼ Ð½ÐµÐ¾Ð±ÑÐ¾Ð´Ð¸Ð¼ C-ÐºÐ¾Ð¼Ð¿Ð¸Ð»ÑÑÐ¾Ñ.";"ÐÑÐ¿Ð¾Ð»Ð½ÐµÐ½Ð¸Ðµ ÑÑÐ¾Ð³Ð¾ Ð¿ÑÐ¸Ð¼ÐµÑÐ° Ð¾ÑÐ¼ÐµÐ½ÐµÐ½Ð¾."], "ÐÑÐ¾Ð±Ð»ÐµÐ¼Ð° Ñ Ð¿ÑÐ¾Ð³ÑÐ°Ð¼Ð¼Ð½ÑÐ¼ Ð¾Ð±ÐµÑÐ¿ÐµÑÐµÐ½Ð¸ÐµÐ¼", 'info'); + return +end + +//-1- ÑÐ¾Ð·Ð´Ð°ÑÑ C-ÐºÐ¾Ð´ Ð² TMPDIR - Ð£ÑÐ°Ð²Ð½ÐµÐ½Ð¸Ðµ ÐÐ°Ð½Ð´ÐµÑÐ¿Ð¾Ð»Ð°, Ð½ÐµÑÐ²Ð½Ð°Ñ ÑÐ¾ÑÐ¼Ð° +code=['#include ' + 'void res22(double *t,double *y,double *yd,double *res,int *ires,double *rpar,int *ipar)' + '{res[0] = yd[0] - y[1];' + ' res[1] = yd[1] - (100.0*(1.0 - y[0]*y[0])*y[1] - y[0]);}' + ' ' + 'void jac22(double *t,double *y,double *yd,double *pd,double *cj,double *rpar,int *ipar)' + '{pd[0]=*cj - 0.0;' + ' pd[1]= - (-200.0*y[0]*y[1] - 1.0);' + ' pd[2]= - 1.0;' + ' pd[3]=*cj - (100.0*(1.0 - y[0]*y[0]));}' + ' ' + 'void gr22(int *neq, double *t, double *y, int *ng, double *groot, double *rpar, int *ipar)' + '{ groot[0] = y[0];}'] +previous_dir = pwd(); +cd TMPDIR; +mputl(code, 't22.c') + +//-2- ÑÐºÐ¾Ð¼Ð¿Ð¸Ð»Ð¸ÑÐ¾Ð²Ð°ÑÑ Ð¸ Ð·Ð°Ð³ÑÑÐ·Ð¸ÑÑ ÐµÐ³Ð¾ +ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c','Makefile','t22loader.sce'); +exec('t22loader.sce') + +//-3- Ð·Ð°Ð¿ÑÑÑÐ¸ÑÑ +rtol=[1.d-6;1.d-6]; +atol=[1.d-6;1.d-4]; +t0=0;y0=[2;0]; +y0d=[0;-2]; +t=[20:20:200]; +ng=1; + +//Ð¿ÑÐ¾ÑÑÐ¾Ðµ Ð¼Ð¾Ð´ÐµÐ»Ð¸ÑÐ¾Ð²Ð°Ð½Ð¸Ðµ +t=0:0.003:300; +yy=dae([y0,y0d],t0,t,atol,rtol,'res22','jac22'); +clf();plot(yy(1,:),yy(2,:)) + +// Ð½Ð°Ð¹ÑÐ¸ Ð¿ÐµÑÐ²ÑÑ ÑÐ¾ÑÐºÑ, Ð³Ð´Ðµ yy(1)=0 +[yy,nn,hotd]=dae("root",[y0,y0d],t0,300,atol,rtol,'res22','jac22',ng,'gr22'); +plot(yy(1,1),yy(2,1),'r+') +xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))) +t=0:0.003:300; +yy=dae([y0,y0d],t0,t,atol,rtol,'res22','jac22'); +clf();plot(yy(1,:),yy(2,:)) + +// Ð³Ð¾ÑÑÑÐ¸Ð¹ Ð¿ÐµÑÐµÐ·Ð°Ð¿ÑÑÐº Ð´Ð»Ñ ÑÐ»ÐµÐ´ÑÑÑÐµÐ¹ ÑÐ¾ÑÐºÐ¸ +t01=nn(1); +[pp,qq]=size(yy); +y01=yy(2:3,qq); +y0d1=yy(3:4,qq); +[yy,nn,hotd]=dae("root",[y01,y0d1],t01,300,atol,rtol,'res22','jac22',ng,'gr22',hotd); +plot(yy(1,1),yy(2,1),'r+') +xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))) +cd(previous_dir); +]]> + ' + 'void res22(double *t,double *y,double *yd,double *res,int *ires,double *rpar,int *ipar)' + '{res[0] = yd[0] - y[1];' + ' res[1] = yd[1] - (100.0*(1.0 - y[0]*y[0])*y[1] - y[0]);}' + ' ' + 'void jac22(double *t,double *y,double *yd,double *pd,double *cj,double *rpar,int *ipar)' + '{pd[0]=*cj - 0.0;' + ' pd[1]= - (-200.0*y[0]*y[1] - 1.0);' + ' pd[2]= - 1.0;' + ' pd[3]=*cj - (100.0*(1.0 - y[0]*y[0]));}' + ' ' + 'void gr22(int *neq, double *t, double *y, int *ng, double *groot, double *rpar, int *ipar)' + '{ groot[0] = y[0];}'] +previous_dir = pwd(); +cd TMPDIR; +mputl(code, 't22.c') +ilib_for_link(['res22' 'jac22' 'gr22'],'t22.c',[],'c','Makefile','t22loader.sce'); +exec('t22loader.sce') +rtol=[1.d-6;1.d-6]; +atol=[1.d-6;1.d-4]; +t0=0;y0=[2;0]; +y0d=[0;-2]; +t=[20:20:200]; +ng=1; +t=0:0.003:300; +yy=dae([y0,y0d],t0,t,atol,rtol,'res22','jac22'); +clf();plot(yy(1,:),yy(2,:)) +[yy,nn,hotd]=dae("root",[y0,y0d],t0,300,atol,rtol,'res22','jac22',ng,'gr22'); +plot(yy(1,1),yy(2,1),'r+') +xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))) +t01=nn(1); +[pp,qq]=size(yy); +y01=yy(2:3,qq); +y0d1=yy(3:4,qq); +[yy,nn,hotd]=dae("root",[y01,y0d1],t01,300,atol,rtol,'res22','jac22',ng,'gr22',hotd); +plot(yy(1,1),yy(2,1),'r+') +xstring(yy(1,1)+0.1,yy(2,1),string(nn(1))); +cd(previous_dir); + ]]> Ð¡Ð¼Ð¾ÑÑÐ¸ÑÐµ ÑÐ°ÐºÐ¶Ðµ