Solving a differential equation using Scicos

Let us assume the following ordinary differential equation (ODE)

y'''-y''+y'-y+2 = 0

where y(x) is unknown function. We are going to solve this equation numerically using Scilab 5.5.2 and its toolbox Xcos. We assume the following initial conditions: y(0) = 1, y'(0) = -1 and y''(0) = 1.
In general you should solve the equation algebraically with respect to its highest derivative

y''' = y'' - y' + y - 2

and build the Scicos diagram, where you construct the input to a set of integrators based on their outputs. Since we deal with a third-order differential equation, we use the series of three integrators. The possible solution looks as follows

Above we have used the following blocks form the Palettes: INTEGRAL_m, EXPRESSION, CSCOPE and CLOCK_c.
In the block EXPRESSION we set the following parameters:

i.e. the number of inputs as 3 and scilab expression, when we put the formula for y''' assuming that u1=y, u2=y' and u3=y''. The output of this block goes to the input of the first integration block. The outputs of the integration blocks are sent to the corresponding inputs of the EXPRESSION block. The output of the first integral block (from the left) is y'', the output form the second one is y' and the output from the last one is y.
Now we set initial conditions in the integral blocks. The initial condition parameter of a given integrator stands for initial value of its output signal. Thus in the first block from the left we set initial value of the second derivative y''(0) = 1:

Then we do the same with the remaining integration blocks, i.e. in the second one we set the initial value of the first derivative y'(0) = -1, while in the last one we put as initial condition y(0) = 1.

We also use blocks CSCOPE and CLOCK_c in order to show results. In the last one we set the following parameters

when Period means the period of data sampling sent to the scope and Initialisation Time is the initial instance when this data is sampled for the first time. In the project window menu we choose Simulation and Setup. In the window Set Parameters we set Final integration time as 30. Note that the simulation time in the Xcos model is now understood as independent variable x of the function y(x). In the block CSCOPE we set Ymin and Ymax parameters in the following way

We save the model and then we start simulation obtaining the following plot

Let us change the constant term of the equation slightly, i.e. let us add 10^10 to 2. We obtain the following result:

i.e. the system is unstable. Similar result we get in the case of subtraction of 10^10 from 2: