summaryrefslogtreecommitdiff
path: root/1736/CH1/EX1.3/Ch01Ex3.sce
blob: 502ad6be2f0518bf94602a965f6fd0f7451a371e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
// Scilab Code Ex1.3 : Page-17 (2006)
clc; clear;
V = 10.58e-29;  // Volume of the unit cell, metre cube
a = poly(0, 'a');      // Declare a variable
a = roots(3*sqrt(3)/2*1.58*a^3-V);        // First lattice parameter, m
c = 1.58*a(3);     // Third lattice parameter, m
printf("\nThe lattice parameters of hcp structure of Ti are:");
printf("\na = %4.2f angstorm, c = %4.2f angstorm", a(3)/1e-010, c/1e-010);

// Result 
// The lattice parameters of hcp structure of Ti are:
// a = 2.95 angstorm, c = 4.67 angstorm