Найти все точки принадлежащие эллипсу

[morg4n]

Задание:

По заданным размерам экрана дисплея в мм и заданной геометрической погрешности вычисляем минимальный геометрической размер пиксела и определяем пиксельное разрешение экрана дисплея.

Осуществляя последовательное прохождение растра от левого верхнего пиксела растра до конца первой строки, затем переходим на вторую строку и т.д. доходим до конца растра (рисунок2).
В каждом пикселе подставляем его геометрическое значение в алгебраическое уравнение и решая его определяем, принадлежит ли данный пиксел нашему эллипсу, если да, то в массив растра записываем 1, если нет, то пишем 0. Запись массива заканчивается в последней точке.

Условия: размер экрана 160 мм ширина, 180мм высота, погрешность 0.13мм (как я понимаю размер пикселя), центр эллипса 80 мм по х , 50мм по y, a 22mm, b 45mm.

#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
double sqr(double some) {return some*some;}
int main()
{
    double x,y;
    int width=160,height=180;
    int x0=80,y0=50,
        a=22,b=45;
    double mis=0.13;

    freopen("rastr_model.txt","w",stdout);
    for(y=height;y>=0;y=y-mis){
        for(x=0;x<=width;x=x+mis){
            if(sqr(x-x0)/sqr(a) + sqr(y-y0)/sqr(b) <1.0+mis
                    && sqr(x-x0)/sqr(a) + sqr(y-y0)/sqr(b) >1.0-mis)
           cout<<"1";
            else
           cout<<"0";
        }
        cout<<"\n";
    }
    return 0;
}

Правильно ли я все понял? Может есть лучшие решения?
И как найти пиксельное разрешение экрана дисплея?

просто получается что-то дикое.. output file over много мб...

[Abstraction]

У Вас в коде x и y не меняются.

[morg4n]

Если Вас не затруднит, подскажите правильную запись цикла.

for(y=height;y>=0;y=y-mis){
        for(x=0;x<=width;x=x+mis){

А так в остальном код нормальный?

[Abstraction]

Если Вас не затруднит, подскажите правильную запись цикла.

Зависит от того, что Вы хотите сделать. Если нужно обойти все пиксели экрана, то скорее так:

for(int y = 0; y<height; ++y){
  for(int x = 0; x<width; ++x){
    //Проверка, принадлежит ли эллипсу пиксель (x,y)
  }
  //Переход на новую строку
}

А так в остальном код нормальный?

Как решение "в лоб" - вроде да. Но сейчас для всех пикселей проводится одинаковая проверка, а ведь большую часть пикселей можно отбросить условиями, которые вычисляются быстрее (к примеру, можно вывести размеры некоторого вписанного в эллипс прямоугольника - принадлежность к нему проверяется в четыре сравнения, без квадратов и делений; можно, напротив, для данного y быстро выделять не больше шести пикселей, которые уже проверять детально).