Precisión
Publicado por Sarita (5 intervenciones) el 17/10/2016 06:57:39
Hola. necesito ayuda con este programa. Lo que pasa que al correrlo, se me indetermina, entonces me comentaron que es un error de aproximación, pero no se como resolverlo. Seria de muy buena ayuda su apoyo. gracias.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
clc;
clear all;
%% CALCULO PARA EL CALOR ESPECIFICO DE UN NANOALAMBRE DE SILICIO %%%
hbar = ((6.62606957e-34)/(2.0*pi)); % J/s = 1.054571628(53)X e^-34 J/s
kb = 1.3806488e-23; % J/K
archivo = 'w_sinw.txt';
filas = 1;columnas = 201;A = fopen(archivo);
M = fscanf(A, '%f', [filas, columnas]); %Valores propios de cada punto
fclose(A);
N=67;for i=3*N
M(i)= 1.0e12*M(i);
end
fid = fopen('Calor_esp_sinw_T250.dat','w');
for jj=1:100
T=jj*1.0d0;
sum=0.0d0;
for j = 2:3*N
Q = (((M(j)*hbar)/(kb*T))*((M(j)*hbar)/(kb*T)))*exp((hbar*M(j))/kb*T);
S = (exp((hbar*M(j))/kb*T)-1.0d0)*(exp((hbar*M(j))/kb*T)-1.0d0);
F = Q/S;sum =sum+F
end
fprintf(fid,'%12.8f, %8.4e, \n',T,sum);
end
fclose(fid);
Valora esta pregunta
0
