Diferenciální rovnice - Matlab – Matematika – Fórum – Programujte.com

Zdravím,
pokouším se už 2. den v matlabu vyřešit problém této rovnice :
u'' + bu' + sin(u) = 0 ;
u(0) =0
u'(0) = lambda

Mým cílem je dostat závislost parametru b na lambda
Pokouším se to vyřešit Eulerovou metodou, avšak rovnice už je převedená z okrajové a nemám tak podle čeho střílet. Z následujícího kodu vycházím, ale v žádé metodě zatím nevychází závislost jak by měla. Nějaká rada?
 

function v = fce(t,y);

global b 

v=[0;0];
v(1) = y(2);
v(2) = -y(2)-b*sin(y(1));

function main
clear all
clc

global b 

b=3;        % tlumení
 

alfa=0;    % poč. odchylka
beta=3;    % poč úhlová rychlost

options=odeset('AbsTol',10^(-8),'RelTol',10^(-8));
[t,y] = ode45(@fce,[0 100],[alfa;beta],options);