A differantial equaton to be solved using matlab

[HP007]

I want to solve the following differential equation using MATLAB and thus want to plot between x & y.
dy/dx=(l*siny-b)/(sinx*cos(y-x))
boundary condition:at x=0;y=sine inverse of(b/l);

 

[da_nang]

Solving it in that form is a bit nasty with Matlab (e.g. with ode45) due to the possibility of discontinuities. I'd suggest multiplying both sides by the denominator and solving it using an implicit method.