Mathematica in the classroom: new tools for exploring precalculus and differential calculus

Prémio de Melhor Artigo de Jovem Investigador atribuído pela empresa Timberlake, apresentado na 1ª Conferência Nacional sobre Computação Simbólica no Ensino e na Investigação - CSEI2012, que decorreu no IST nos dias 2 e 3 de Abril.

A time frequency approach to blind deconvolution in multipath underwater channels

Blind deconvolution is studied in the underwater acoustic channel context, by time-frequency (TF) processing. The acoustic propagation environment is modelled by ray tracing and mathematically described by a multipath propagation channel. Representation of the received signal by means of a signal-dependent TF distribution (radially Gaussian kernel distribution) allowed to visualize the resolved replicas of the emitted signal, while signi cantly attenuating the inherent interferences of classic quadratic TF distributions. The source signal instantaneous frequency estimation was the starting point for both source and channel estimation. Source signature estimation was performed by either TF inversion, based on the Wigner-Ville distribution of the received signal, or a subspace- -based method. The channel estimate was obtained either via a TF formulation of the conventional matched- lter, or via matched- - ltering with the previously obtained source estimate. A shallow water realistic scenario is considered, comprising a 135-m depth water column and an acoustic source located at 90-m depth and 5.6-km range from the receiver. For the corresponding noiseless simulated data, the quality of the best estimates was 0.856 for the source signal, and 0.9664 and 0.9996 for the amplitudes and time-delays of the impulse response, respectively. Application of the proposed deconvolution method to real data of the INTIMATE '96 sea trial conduced to source and channel estimates with the quality of 0.530 and 0.843, respectively. TF processing has proved to remove the typical ill-conditioning of single sensor deterministic deconvolution techniques.

A comunicação matemática escrita no 10º ano de escolaridade em contexto de trabalho de grupo

Relatório da Prática de Ensino Supervisionada, Mestrado em Ensino da Matemática, Universidade de Lisboa, 2015

A new MILP based approach for unit commitment in power production planning

This paper presents a complete, quadratic programming formulation of the standard thermal unit commitment problem in power generation planning, together with a novel iterative optimisation algorithm for its solution. The algorithm, based on a mixed-integer formulation of the problem, considers piecewise linear approximations of the quadratic fuel cost function that are dynamically updated in an iterative way, converging to the optimum; this avoids the requirement of resorting to quadratic programming, making the solution process much quicker. From extensive computational tests on a broad set of benchmark instances of this problem, the algorithm was found to be flexible and capable of easily incorporating different problem constraints. Indeed, it is able to tackle ramp constraints, which although very important in practice were rarely considered in previous publications. Most importantly, optimal solutions were obtained for several well-known benchmark instances, including instances of practical relevance, that are not yet known to have been solved to optimality. Computational experiments and their results showed that the method proposed is both simple and extremely effective.

Mono(Eta(5)-cyclopentadienyl)metal(II) complexes with thienyl acetylide chromophores: synthesis, electrochemical studies, and first hyperpolarizabilities

A series of mono(eta(5)-cyclopentadienyl)metal-(II) complexes with nitro-substituted thienyl acetylide ligands of general formula [M(eta(5)-C5H5)(L)(C C{C4H2S}(n)NO2)] (M = Fe, L = kappa(2)-DPPE, n = 1,2; M = Ru, L = kappa(2)-DPPE, 2 PPh3, n = 1, 2; M = Ni, L = PPh3, n = 1, 2) has been synthesized and fully characterized by NMR, FT-IR, and UV-Vis spectroscopy. The electrochemical behavior of the complexes was explored by cyclic voltammetry. Quadratic hyperpolarizabilities (beta) of the complexes have been determined by hyper-Rayleigh scattering (HRS) measurements at 1500 nm. The effect of donor abilities of different organometallic fragments on the quadratic hyperpolarizabilities was studied and correlated with spectroscopic and electrochemical data. Density functional theory (DFT) and time-dependent DFT (TDDFT) calculations were employed to get a better understanding of the second-order nonlinear optical properties in these complexes. In this series, the complexity of the push pull systems is revealed; even so, several trends in the second-order hyperpolarizability can still be recognized. In particular, the overall data seem to indicate that the existence of other electronic transitions in addition to the main MLCT clearly controls the effectiveness of the organometallic donor ability on the second-order NLO properties of these push pull systems.

Localização Ótima de Aparelhos de Corte Normalmente Abertos e Normalmente Fechados em Redes de Distribuição

