- /
- 
        teal dive
        on 6 Nov 2023
        
        
 
    - 6
- 20
- 0
- 0
- 260
function drawframe(f)
    f = f^0.95
    screen = zeros(200, 200, 3);
        for i = 1:200
            for j = 1:200
                d1 = getDist(i,j,100,100);
                if d1 < d1*f^1.3
                    screen(i,j,3) = sind(d1*f^1.2);
                    screen(i,j,2) = cosd(d1*f^1.4);
                end
            end
        end
    imshow(screen)
end
function d = getDist(x1, y1, x2, y2)
   d = ((x1-x2).^2 + (y1-y2).^2).^0.5;
end


 
 
           