G3 FC1N T1 novo

Universidade Federal do Pará
Centro de Ciências Exatas e Naturais
Departamento de Física
Física Computacional

TRABALHO 1
Exercícios 1.1 e 1.2 do Capítulo 1 da Apostila

Belém
Abril/2014

Universidade Federal do Pará
Centro de Ciências Exatas e Naturais
Departamento de Física
Física Computacional

TRABALHO 1
Exercícios 1.1 e 1.2 do Capítulo 1 da Apostila

Grupo 3
Kleber Barbosa Leal

Matrícula: 12094001201

Rogério de Oliveira Pinto

Matrícula: 12094002901

Wladison W. T. Aboim

Matrícula: 12094002001

Ícaro Pacheco

Matrícula: 07029003201

Elielzer de Souza Nuayed

Matrícula: 09094004101

Vanessa Silva Ferreira

Matrícula: 09094001129

Belém
Abril/2014

1.

Trabalho 1, parte (1.1)

Edite um arquivo de nome “tabela.m” copiando o script (código fonte) para MATLAB apresentado na Figura1.2. No ambiente MATLAB execute o código e verifique a presença do arquivo TAB.DAT no diretório de trabalho, após a execução do script. Abra o arquivo
TAB.DAT com o editor e verifique se número de linhas está correto.
2.1.

Script
%Código fonte do Exemplo 1 em MATLAB/OCTAVE
% tabela.m
% cria TAB.DAT c/ t v e x
Nlin = 1000; dt = 0.1; a =input('Entre c/ a :'); to=input('Entre c/ to:'); vo=input('Entre c/ vo:'); xo=input('Entre c/ xo:'); fid = fopen('TAB.DAT','w');
N = 0; while N < Nlin; t = to + N * dt; v = vo + a * t; x = xo + vo *t + 0.5* a*t^2; fprintf(fid,'%e %e %e \n',t,v,x);
N = N +1; end; fclose(fid); disp('FIM DA EXECUÇÃO')

2.2.

Execução do script

2.3.

Diretório de trabalho

2.4.

Verificação do número de linhas no arquivo TAB.DAT

Resultado: 1000 linhas.

2.
a)

Trabalho 1, parte (1.2)
Ao final do código (script) “tabela.m” adicione a linha load TAB.DAT e abaixo outra linha plot(TAB(:,1), TAB(:,3),’b–’); grid; salve, execute o código e descreva o que aconteceu.
%Código fonte do Exemplo 1 em MATLAB/OCTAVE
% tabela.m
% cria TAB.DAT c/ t v e x
Nlin = 1000; dt = 0.1; a =input('Entre c/ a :'); to=input('Entre c/ to:'); vo=input('Entre c/ vo:'); xo=input('Entre c/ xo:'); fid =