Tipicamente as redes elétricas de distribuição apresentam uma topologia parcialmente malhada e são exploradas radialmente. A topologia radial é obtida através da abertura das malhas nos locais que otimizam o ponto de operação da rede, através da instalação de aparelhos de corte que operam normalmente abertos. Para além de manterem a topologia radial, estes equipamentos possibilitam também a transferência de cargas entre saídas, aquando da ocorrência de defeitos. As saídas radiais são ainda dotadas de aparelhos de corte que operam normalmente fechados, estes têm como objetivo maximizar a fiabilidade e isolar defeitos, minimizando a área afetada pelos mesmos. Assim, na presente dissertação são desenvolvidos dois algoritmos determinísticos para a localização ótima de aparelhos de corte normalmente abertos e fechados, minimizando a potência ativa de perdas e o custo da energia não distribuída. O algoritmo de localização de aparelhos de corte normalmente abertos visa encontrar a topologia radial ótima que minimiza a potência ativa de perdas. O método é desenvolvido em ambiente Matlab – Tomlab, e é formulado como um problema de programação quadrática inteira mista. A topologia radial ótima é garantida através do cálculo de um trânsito de potências ótimo baseado no modelo DC. A função objetivo é dada pelas perdas por efeito de Joule. Por outro lado o problema é restringido pela primeira lei de Kirchhoff, limites de geração das subestações, limites térmicos dos condutores, trânsito de potência unidirecional e pela condição de radialidade. Os aparelhos de corte normalmente fechados são localizados ao longo das saídas radiais obtidas pelo anterior algoritmo, e permite minimizar o custo da energia não distribuída. No limite é possível localizar um aparelho de corte normalmente fechado em todas as linhas de uma rede de distribuição, sendo esta a solução que minimiza a energia não distribuída. No entanto, tendo em conta que a cada aparelho de corte está associado um investimento, é fundamental encontrar um equilíbrio entre a melhoria de fiabilidade e o investimento. Desta forma, o algoritmo desenvolvido avalia os benefícios obtidos com a instalação de aparelhos de corte normalmente fechados, e retorna o número e a localização dos mesmo que minimiza o custo da energia não distribuída. Os métodos apresentados são testados em duas redes de distribuição reais, exploradas com um nível de tensão de 15 kV e 30 kV, respetivamente. A primeira rede é localizada no distrito do Porto e é caraterizada por uma topologia mista e urbana. A segunda rede é localizada no distrito de Bragança e é caracterizada por uma topologia maioritariamente aérea e rural.

Otimização de esquemas terapêuticos do carcinoma da próstata envolvendo braquiterapia de baixa taxa de dose e radioterapia externa : abordagens física e radiobiológica

