Matematické Fórum

Dobrý den, jsem začátečník v MatLabu a mám problémy s vykreslením grafu v následujicím skriptu. Zřejmě je problém s geometrií a goniometrickými funkcemi. Stále se dopočítávám k hodnotám "NaN". Byl bych rád, kdyby mi někdo opravil chybu, která mi znemožňuje vykreslení grafu. K výpočtům jsem použil tuto stránku: http://www.spszr.cz/~blazicek/Projekt/k … is_kyv.htm a zadání: http://www.sdilej.eu/pics/a6b017e87ada2 … 194425.jpg Předem děkuji.


function kulmech(omega)
r = 0.5;
a = 2;
b = 3;
c = 2;
pole = zeros(360, 3);

for i=0:1:360;
    L=sqrt(r*r+a*a-2*r*a*cosd(i));
    beta = asind(r * (sind(i)/L));
    R = b * tand(beta);
   
    %rychlost
    v0=omega*r;
    va = v0 * sind(i + beta);
    vb = va * (R/L);
    v = vb * cosd(i);
    pole(i + 1, 1) = v;
   
    %zrychleni
    omega1 = vb/R;
    anl = omega1 * omega1 * R;
    a = anl * sind(beta);
    pole(i + 1, 2) = a;
   
    %poloha
    h = (-1) * R * sind(beta);
    pole (i + 1, 3) = h;
end

plot(pole)