# Difference between revisions of "Systems of differential equations"

Jump to navigationJump to search

Display solutions of the ordinary differential equation

$y_1'= f_1(x,y_1,y_2)$
$y_2'= f_2(x,y_1,y_2)$

with initial values $(x_0,c_1)$, $(x_0,c_2)$.

f1(x,y1,y2)=
f2(x,y1,y2)=

### The underlying JavaScript code

<form>
f<sub>1</sub>(x,y1,y2)=<input type="text" id="odeinput1" value="y1+y2"><br />
f<sub>2</sub>(x,y1,y2)=<input type="text" id="odeinput2" value="y2+1"><input type=button value="ok" onclick="doIt()">
</form>
var brd = JXG.JSXGraph.initBoard('jxgbox', {axis:true, boundingbox:[-11,11,11,-11]});
var N = brd.create('slider',[[-7,9.5],[7,9.5],[-15,10,15]], {name:'N'});
var P1 = brd.create('point',[1,-1], {name:'(x_0,c_1)'});
var line = brd.create('line',[function(){return -P1.X();},function(){return 1;},function(){return 0;}],{visible:false});
var P2 = brd.create('glider',[1,-0.5,line], {name:'(x_0,c_2)'});

function doIt() {
var txt1 = document.getElementById("odeinput1").value;
var txt2 = document.getElementById("odeinput2").value;

var snip1 = brd.jc.snippet(txt1, true, 'x, y1, y2');
var snip2 = brd.jc.snippet(txt2, true, 'x, y1, y2');
f = function (x, yy) {
return [snip1(x, yy[0], yy[1]), snip2(x, yy[0], yy[1])];
}
brd.update();
}

function ode() {
return JXG.Math.Numerics.rungeKutta('heun', [P1.Y(),P2.Y()], [P1.X(), P1.X()+N.Value()], 200, f);
}

var g1 = brd.create('curve', [[0],[0]], {strokeColor:'red', strokeWidth:2, name:'y_1', withLabel:false});
var g2 = brd.create('curve', [[0],[0]], {strokeColor:'black', strokeWidth:2, name:'y_2', withLabel:false});
g1.updateDataArray = function() {
var data = ode();
var h = N.Value()/200;
var i;

this.dataX = [];
this.dataY = [];
for(i=0; i<data.length; i++) {
this.dataX[i] = P1.X()+i*h;
this.dataY[i] = data[i][0];
}
};
g2.updateDataArray = function() {
var data = ode();
var h = N.Value()/200;
var i;

this.dataX = [];
this.dataY = [];
for(i=0; i<data.length; i++) {
this.dataX[i] = P2.X()+i*h;
this.dataY[i] = data[i][1];
}
};
doIt();