RESUMO: Este trabalho teve como objetivo a determinação de esquemas de tratamento alternativos para o carcinoma da próstata com radioterapia externa (EBRT) e braquiterapia de baixa taxa de dose (LDRBT) com implantes permanentes de Iodo-125, biologicamente equivalentes aos convencionalmente usados na prática clínica, com recurso a modelos teóricos e a métodos de Monte Carlo (MC). Os conceitos de dose biológica efetiva (BED) e de dose uniforme equivalente (EUD) foram utilizados, com o modelo linear-quadrático (LQ), para a determinação de regimes de tratamento equivalentes. Numa primeira abordagem, utilizou-se a BED para determinar: 1) esquemas hipofracionados de EBRT mantendo as complicações retais tardias de regimes convencionais com doses totais de 75,6 Gy, 77,4 Gy, 79,2 Gy e 81,0 Gy; e 2) a relação entre as doses totais de EBRT e LDRBT de modo a manter a BED do regime convencional de 45 Gy de EBRT e 110 Gy de LDRBT. Numa segunda abordagem, recorreu-se ao código de MC MCNPX para a simulação de distribuições de dose de EBRT e LDRBT em dois fantomas de voxel segmentados a partir das imagens de tomografia computorizada de pacientes com carcinoma da próstata. Os resultados das simulações de EBRT e LDRBT foram somados e determinada uma EUD total de forma a obterem-se: 1) esquemas equivalentes ao tratamento convencional de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT; e 2) esquemas equivalentes a EUD na próstata de 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy e 110 Gy. Em todos os resultados nota-se um ganho terapêutico teórico na utilização de esquemas hipofracionados de EBRT. Para uma BED no reto equivalente ao esquema convencional, tem-se um aumento de 2% na BED da próstata com menos 5 frações. Este incremento dá-se de forma cada vez mais visível à medida que se reduz o número de frações, sendo da ordem dos 10-11% com menos 20 frações e dos 35-45% com menos 40 frações. Considerando os resultados das simulações de EBRT, obteve-se uma EUD média de 107 Gy para a próstata e de 42 Gy para o reto, com o esquema convencional de 110 Gy de LDRBT, seguidos de 25 frações de 1,8 Gy de EBRT. Em termos de probabilidade de controlo tumoral (igual EUD), é equivalente a este tratamento a administração de EBRT em 66 frações de 1,8 Gy, 56 de 2 Gy, 40 de 2,5 Gy, 31 de 3 Gy, 20 de 4 Gy ou 13 de 5 Gy. Relativamente à administração de 66 frações de 1,8 Gy, a EUD generalizada no reto reduz em 6% com o recurso a frações de 2,5 Gy e em 10% com frações de 4 Gy. Determinou-se uma BED total de 162 Gy para a administração de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT. Variando-se a dose total de LDRBT (TDLDRBT) em função da dose total de EBRT (TDEBRT), de modo a garantir uma BED de 162 Gy, obteve-se a seguinte relação:.......... Os resultados das simulações mostram que a EUD no reto diminui com o aumento da dose total de LDRBT para dose por fração de EBRT (dEBRT) inferiores a 2, Gy e aumenta para dEBRT a partir dos 3 Gy. Para quantidades de TDLDRBT mais baixas (<50 Gy), o reto beneficia de frações maiores de EBRT. À medida que se aumenta a TDLDRBT, a EUD generalizada no reto torna-se menos dependente da dEBRT. Este trabalho mostra que é possível a utilização de diferentes regimes de tratamento para o carcinoma da próstata com radioterapia que possibilitem um ganho terapêutico, quer seja administrando uma maior dose biológica com efeitos tardios constantes, quer mantendo a dose no tumor e diminuindo a toxicidade retal. A utilização com precaução de esquemas hipofracionados de EBRT, para além do benefício terapêutico, pode trazer vantagens ao nível da conveniência para o paciente e economia de custos. Os resultados das simulações deste estudo e conversão para doses de efeito biológico para o tratamento do carcinoma da próstata apresentam linhas de orientação teórica de interesse para novos ensaios clínicos. --------------------------------------------------ABSTRACT: The purpose of this work was to determine alternative radiotherapy regimens for the treatment of prostate cancer using external beam radiotherapy (EBRT) and low dose-rate brachytherapy (LDRBT) with Iodine-125 permanent implants which are biologically equivalent to conventional clinical treatments, by the use of theoretical models and Monte Carlo techniques. The concepts of biological effective dose (BED) and equivalent uniform dose (EUD), together with the linear-quadratic model (LQ), were used for determining equivalent treatment regimens. In a first approach, the BED concept was used to determine: 1) hypofractionated schemes of EBRT maintaining late rectal complications as with the conventional regimens with total doses of 75.6 Gy, 77.4 Gy, 79.2 Gy and 81.0 Gy; and 2) the relationship between total doses of EBRT and LDRBT in order to keep the BED of the conventional treatment of 45 Gy of EBRT and 110 Gy of LDRBT. In a second approach, the MC code MCNPX was used for simulating dose distributions of EBRT and LDRBT in two voxel phantoms segmented from the computed tomography of patients with prostate cancer. The results of the simulations of EBRT and LDRBT were added up and given an overall EUD in order to obtain: 1) equivalent to conventional treatment regimens of 25 fraction of 1.8 Gy of EBRT in combination with 110Gy of LDRBT; and 2) equivalent schemes of EUD of 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy, and 110Gy to the prostate. In all the results it is noted a therapeutic gain using hypofractionated EBRT schemes. For a rectal BED equivalent to the conventional regimen, an increment of 2% in the prostate BED was achieved with less 5 fractions. This increase is visibly higher as the number of fractions decrease, amounting 10-11% with less 20 fractions and 35-45% with less 20 fractions. Considering the results of the EBRT simulations an average EUD of 107 Gy was achieved for the prostate and of 42 Gy for the rectum with the conventional scheme of 110 Gy of LDRBT followed by 25 fractions of 1.8 Gy of EBRT. In terms of tumor control probability (same EUD) it is equivalent to this treatment, for example, delivering the EBRT in 66 fractions of 1.8 Gy, 56 fractions of 2 Gy, 40 fractions of 2.5 Gy, 31 fractions of 3 Gy, 20 fractions of 4 Gy or 13 fractions of 5 Gy. Regarding the use of 66 fractions of 1.8 Gy, the rectum EUD is reduced to 6% with 2.5 Gy per fraction and to 10% with 4 Gy. A total BED of 162 Gy was achieved for the delivery of 25 fractions of 1.8 Gy of EBRT in combination with 110 Gy of LDRBT. By varying the total dose of LDRBT (TDLDRBT) with the total dose of EBRT (TDEBRT) so as to ensure a BED of 162 Gy, the following relationship was obtained: ....... The simulation results show that the rectum EUD decreases with the increase of the TDLDRBT, for EBRT dose per fracion (dEBRT) less than 2.5 Gy and increases for dEBRT above 3 Gy. For lower amounts of TDLDRBT (< 50Gy), the rectum benefits of larger EBRT fractions. As the TDLDRBT increases, the rectum gEUD becomes less dependent on the dEBRT. The use of different regimens which enable a therapeutic gain, whether deivering a higher dose with the same late biological effects or maintaining the dose to the tumor and reducing rectal toxicity is possible. The use with precaution of hypofractionated regimens, in addition to the therapeutic benefit, can bring advantages in terms of convenience for the patient and cost savings. The simulation results of this study together with the biological dose conversion for the treatment of prostate cancer serve as guidelines of interest for new clinical trials.

