%structure for solving the ODE for Lorenz Attractor by Ankur Kumar
%source: wikipedia

function a = sol(t,x)
%condition: o,b,p>0
o = 10; %prandtl number
b = 8/3; %chaotic number
p = 99.96; %rayleigh number
dxdt = o*(x(2)-x(1));
dydt = x(1)*p-x(1)*x(3)-x(2);
dzdt = x(1)*x(2) - b*x(3);
a = [dxdt ;dydt; dzdt];