Брюсселятор
independent: { t: 0 }, dependent: { // x: [2,4,3, -1,-0.2,-0.01,-0.001, 0.001,0.01,0.1, -1,-3,-5, -5 ], // y: [2,4,3, -5,-5,-5,-5, -5,-5,-5, 10,10,10, 7] // x: [4, 4, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0, 0.5], // y: [0, 1, 2, 3, 3.5, 0.05, 0.1, 0.5, 1, 2, 3, 4, 4 ] // x: [0, 2], // y: [0, 4] x: [4], y: [4] }, parameters: { A: { min: 1, max: 2, value: 1.5 }, B: { min: 0.01, max: 4.9, value: 3 }, }, derivatives: { dx_dt: 'A + x*x*y-(B+1)*x', dy_dt: 'B * x - x*x*y' }, param_charts: [ ], solver: { step: 0.01, accuracy: 0.001, autostop: 30, stop_any: true, }, plot_options: { // 'x': { min: -5, max: 10, min_stop: -10, max_stop: 30,}, // 'y': { min: -5, max: 10, min_stop: -10, max_stop: 30,}, 'x': { min: 0, max: 10,}, 'y': { min: 0, max: 10,}, }, plots: { 'x(t)': true, 'y(x)': true, 'y(t)': true, }, init_hook: function() { var a = calc.sys_state['B']-1; var b = calc.sys_state['A']*calc.sys_state['A']; var c = -calc.sys_state['B']; var d = -calc.sys_state['A']*calc.sys_state['A']; var s_2 = (a+d)/2; var delta = a*d-b*c; var s2_4 = s_2*s_2; if (s2_4 >= delta) { var l1 = '' + (s_2 - Math.sqrt(s2_4-delta)).toPrecision(3); var l2 = '' + (s_2 + Math.sqrt(s2_4-delta)).toPrecision(3); } else { var l1 = '' + s_2.toPrecision(3) + '+ ' + Math.sqrt(delta-s2_4).toPrecision(3) + ' \\cdot i'; var l2 = '' + s_2.toPrecision(3) + '- ' + Math.sqrt(delta-s2_4).toPrecision(3) + ' \\cdot i'; } $('trax_controller_postfix').innerText = '$$ \\begin{aligned} \\lambda_1 &= {'+l1+'} \\\\ \\lambda_2 &= {'+l2+'} \\end{aligned} $$'; MathJax.Hub.Queue(["Typeset",MathJax.Hub,"trax_controller_postfix"]); },
Счёт:
Запускать по щелчку
Редактировать модель