Associations between alcohol consumption and selected cytokines in a Swiss population-based sample (CoLaus study).

To assess the associations between alcohol consumption and cytokine levels (interleukin-1beta - IL-1β; interleukin-6 - IL-6 and tumor necrosis factor-α - TNF-α) in a Caucasian population. Population sample of 2884 men and 3201 women aged 35-75. Alcohol consumption was categorized as nondrinkers, low (1-6 drinks/week), moderate (7-13/week) and high (14+/week). No difference in IL-1β levels was found between alcohol consumption categories. Low and moderate alcohol consumption led to lower IL-6 levels: median (interquartile range) 1.47 (0.70-3.51), 1.41 (0.70-3.32), 1.42 (0.66-3.19) and 1.70 (0.83-4.39) pg/ml for nondrinkers, low, moderate and high drinkers, respectively, p<0.01, but this association was no longer significant after multivariate adjustment. Compared to nondrinkers, moderate drinkers had the lowest odds (Odds ratio=0.86 (0.71-1.03)) of being in the highest quartile of IL-6, with a significant (p<0.05) quadratic trend. Low and moderate alcohol consumption led to lower TNF-α levels: 2.92 (1.79-4.63), 2.83 (1.84-4.48), 2.82 (1.76-4.34) and 3.15 (1.91-4.73) pg/ml for nondrinkers, low, moderate and high drinkers, respectively, p<0.02, and this difference remained borderline significant (p=0.06) after multivariate adjustment. Moderate drinkers had a lower odds (0.81 [0.68-0.98]) of being in the highest quartile of TNF-α. No specific alcoholic beverage (wine, beer or spirits) effect was found. Moderate alcohol consumption is associated with lower levels of IL-6 and (to a lesser degree) of TNF-α, irrespective of the type of alcohol consumed. No association was found between IL-1β levels and alcohol consumption.

Admissions of patients with alcohol intoxication in the Emergency Department: a growing phenomenon.

QUESTION UNDER STUDY: To investigate the change over time in the number of ED admissions with positive blood alcohol concentration (BAC) and to evaluate predictors of BAC level. METHODS: We conducted a single site retrospective study at the ED of a tertiary referral hospital (western part of Switzerland) and obtained all the BAC performed from 2002 to 2011. We determined the proportion of ED admissions with positive BAC (number of positive BAC/number of admissions). Regression models assessed trends in the proportion of admissions with positive BAC and the predictors of BAC level among patients with positive BAC. RESULTS: A total of 319,489 admissions were recorded and 20,021 BAC tests were performed, of which 14,359 were positive, divided 34.5% female and 65.5% male. The mean (SD) age was 41.7(16.8), and the mean BAC was 2.12(1.04) permille (g of ethanol/liter of blood). An increase in the number of positive BAC was observed, from 756 in 2002 to 1,819 in 2011. The total number of admissions also increased but less: 1.2 versus 2.4 times more admissions with positive BAC. Being male was independently associated with a higher (+0.19 permille) BAC, as was each passing year (+0.03). A significant quadratic association with age indicated a maximum BAC at age 53. CONCLUSION: We observed an increase in the percentage of admissions with positive BAC that was not limited to younger individuals. Given the potential consequences of alcohol intoxication, and the large burden imposed on ED teams, communities should be encouraged to take measures aimed at reducing alcohol intoxication.

