Tutorial4_ODE_solvers

5 files changed, 63 insertions, 4 deletions

diff --git a/Tutorial4_ODE/Scilab_code/Tutorial4_ode_matrix.sce b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_matrix.sce
new file mode 100644
index 0000000..fff34eb
--- /dev/null
+++ b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_matrix.sce
@@ -0,0 +1,19 @@
+//This script demonstrates the use of ODE solver when matrix is involved
+
+//Definition of the function
+function ydot = func(t,y)
+  ydot = A*y + ones(3,1);
+endfunction
+
+//Defining the matrix
+A = [1 2 1;9 6 1;1 3 4];
+//Initial condition
+y0 = [-1;2;3];
+//Start time
+t0= 0;
+//The time at which the solution is computed
+t = 0.2;
+//Calling the ode solver 
+sol=ode(y0,t0,t,func);
+//Displaying the solution
+disp(sol,"y(t) at t = 0.2");
diff --git a/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple.sce b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple.sce
new file mode 100644
index 0000000..ad42d4d
--- /dev/null
+++ b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple.sce
@@ -0,0 +1,19 @@
+//This script demonstrates the use of ODE solver
+
+//Definition the function. 
+//The function is dy/dt = cos(t)*sin(t) - tan(t) + 1
+function ydot = func(t,y)
+  ydot = cos(t)*sin(t) - tan(t) + 1
+endfunction
+
+//Initial condition of the problem, a scalar or vector
+y0 = -1;
+//Initial time, a real scalar
+t0= 0;
+//The time at which the solution is computed
+t = 0.2;
+//Calling the ode solver 
+sol=ode(y0,t0,t,func);
+//Displaying the solution
+disp(sol,"y(t) at t = 0.2");
+
diff --git a/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple_plot.sce b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple_plot.sce
new file mode 100644
index 0000000..4ecf87e
--- /dev/null
+++ b/Tutorial4_ODE/Scilab_code/Tutorial4_ode_simple_plot.sce
@@ -0,0 +1,19 @@
+//This script demonstrates the use of ODE solver and computing the 
+//solution at different times. The output is in the form of a plot 
+//of time versus computed solution
+
+//Definition the function. 
+function ydot = func(t,y)
+  ydot = t^2*exp(-2*t) + y
+endfunction
+
+//Initial condition
+y0 = -1;
+//Start time
+t0= 0;
+//The array of time where the solution is computed
+t = 4:0.01:10;
+//Calling the ODE solver
+sol=ode(y0,t0,t,func);
+//Plotting the result
+plot(t,sol)
diff --git a/Tutorial5_Nonlinear_equation/Scilab_code/Tutotial5_nonlinear_equation.sce b/Tutorial5_Nonlinear_equation/Scilab_code/Tutotial5_nonlinear_equation.sce
index 6d38c84..5da2894 100644
--- a/Tutorial5_Nonlinear_equation/Scilab_code/Tutotial5_nonlinear_equation.sce
+++ b/Tutorial5_Nonlinear_equation/Scilab_code/Tutotial5_nonlinear_equation.sce
@@ -5,16 +5,14 @@ clc
 
 
 //The nonlinear equation for which solutions are sought
-function y = nonlinear_func(x)
-  y = cos(x)*sin(x) - tan(x) + 1
-endfunction
+exec func.sci;
 
 
 //Initial guess for the solution 
 x0 = 0;
 
 //Computation of solution using fsolve
-y_result = fsolve(x0,nonlinear_func);
+y_result = fsolve(x0,func);
 //Display the solution in command window
 disp(y_result,'Solution of the equation')
 
diff --git a/Tutorial5_Nonlinear_equation/Scilab_code/func.sci b/Tutorial5_Nonlinear_equation/Scilab_code/func.sci
new file mode 100644
index 0000000..91b9fcc
--- /dev/null
+++ b/Tutorial5_Nonlinear_equation/Scilab_code/func.sci
@@ -0,0 +1,4 @@
+//This is script to define a fuction
+function y = func(x)
+  y = cos(x)*sin(x) - tan(x) + 1
+endfunction