Subclinical thyroid dysfunction and the risk of heart failure events: an individual participant data analysis from 6 prospective cohorts.

BACKGROUND: American College of Cardiology/American Heart Association guidelines for the diagnosis and management of heart failure recommend investigating exacerbating conditions such as thyroid dysfunction, but without specifying the impact of different thyroid-stimulation hormone (TSH) levels. Limited prospective data exist on the association between subclinical thyroid dysfunction and heart failure events. METHODS AND RESULTS: We performed a pooled analysis of individual participant data using all available prospective cohorts with thyroid function tests and subsequent follow-up of heart failure events. Individual data on 25 390 participants with 216 248 person-years of follow-up were supplied from 6 prospective cohorts in the United States and Europe. Euthyroidism was defined as TSH of 0.45 to 4.49 mIU/L, subclinical hypothyroidism as TSH of 4.5 to 19.9 mIU/L, and subclinical hyperthyroidism as TSH <0.45 mIU/L, the last two with normal free thyroxine levels. Among 25 390 participants, 2068 (8.1%) had subclinical hypothyroidism and 648 (2.6%) had subclinical hyperthyroidism. In age- and sex-adjusted analyses, risks of heart failure events were increased with both higher and lower TSH levels (P for quadratic pattern <0.01); the hazard ratio was 1.01 (95% confidence interval, 0.81-1.26) for TSH of 4.5 to 6.9 mIU/L, 1.65 (95% confidence interval, 0.84-3.23) for TSH of 7.0 to 9.9 mIU/L, 1.86 (95% confidence interval, 1.27-2.72) for TSH of 10.0 to 19.9 mIU/L (P for trend <0.01) and 1.31 (95% confidence interval, 0.88-1.95) for TSH of 0.10 to 0.44 mIU/L and 1.94 (95% confidence interval, 1.01-3.72) for TSH <0.10 mIU/L (P for trend=0.047). Risks remained similar after adjustment for cardiovascular risk factors. CONCLUSION: Risks of heart failure events were increased with both higher and lower TSH levels, particularly for TSH ≥10 and <0.10 mIU/L.

The surface electromyography-force relationship during isometric dorsiflexion in males and females

This study evaluated sex-related differences in the tibialis anterior (TA) surface electromyography (EMG) to force relationship. One-hundred participants (50 males and 50 females) performed three isometric contractions at 20, 40, 60, 80, and 100% of maximal voluntary contraction (MVC) in an apparatus designed to isolate the action of the dorsiflexors. The surface EMG signal was amplified (lOOOx), band-pass filtered (10-500Hz), and sampled at 2048 Hz. The load cell signal was low-passed filtered at 100 Hz and sampled at the same rate. Males were stronger than females {P <0.05). However, there was no significant difference in root-mean-square (RMS) values between sexes {P <0.05). Both sexes exhibited a quadratic increase in RMS across force levels (P <0.05). The mean power frequency (MNF) for males was greater than for females {P <0.05). Males and females exhibited a linear increase in both frequency measures up to 80% of MVC (P <0.05). Between 80 and 100% MVC, the frequency values for the females plateaued while males showed a decrease {P <0.05). The magnitude of the difference in MNF between males and females was consistent with sex-specific TA physiology. In general, the pattern of means for RMS and MNF between males and females revealed no sex-related differences in the surface EMG/force relationship. We therefore conclude that there are no sex-related differences in the gradation of muscle force.

Cognitive set, absorption, and achievement as moderators of the perceptual closure effect in picture priming /

The present study examined individual differences in Absorption and fantasy, as well as in Achiievement and achievement striving as possible moderators of the perceptual closure effect found by Snodgrass and Feenan (1990). The study also examined whether different instructions (experiential versus instrumental) interact with the personality variables to moderate the relationship between priming and subsequent performance on a picture completion task. 1 28 participants completed two sessions, one to fill out the MPQ and NEO personality inventories and the other to complete the experimental task. The experimental task consisted of a priming phase and a test phase, with pictures presented on a computer screen for both phases. Participants were shown 30 pictures in the priming phase, and then shovm the 30 primed pictures along with 30 new pictures for the test phase. Participants were randomly assigned to receive one of the two different instruction sets for the task. Two measures of performance were calculated, most fragmented measure and threshold. Results of the present study confirm that a five-second exposure time is long enough to produce the perceptual closure effect. The analysis of the two-way interaction effects indicated a significant quadratic interaction of Absorption with priming level on threshold performance. The results were in the opposite direction of predictions. Possible explanations for the Absorption results include lack of optimal conditions, lack of intrinsic motivation and measurement problems. Primary analyses also revealed two significant between-subject effects of fantasy and achievement striving on performance collapsed across priming levels. These results suggest that fantasy has a beneficial effect on performance at test for pictures primed at all levels, whereas achievement striving seems to have an adverse effect on performance at test for pictures primed at all levels. Results of the secondary analyses with a revised threshold performance measure indicated a significant quadratic interaction of Absorption, condition and priming level. In the experiential condition, test performance, based on Absorption scores for pictures primed at level 4, showed a positive slope and performance for pictures primed at levels 1 and 7 based on Absorption showed a negative slope. The reverse effect was found in the instrumental condition. The results suggest that Absorption, in combination with experiential involvement, may affect implicit memory. A second significant result of the secondary analyses was a linear three-way interaction of Achievement, condition and priming level on performance. Results suggest that as Achievement scores increased, test performance improved for less fragmented primed pictures in the instrumental condition and test performance improved for more highly fragmented primes in the experiential condition. Results from the secondary analyses suggest that the revised threshold measure may be more sensitive to individual differences. Results of the exploratory analyses with Openness to Experience, Conscientiousness and agentic positive emotionality (PEM-A) measures indicated no significant effects of any of these personality variables. Results suggest that facets of the scales may be more useful with regard to perceptual research, and that future research should examine narrowly focused personality traits as opposed to broader constructs.

Analysis Of The A11a22 X 11a11 And A 33a22 X 1 electronic transitions in selenocarbonyl difluoride

The absorption spectrum of F2CSe in the 18800-21900 cm-1 region has been recorded at -770 C and 220 C under the conditions of medium resolution. The responsible electronic promotion is TI* + n excitation which leads to 3A2 and lA2 excited states. Progressions in vI', v2', v3" v4' and v4" have been identified in the spectrum and have been analyzed in terms of vibronic transitions between a planar ground state and a nQnplanar excited state. The - 3 - 1 - 1 - 1 origins of the a A2 + X Al and A A2 + X Al systems were assigned to the bands at 19018 cm-l and 19689 cm-l . This has given a singlet-triplet splittl. n g lA2 - 3A2 P f 671 cm -1 The out-of-plane wagging levels were found to be anharmonic. 1 -1 Barrier heights of 2483 cm- and 2923 cm were obtained for the lA2 and 3A2 upper states from a fitting of the energy levels of a Lorentzian-quadratic function to the observed levels in the out-of-plane wagging modes. 1 3 For the A2 and A2 states nonplanar equilibrium angles of 30.10 and 31.40 have been evaluated respectively. i

ReAlM - a system to manipulate relations

Given a heterogeneous relation algebra R, it is well known that the algebra of matrices with coefficient from R is relation algebra with relational sums that is not necessarily finite. When a relational product exists or the point axiom is given, we can represent the relation algebra by concrete binary relations between sets, which means the algebra may be seen as an algebra of Boolean matrices. However, it is not possible to represent every relation algebra. It is well known that the smallest relation algebra that is not representable has only 16 elements. Such an algebra can not be put in a Boolean matrix form.[15] In [15, 16] it was shown that every relation algebra R with relational sums and sub-objects is equivalent to an algebra of matrices over a suitable basis. This basis is given by the integral objects of R, and is, compared to R, much smaller. Aim of my thesis is to develop a system called ReAlM - Relation Algebra Manipulator - that is capable of visualizing computations in arbitrary relation algebras using the matrix approach.

Extending relAPS to first order logic

RelAPS is an interactive system assisting in proving relation-algebraic theorems. The aim of the system is to provide an environment where a user can perform a relation-algebraic proof similar to doing it using pencil and paper. The previous version of RelAPS accepts only Horn-formulas. To extend the system to first order logic, we have defined and implemented a new language based on theory of allegories as well as a new calculus. The language has two different kinds of terms; object terms and relational terms, where object terms are built from object constant symbols and object variables, and relational terms from typed relational constant symbols, typed relational variables, typed operation symbols and the regular operations available in any allegory. The calculus is a mixture of natural deduction and the sequent calculus. It is formulated in a sequent style but with exactly one formula on the right-hand side. We have shown soundness and completeness of this new logic which verifies that the underlying proof system of RelAPS is working correctly.